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The adaptation of Chree’s method of analyzing geomagnetic data, as carried out by Monk 
and Compton for their Mexican cosmic-ray data, has been applied to cosmic-ray intensity 
measurements at Boulder, Colorado. Secondary pulses in the cosmic-ray intensity are apparent 
in the difference curves. These occur at approximately 28-day intervals both before and after 
the primary pulse, in good agreement with the findings of Monk and Compton, but the Boulder 
pulses are somewhat smaller and more irregular. Within the four (previous and subsequent) 
pulse intervals investigated, magnitudes of the secondary pulses do not appear to depend upon 


the interval from the primary pulse. 


OLLOWING a preliminary application of 
Chree’s' method of analysis by Gill,? Monk 
and Compton’ made an extended investigation 
of cosmic-ray data from a Carnegie Institution 
Model C meter at Teoloyucan, Mexico. For the 
period of investigation, February, 1937, through 
December, 1938, they found intensity variations 
of considerable regularity at about 28-day in- 
tervals from the primary pulse. The secondary 
pulses, of about 0.3-percent amplitude, did not 
display the decrease in amplitude with increase 
of pulse number which characterized Chree’s 
magnetic character pulses. 

Having available cosmic-ray data obtained at 
Boulder by means of locally designed and con- 
structed equipment, it was thought desirable to 
investigate whether these displayed similar re- 


* Presented at the 1940 annual meeting of the American 
Physical Society in Philadelphia. 

1C. Chree, Phil. Trans. Roy. Soc. A212, 75 (1913); 
A213, 245 (1914); C. Chree and J. M. Stagg, Phil. Trans. 
Roy. Soc. A227, 21 (1928). 

?P. S. Gill, Phys. Rev. 55, 429 (1939). 

3A. T. Monk and A. H. Compton, Rev. Mod. Phys. 11, 
173 (1939). 


currences of intensity fluctuations. The observa- 
tions were made by Dr. V. A. Long and Mr. R. 
M. Whaley‘ for the purpose of investigating 
bursts, diurnal and seasonal variations in in- 
tensity, and various geophysical influences. Be- 
cause they will give a more detailed discussion 
of the apparatus when their findings are prepared 
for publication, only a brief description is in- 
cluded here. 

The writer’s’ 13.8-liter chamber with air at 
about 160 atmospheres was used. The 3.7-inch 
diameter sphere used in an earlier investigation® 
was incorporated as a symmetrical central elec- 
trode. Equipped with a 5-inch lead shield, the 
chamber was located in a room in the basement 
of the same building where the earlier visual 
observations were made. Guard tubes were made 
small. In the sealed compensation condenser (of 
the usual resistance-capacitance bridge arrange- 
ment for application of the collecting field) was 


a rs. ag Long and R. M. Whaley, Phys. Rev. 59, 470A 
1). 
5]. W. Broxon, Phys. Rev. 37, 1320 (1931). 
6 J. W. Broxon, Phys. Rev. 42, 321 (1932). 
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Difference Scale in Percent of General Average I, 


Days Following Selected High-and Low-I, Days; n 


Fic. 1. Primary and subsequent pulses in cosmic-ray 
intensity at Boulder, Colorado. The arrows show the 
largest of the probable errors of the means of the 63 to 
75 values for each of a dozen n-days selected at random. 


placed a capsule of radium salt which produced 
a compensation current slightly in excess of the 
average cosmic-ray current. The indicator was a 
quadrant electrometer automatically earthed for 
2 minutes at intervals of 1 hour. Deflections were 
recorded on a 20’’X8” photographic paper on a 
drum revolving once in 25 hours. Sensitivity was 
checked daily. The room temperature was regu- 
lated, and all the equipment except the photo- 
graphic apparatus was contained in an insulating 
box. The temperature of the air in the box 
varied about 1°C during the 18 months of 
observation. 

Data for the period from May 25, 1938, 
through December 1, 1939, were used. Deflec- 
tions for the individual hours, corrected for 
bursts, were supplied by Dr. Long and Mr. 
Whaley. These were arranged according to 
Greenwich days, daily averages found, and 
barometric coefficients determined for every 
month. Because these showed no regular seasonal 
variation, an average barometric coefficient was 
used in reducing to the average barometric 
pressure of 24.75 in. Hg. Corrections for the 
slight decay of the Ra were hardly necessary, 
but were applied. Application of the instrument 
constants indicated an average cosmic-ray ioniza- 
tion current of 38.16 ions per cc per sec. 

Unfortunately, complete records were not ob- 
tained for all days. On this account, 14 days of 


the period were regarded as blank because 
records for less than 8 hours per day were avail- 
able. However, 67 imperfect daily records were 
retained. Of these, 5 had records for 8 but less 
than 12 hr., 16 had records for 12 but less than 
18 hr., 27 had records for 18 but less than 23 hr., 
and 19 had records for 23 hr. In view of the 
method of analysis, it is significant that the 14 
recordless days among the 556 were distributed 
in such a manner that not more than 4 occurred 
in any single calendar month. 

The procedure of Monk and Compton was 
followed closely. The five days with greatest 
cosmic-ray intensity in each month from June, 
1938, through November, 1939, were selected as 
the zero-days for the positive-pulse curves. The 
five with least intensity in each month served as 
zero-days for the negative-pulse curves. For any 
particular zero-day, the subsequent days were 
assigned positive numbers and the previous 
days, negative numbers, in the order of their 
occurrence after or before the corresponding 
zero-day. 

Tables of cosmic-ray ionizations were then 
formed with a column for each of the zero-days 
during the first 15 months and for each corre- 
sponding day with positive number to »=135. 
Similar tables were formed for the zero-days of 
the last 15 months and the corresponding days 
with negative numbers to m= — 135. Full columns 
had 75 items; none had less than 63. 


Difference Scale in Percent of General Average I, 


Doys Preceding Selected High-and Low-1, Days 


Fic. 2. Primary and previous pulses in cosmic-ray 
intensity at Boulder, Colorado. 
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RECURRENCE PHENOMENA IN COSMIC-RAYS 


The averages of the values in the several 
columns were plotted as ordinates with the 
corresponding day numbers as abscissae to obtain 
the upper pairs of curves of Fig. 1 and Fig. 2. 
The full-line ‘positive’ curve in each case has as 
its zero-days the five per month of greatest 
ionization. The broken-line ‘‘negative’’ curves 
have days of least ionization as zero-days. In 
Fig. 1 are the subsequent-pulse curves corre- 
sponding to positive day numbers, while Fig. 2 
contains the previous-pulse curves corresponding 
to negative day numbers. The lowest curve in 
each diagram, labeled ‘‘difference,” has ordinates 
obtained by subtracting from the ordinate of 
the ‘‘positive”’ curve for a particular day number, 
the ordinate of the ‘‘negative’’ curve for the 
same day number. 

The “primary pulse’ variations corresponding 
to zero-days for the entire 18 months (not shown 
in the diagrams) were about 0.7 percent of the 
average for the period. The corresponding pri- 
mary pulse magnitude was about one percent for 
the Mexican data. The secondary pulses shown 
in both the positive and negative curves for both 
positive and negative day numbers are quite 
irregular, and the season-like variation is con- 
siderable. The unidirectional trend is canceled 
out in the difference curves, however, and rather 
distinct though quite irregular pulses are ap- 
parent. These appear to be centered fairly well 
at the 28-day intervals designated by the vertical 
lines. It should be borne in mind that the 
deviations or “‘pulses’’ represented in the differ- 
ence curves must be regarded as of double the 
magnitude of the “‘pulses”’ in the original curves. 
On this basis the ‘amplitude’ or magnitude of 
the pulse deviations in the original curves must 
be regarded as rather less than 0.2 percent of 
the mean. Neither subsequent nor previous 
pulses show any regular tendency to decrease 
with distance from the primary pulse. In this 
respect they confirm the observations of Monk 
and Compton who have found recently’ that the 
pulses in the intensity at Teoloyucan persist for 
ten 27- or 28-day intervals without decrease in 
amplitude. (Professor Compton informed the 
writer that when extended to ten intervals the 


7A. H. Compton and A. T. Monk, Phys. Rev. 59, 112 
(1941). 


Difference Scalé in Percent of General Averoge I, 
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Fic. 3. Combination of difference pulses of Fig. 1 and 
Fig. 2. Ordinates are averages of differences for the same 
day number, regardless of sign. 


persistence of the pulse amplitude was shown 
clearly, although the earlier extension to eight 
pulses had indicated a decrease of amplitude.) 

In Fig. 3 is shown the combination difference 
curve. The ordinates of this curve represent the 
average of the ordinates for positive and negative 
day numbers of the same magnitude. Irregu- 
larities are not removed in this curve as in the 
corresponding curve by Monk and Compton. 
Although quite as irregular as in the individual 
difference curves, the pulses in this combination 
curve are of about the same magnitude and 
frequency as in the difference curves of Fig. 1 
and Fig. 2. Perhaps it is quite as significant that 
they are not less regular or less definite. 

The smaller magnitude and the greater irregu- 
larity of the pulses shown here in comparison 
with those found by Monk and Compton may 


‘be due in part to the differences in geographical 


location and in part to the fact that a shorter 
period was considered in the present instance, 
but they may also be due in part to differences in 
meteorological and other conditions. In both 
cases corrections were made for barometric 
changes, but none were applied for temperature 
or other effects. The Boulder data display a 
pronounced long time outdoor temperature effect, 
and a somewhat smaller short time temperature 
effect, as well as a magnetic effect. While the 
temperature at Teoloyucan is very uniform (and 
the temperature coefficient is small according to 
Gill?) the monthly average temperature at Boulder 
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varied about 50°F. It should also be mentioned 
that no ‘“‘absolute”’ calibration of the apparatus 
by comparison with one whose residual ionization 
has been measured under an enormous shield 
has ever been made. However, the cumulative 
observations made with the apparatus and the 
magnitudes of the several “‘coefficients’’ obtained 
from the data it yields, have satisfied the writer 
that the radioactive contamination of the cham- 
ber must be very small. 

While Monk and Compton consider the pulses 
likely due to a solar influence, they consider that 
the persistence of the pulse amplitude of the 


cosmic-ray fluctuations and the decrease in 
Chree’s magnetic pulses with distance from the 
primary pulse indicate the two varieties of dis- 
turbance are of different origin. However, the 
similarities of the frequencies of the pulses ob- 
tained by the same method of analysis combined 
with Graziadei’s* association of his approximately 
0.4-percent cosmic-ray intensity fluctuations of 
27.2-day periodicity with solar disturbances, 
indicate the desirability of further examination 
of this point. 


8H. T. Graziadei, Akad. Wiss. Wien. [Ila] 145, 495 
(1936). 
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The scattering of particles of spin $ and mesotronic mass m by heavy nuclear constituents is 
investigated with the pair theory of Critchfield and Teller, in which a nucleus appears as an 
extended source of mesotrons. A rigorous formula, (19), is found for the cross section, by 
obtaining explicitly the normal coordinates of the perturbed mesotron field. When the ‘‘con- 
stants” are adjusted to give a model of nuclear forces with approximately the right range and 
saturation properties, o—°32(mc?/E)*(h/mc)?. Although, in other models, ¢ may decrease more 
rapidly with the energy E, it is always about 10°‘ cm? at low energies, as contrasted with the 
cross section of about 10~*8 cm? observed for cosmic-ray mesotrons. These results are discussed 
in relation to other accounts of the pair field theory. 


RITCHFIELD, Teller, Wigner and Lamb'~* 
have developed a field theory which ex- 
plains, in a qualitative way, the saturation and 
spin dependence of nuclear forces. The field, 
when unperturbed by nuclei, is assumed to be a 
Fermi gas of charged particles with the mass of 
the mesotron and described by the Dirac equa- 
tion. A proton or neutron interacts with this 
field by emission and absorption of neutral 
pairs of ‘“‘mesotrons,’’ and such processes can 
occur not only at the point which specifies the 
position of the nucleus, but over a finite region, 


(1937). L. Critchfield and E. Teller, Phys. Rev. 51, 289 
(aos = Critchfield and Teller, Phys. Rev. 56, 530 
3C. L. Critchfield, Phys. Rev. 56, 540 (1939). 
*C. L. Critchfield and W. E. Lamb, Jr., Phys. Rev. 58, 
46 (1940). 


spherically symmetrical about that point, and of, 
roughly, the nuclear radius. 

-Past discussions have been concerned with 
bound states of the pair field and with nuclear 
forces; but here we shall find how mesotrons 
are scattered by nuclei. It is of some interest to 
compare our conclusions with the remarkably 
small experimental cross section for the scatter- 
ing of cosmic-ray mesotrons by nuclei: for meso- 
tron energies of =350 Mev, Wilson®> has ob- 
tained an upper limit of 10-** cm? for the cross 
section. Our results seem to be about a hundred 
times larger.® 

. G. Wilson, Proc. Roy. Soc. 174, 73 (1940). 

E. Marshak and V. F. Weisskopf [Phys. Rev. 59, 
130 (1941), have found agreement with experiment witha 
similar pair-field theory, which, however, pictures the 


nucleus as a point source. The divergent nuclear interaction 
characteristic of such a model is fitted to experiment by 
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PAIR THEORY OF 


Although it is not certain that mesotrons of 
half-integer spin are present in cosmic rays or 
that they intervene in 8-processes, it is likely, 
indeed, that a nucleus is better described as an 
extended source of charged particles, rather 
than a point source. In theories embodying such 
a description, the formal methods of this note 
may find further application. 


I 


The notation used is that of Lamb and Critch- 
field;* the coupling constant is the nuclear 
source distribution is given by u(x) satisfying 


fucyax=t; 


and the total field Hamiltonian may be ex- 
pressed in terms of the quantized amplitude 
y(x), a 4-rowed one-column matrix in the spin- 
variables: 


(1) 


f Vt(x)[8—i(@, ¥) W(x)dx 


f f | (2) 


8, w are the usual Dirac matrices, h=y=c=1 in 
these units; and the commutation laws charac- 
teristic of the Fermi statistics are assumed to 
hold among the components of y and yf. 

By resolving ¥(x) into its normal vibrations 
We(x) we may transform H to diagonal form. 
The equation governing Wz(x) is self-adjoint ; 


[E+i(a, ¥)—BWe(x) 
+nu(x)B f u(x’We(x’)dx’=0; (3) 


and, therefore, the spectrum of E is real, and 
the solutions of yz are orthogonal in the sense: 


(E-E’) f (4) 


“cutting off” the nuclear forces. The smallness of the 
coupling constant obtained in this way is chiefly responsible 
for their small cross section. We wish to take this oppor- 
tunity to thank Drs. Marshak and Weisskopf for their 
kindness in sending us their manuscript before publication. 
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777 
Dividing (3) by the operator [E+i(a, ¥)—8 ], 


nu(x) 


¥e(x) +——. 
E+1(a, v) 


B f 
Yo(x) ’ (5) 


where Yo(x) must satisfy the Dirac equation and 
the boundary conditions appropriate to a wave 
function, 


[E+i(a, ¥)—8 =0. (6) 
(5) may be rewritten in the form 
=Yo(x), (7 
where (7) 
(8) 


By multiplication with n8u(x) and integration 
over all x, 


[1+n(1-+BE) F(q) f u(x)po(x)dx. (9) 


In this equation, 
F(q)= f u(x)[A+q? }-'u(x)dx 
(10) 
0 


where v(p)=v(p) is the Fourier component of 
u(x); i.e., 


v(p) = (27) f dxu(x) exp[ —i(p, x) ], 
. (11) 
arf v(p)*p*dp=1. 


For every value of E, such that gq is real, | | 21, 
¥o(x) is a free-particle solution of the Dirac 
equation with propagation vector q, g’=(q, q): 


¥o(x) = ¢(q) exp[7(q, x) J, (12) 


where ¢(q) is a matrix depending on q and the 
sign of E. If |E| <1, on the other hand, yo 
vanishes identically; and then (9) becomes a 
secular equation for E. Normal vibrations corre- 
sponding to the discrete spectrum of solutions 
are quadratically integrable and cannot con- 
tribute to the scattering problem. 
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In (10), therefore, g must be taken real. The 
integrand of F(q) has a pole at p=q; and F(q) 
is defined by deforming the contour into the 
lower half-plane. This makes u(x)/(A+q’) well- 
defined and asymptotic to an outgoing wave, a 
condition necessary for scattering : 


F(q) = f(q) — 2n*iqu(q)?, 
(q) = f(q) — 13) 


f(q)=4 dplp*v(p)? — ]/(q?— 


A useful form of these relations is obtained by 
direct application of the Green’s function for 
outgoing waves, and explicit integration over 
angles: 


faxfax 


4r|x—x’| 


x (14) 


The real part of F(qg) is given as : 


+o 
=+4r] d dx'xx’ 
f@ J x xx’u(x) 
Xu(x’)(2g)~! sin(g|x—x’|). (15) 
We may now solve (9) for {zg and (7) for Wz: 


_ (27)!nv(q) 
1+2nF—(nFq)* 
ve(x) —$(q) exp[i(q, x) ¥) 


—(E—B)F(q) ]é(q), (16) 


, 


x x |] 
~[E+B-—i(a, ¥) 
x f exp -ia(a, x)] 
~ 0) (17) 


where n=x/x. 

The expectation value of the particle-density 
operator in a stationary state of the field is a 
sum of terms of the form Wzf(x)Wze(x) over all 


occupied normal vibrations; and therefore each 
normal vibration has its own scattering coeffi- 
cient per unit solid angle: 


1+2nF—(ngF)?|? 
X [8 —n(E—8) F*(g) J[E+6+¢(a, n) 
x [8—n(E—8) F(q) ]é(q). (18) 


st(q) 


We may set 
=exp[ —8(a, q)(2g)~! 


and average over the spin orientation of £(0) to 
obtain a result, independent of the sign of E: 


|1+2nF—(ngF)?\? 
cos*3d], (19) 


where cos? =(q, x)/gx. The cross section be- 
comes independent of 7 as n— «: 


167‘v(g)* cos?}0 
limog = (20) 


The functions which appear in the foregoing 
analysis may be discussed generally in terms of a 
mean width a of the source-distribution func- 
tion u: 


Thus, 

f(0)=—-a*, f(a")=0, fq)~q* 
for (22) 

and 


[gv(q)?]=0 for g=0, av(a~)*=max., 
(23) 


for g>a-'. In order that a be approximately the 
nuclear radius it should be taken 1 in our 
units; and one may observe that og is measured 
in units of 5.13 cm’. 

With a particular choice of u(x), it is not 
difficult to obtain explicit formulae and numerical 
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19) 


PAIR THEORY OF 


values for the cross section; e.g. for 
u(x) 


= 


(24) 
f(q) = 
F(q) = —a*(1—iaqg)~*, 
limog = cos*30. (25) 


A more appropriate choice of u might be 
u(x) = exp(—x*/a’), 
v(p) = exp(—a*p?/4), (26) 
= Fi(1, 2, —a°g?/2); 
and 
limog = F1(—}, 3, a?g?/2)? 
+7a*gq?/2}-' cos?}3. (27) 
The experimental values, which make ESS, 


can be inserted, to obtain the total cross sections 
for the respective functions (24) and (26): 


5X 10-4 cm? >a 
cm? cm’, 


the larger value corresponding to E=1. 

In the second case, the cross section is too 
large for ES4, a large fraction of the experi- 
mental range in these experiments. 

We conclude that this theory does not account 
for the small scattering of cosmic-ray mesotrons. 


II 


The equation (19) for oz shows what is 
apparently resonance whenever 


Real (28) 
and when this condition obtains, 
Loz cosd 


which, when integrated over all angles gives the 
expected [Gr ]res=4aq~. To satisfy (28) one 
must take 


na*=} for for g>a"; (29) 


but for g<a~' one would need imaginary values 
of ». The apparent resonance, in both cases, 
actually makes the value of og very large com- 
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pared to that when g is displaced by an amount 

We may understand the origin of the low 
energy resonance by studying the discrete spec- 
trum of our fundamental equation (3), which 
has been defined by ¥o=0, | E! <1, g=iQ in (9): 


[ —4rn(1 +88) =0. 


Thus ¢¢ must be an eigenvector of 8, satisfying 
B&e=+fe, which gives rise to the secular equa- 
tion for Q: 


f dpp*v(p)?/ 
(p?+Q2)=+1. (30) 


The alternative sign results in there being al- 
ways but one solution for Q as a function of », 
and two solutions for E=+|(1—@Q*)!|. (In 
the particular cases of (24) and (26) the secular 
equation has the respective forms: 


F, 
and 


na*[1+(1—Q*)! 
with 


aQ 
’=(aQ) = f exp(—/°/2)dt.) 


As 7 is increased from zero, a solution of the 
secular equation first becomes possible when 
Q=0, E=+1 and na*=} according to (30) 
and (21). This means that there exist two dis- 
crete levels, each doubly degenerate because of 
the spin, one of which separates itself from the 
continuum at E=+1, the other at E=—1. 
As 7 increases, they approach, cross at E=0; 
and then each moves asymptotically toward the 
point of origin of the other. The condition for 
low energy resonance in gg is precisely that for 
the appearance of discrete levels at the edges of 
the continuum. 

The pair-field theory of nuclear forces would 
require » to be quite large compared to a™, 
however, and the foregoing resonance could 
never occur. Instead, the high energy resonance 
for gn would be present; and, for this, there 
seems to be no familiar physical’ picture. Its 
presence appears to be connected with the exist- 
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ence of the bound, lowest, state of the pair- 
field considered in references (2), (4). 

To see how this comes about, and also to 
bring our account of the behavior of the pair- 
field into closer connection with that given in 
references (2) and (4), we present a_phase- 
shift analysis of the scattering process. 

In polar coordinates, following Dirac,’ we 
may rewrite (3): 


[E+ 
+nu(x)B f (31) 


where j=6[(e, L)+1], the ‘angular momen- 
tum”’ operator, «=(a, x)/x, L= —i[xxV], and 
Wi=k). with k any nonvanishing integer. 
The perturbation term affects only those func- 
tions with spherically symmetric parts, i.e., 
those with k=+1, because of the spherical 
symmetry of u(x). For all values of |k| >1, the 
integral vanishes, and we have to deal with the 
equation for a free particle, the solutions of 
which show no scattering. In what follows, 
therefore, we shall confine ourselves to k= +1. 

The angular dependence of ¥;,(x) now occurs 
only in the form of the first spherical harmonic; 
and we may eliminate the angular coordinates 
through the formal introduction of ¢: 


=[1 — (te/g*)(E—B) (d/dx) Sx(x)/x], 


where S;,(x) is a scalar function of the radius 
alone, and ¢; is a 4-rowed, one-column matrix 
satisfying the characteristic equation Bf, =k¢;,. 
There results for S;(x), the relation: 


(d/dx)*}S;(x) +42n(1+kE) 
This is quite comparable to (3) and may be 


solved in similar fashion with the aid of the 
Green's function 


+0 
f sing |x—x’|. (33) 


The solutions of (3) and of (31) are brought into 


™P. A. M. Dirac, Principles of Quantum Mechanics 
(Oxford University Press, 1935), §73. 


agreement by observing that 
0 


in accordance with (33) and (15). In this way, 
one obtains: 


2n*qu(q)?n(1+kE) 
and the phase shift 6, is defined by 
1+n(1+kE) f(g) 


When 6,/z is given the customary interpreta- 
tion of the fractional momentum shift of a 
vibration in the continuum, due to the perturba- 
tion, this becomes exactly Eq. (9) of reference 
4;5 and it leads directly to the value obtained 
therein, for the energy of the bound state of the 


field : 
AE=—(2/n) f (36) 
0 


osgx, (34) 


S).(x) ~ singx — 


(35) 


tané,=— 


(This is, of course, not quite the Jowest energy 
of the field; for one could also fill the two dis- 
crete states of energy ——1. Their contribution 
to AE would be small, however, when 7 is taken 

For values of E<0, with g>a™", 64; remains 
but for g<n and for g>n, 
varying through the resonance value of 2/2 
in the neighborhood of g=n. The resonance of 
the cross section at g=7n can thus be connected 
with the limitation of AE to finite value, of 
order —7n. 

The problem of scattering of mesotrons when 
the coupling between the nucleus and the pair- 
field is spin dependent, cannot be solved in the 
same manner as the problem treated here, be- 
cause the reaction of emission and absorption pro- 
cesses on the nuclear spin plays an essential part. 

I wish to thank Professor J. R. Oppenheimer 
for suggesting this problem and for many 
stimulating discussions concerning its solution. 
It is a pleasure, too, to recall Dr. Critchfield’s 
helpful and very cordial letters. 


8’ Dr. Critchfield has independently arrived at the 
formula (19) by applying this interpretation to his Eq. (9), 
reference 4. With his friendly agreement, this note appears 
in its original form. 
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The excitation curves for the emission of y-rays, neu- 
trons, and protons from the disintegration of C® by deu- 
terons were investigated. Resonances for y-ray emission 
were found at 0.92, 1.16, 1.30, 1.43 and 1.74 Mev. The 
resonance at 1.43 Mev has a half-width less than 10 kev 
while the other resonances are considerably broader. 
Resonances for neutron emission were found at 0.92, 1.16, 
1.30, 1.74 and 1.82 Mev. Resonances for proton emission 
are found at 0.92, 1.16, 1.23 and 1.74 Mev. The y-rays 
from C® were shown to come from the reaction C"(d, p)C'*. 
The short range protons from this reaction were observed 


and had a range of 1.1 cm. The Q value for this 1.1-cm 
proton group is —0.52+0.07 Mev. The excitation curves 
for the emission of y-rays and protons from the disintegra- 
tion of C® by deuterons were investigated. Both curves 
showed a resonance at 1.55-Mev deuteron energy. The 
y-rays from C™ have an energy of 5.5 Mev and come from 
the reaction C"(d, n)N“*. The low energy group of neu- 
trons from this reaction were found when a 23-percent C™ 
target was used. The Q value for this group of neutrons is 
0.40 Mev. The relative and absolute yields of the various 
reactions were determined. 


HE protons from the disintegration of 

carbon by deuterons were observed by 
Cockcroft and Walton! in 1934. When they 
bombarded carbon with 0.5-Mev deuterons they 
observed protons with a range of 14 cm which 
they attributed to the reaction: 


(1) 


This reaction has been extensively studied since 
that time and the value of Q; is now accurately 
known to be 2.71+0.05 Mev.’ 

Radioactive nitrogen from the bombardment 
of carbon by deuterons was simultaneously ob- 
served by Crane and Lauritsen’ and by Hender- 
son, Livingston and Lawrence.‘ The production 
of radioactive nitrogen was attributed to the 
reaction : 


(2) 


Tuve and Hafstad® found that the neutrons pro- 
duced in this reaction were of low energy. Later 
Bonner and Brubaker® obtained an accurate 


* Now at Bartol Research Foundation. 
1J. D. Cockcroft and E. T. S. Walton, Proc. Roy. Soc. 
144, 704 (1934). 
2]. D. Cockcroft and W. B. Lewis, Proc. Roy. Soc. 154, 
261 (1936). 
1958). R. Crane and C. C. Lauritsen, Phys. Rev. 45, 430 
* M. C. Henderson, M. - Livingston and E. O. Lawrence, 
Phys. Rev. 45, 428 (1934 
ag 35). A. Tuve and L. R. Hafstad, Phys. Rev. 48, 106 
*T. W. Bonner and W. M. Brubaker, Phys. Rev. 50, 
308 (1936). 


measurement of the range of the recoil protons 
produced by these neutrons in a cloud chamber. 
From the range-energy curve for protons which 
was available at that time a value of Q2= —0.37 
Mev was obtained. Later, when a more accurate 
range-energy curve for protons was obtained 
from the experiments of Parkinson, Herb, Bel- 
lamy and Hudson,’ a recalculation’ from these 
data gave a value of Q2=—0.25+0.03 Mev. 
Another value of Q2 was obtained by Lewis and 
Burcham$ from the observation that the forma- 
tion of N® begins at 0.32 Mev. This gives a 
value of Q2= —0.28 Mev which is in good agree- 
ment with the value of Bonner and Brubaker. 
Bonner and Brubaker also found a neutron 
group representing about 1 percent of all the 
neutrons from carbon for which a Q value of 
5.2 Mev was determined. They also observed 
another weak group of neutrons corresponding 
to a Q value of 1.2 Mev. The production of these 
two high energy groups of neutrons was at- 
tributed to the C"™ isotope according to reaction: 


(3) 


Gamma-rays from the bombardment of carbon 
by 0.8-Mev deuterons were first reported by 
Lauritsen and Crane.® From absorption experi- 


7D. B. Parkinson, R. G. Herb, J. C. Bellamy and C. M. 
Hudson, Phys. Rev. 52, 75 (1937). 

tT, W. Bonner, Phys. Rev. 53, 496 (193 
< E. Bretscher, Kernphysik (J. Springer, Berlin, 1937), p 

°C. C. Lauritsen and H. R. Crane, Phys. Rev. 45, 345 
(1934). 
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Fic. 1. Diagram of the observational end of the 
Rice pressure Van de Graaff generator. 


ments they estimated that rays of 3.5 Mev 
were produced. The y-rays were later studied by 
McMillan’ who carried out absorption measure- 
ments in various elements. From these measure- 
ments he deduced a y-ray energy of 3 Mev. 
Other experiments on the y-rays from carbon 
bombarded by deuterons were carried out by 
Tuve and Hafstad."! From cloud-chamber studies 
of the Compton electrons, they obtained evi- 
dence for a strong line at about 2.7 Mev and a 
much weaker one extending up to 4 Mev or 
higher. It was at first thought that the y-rays 
from carbon might accompany the emission of 
the 14-cm protons of reaction (1). At this time 
the mass of C" was considerably in error and 
the energy of the y-ray had to be added to the 
energy of the 14-cm protons to give a Q value 
consistent with the isotopic masses. Later, in 
calculating the masses from disintegration data, 
Bethe found” that a consistent set of masses 
was obtained only if the y-ray energy was not 
added to the proton energy. When the masses 
of C® and C" were accurately determined by 
Bainbridge and Jordan," it was certain that the 
energy of the 14-cm protons alone gave the 
proper value of Qi. 

The weak y-rays of 4 Mev (or higher) were 
attributed by Bonner and Brubaker to reaction 
(3) involving C™. There is sufficient energy 
available in this reaction for these y-rays to 
come from an excited N"™ after the emission of a 
low energy group of neutrons. 


1 E. McMillan, Phys. Rev. 46, 868 (1934). 

1935), A. Tuve and L. R. Hafstad, Phys. Rev. 48, 106 
( 

2H. A. Bethe, Phys. Rev. 47, 633 (1935). 

18K, T. Bainbridge and E. B. Jordan, Phys. Rev. 51, 384 


(1937). 
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The discovery of a weak group of 48-cm 
protons from carbon bombarded by 0.8-Mev 
deuterons was made by Bower and Burcham." 
The protons come from the reaction: 


C834 (4) 


This reaction was also independently observed by 
Pollard" at higher bombarding voltage. 

Recently we have reported in a Letter to the 
Editor in The Physical Review that the excitation 
curve for the emission of y-rays from carbon 
bombarded by deuterons shows resonances.'* 
Resonances were reported at 0.92, 1.16, 1.30, 
1.43 and 1.74 Mev. Furthermore it was found!’ 
that the neutrons from reaction (2) showed some 
of these same resonances; namely, those at 0.92, 
1.16, and 1.30 Mev and an additional resonance 
at 1.80 Mev. Further experiments'*® showed that 
the 14-cm protons from reaction (1) also showed 
the resonance at 0.92 Mev. A later Letter'® gave 
results on the yield of the high energy protons 
from C, The present paper is a report of ex- 
tended results on the yield of the products from 
carbon bombarded by deuterons. 


YIELD 


PROTON ENERGY IN KV 


Fic. 2. Gamma-ray resonances from F!*+H! which were 
used in calibrating the electrostatic voltmeter of the Van 
de Graaff generator. The observed 8.5-kev half-width is 
thought to be chiefly experimental and represents the mini- 
mum width which could be obtained for other resonances. 


. C. Bower and W. E. Burcham, Proc. Roy. Soc. 
A173, 379 (1939). 
% E. Pollard, Phys. Rev. 56, 1168 (1939). 
(1940). E. Bennett and T. W. Bonner, Phys. Rev. 58, 183 
17 T. W. Bonner, E. Hudspeth and W. E. Bennett, Phys. 
Rev. 58, 185 (1940). 
18M. M. Rogers, W. E. Bennett, T. W. Bonner and E. 
Hudspeth, Phys. Rev. 58, 186 (1940). 
19 W. E. Bennett, T. W. Bonner, E. Hudspeth and B. E. 
Watt, Phys. Rev. 58, 478 (1940). 
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EXPERIMENTAL ARRANGEMENT 


The Rice Institute pressure Van de Graaff 
machine was used in the experiments on the 
disintegration of carbon by deuterons. A com- 
plete description of the apparatus will be pub- 
lished elsewhere. For a description of the present 
experiment, it is sufficient to say that the 
machine is of the same general type as the 
pressure machine developed by Herb and col- 
laborators”® at the University of Wisconsin. 

A schematic diagram of the observational end 
of the high voltage apparatus is given in Fig. 1. 
The ion beam passes through a gate valve after 
emerging from the accelerating tube. Then the 
beam passes through a box which contains an elec- 
trically operated shutter. After passing through 
the shutter box the beam is bent through 90° 
by a large electromagnet, which resolves the 
beam into different components, corresponding 
to the different values of e/m for the positive 


YIELD 


1370 1410 1450 


DEUTERON ENERGY IN KV 


Fic. 3. The y-ray resonance at 1.43 Mev obtained by 
using a carbon target 3} as thick as the thin one used for the 
excitation curve of reference 16. While this shows that 
this 10-kev half-width is not due to target thickness, 
comparison with the half-width obtained in Fig. 2 for the 
fluorine y-ray resonances indicates that this may really be 
a much narrower excitation level. 


ions. The current then enters a Faraday cage 
before falling on the target. The ion current on 
the target is measured by a current integrator of 
the neon discharge type.** The voltage at the 
high potential center electrode is measured by 
means of an electrostatic voltmeter. This volt- 
meter was calibrated by observing the resonances 


20R. G. Herb, D. B. Parkinson and D. W. Kerst, Phys. 


Rev. 51, 75 (1937). 
1B. E. Watt, Rev. Sci. Inst. in press. 
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Fic. 4. Sample data of the detailed study which was made 
of the y-ray resonances of the curve of reference 16. 


for the emission of y-rays from the bombardment 
of fluorine by protons. The y-rays were detected 
by a Wulf type electroscope filled with argon at 
a pressure of 70 atmospheres. Figure 2 gives the 
results that were obtained when a thin target of 
CaF, was bombarded with protons. The reso- 
nances at 0.862 and 0.927 Mev were used to 
obtain the value of the constant of the electro- 
static voltmeter. The resonances” at 0.660 and 
1.363 Mev were also observed and each resonance 
voltage gave very nearly the same voltmeter 
constant. The subsequent readings of the volt- 
meter were thought to be reliable to at least 
20 kev. In order to stabilize the voltage and to 
facilitate making small changes in the voltage, 
a corona current of from 50 to 100 microamperes 
from the central electrode to an adjustable 
probe was used. This current passed from the 
probe to ground through a variable high re- 
sistance balanced against the e.m.f. of a battery. 
Small changes in the voltage could be noted by 
observing the change in the corona current from 
the high potential electrode. Over a region of 
about 0.10 Mev the change of corona current 
was found to be proportional to the change in 
voltage, and so could be used to measure the 
change in voltage. 


2 E. J. Bernet, R. G. Herb and D. B. Parkinson, Phys. 
Rev. 54, 398 (1938). 
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Fic. 5. Absorption curves taken with coincidence Geiger counters for the Compton electrons produced by the y-rays 
from C+H?. Curve (A) is for an ordinary carbon target. Curve (B) is for a 23-percent C™ target. Curve (C) is the same 
as curve (B) except that the ordinates have been multiplied by 5 to aid in determining the end point. 


y-RAays FROM C+H? 


The excitation curves for the y-rays have been 
obtained by using thin carbon targets which 
were usually made by the evaporation of paraffin 
onto silver disks. The thickness of most of the 
targets was only a few thousand volts. The 
y-rays were detected by means of coincidences 
of Geiger counters and by measurements with 
the Wulf type electroscope filled with argon at 70 
atmospheres. Since the electroscope gave a 
greater accuracy for a given time of observation 
than the counters, it was used in most of the 
work. Figure 1 of reference 16 gives the excitation 
curve for the y-rays from 0.55 up to 2.0 Mev. 
This curve shows resonance at 0.92, 1.16, 1.30, 
1.43 and 1.74 Mev. The experimental half- 
widths of the resonances vary from about 0.25 
Mev for the resonance at 1.74 Mev to 10 kev 
for the resonance at 1.43 Mev. The narrow 
resonance at 1.43 Mev was investigated further 
by using a target } as thick as that used to get 


the curve of reference 16. The results obtained 
with this thinner target are given in Fig. 3. 
The half-width of the resonance remains 10 kev, 
and so it is shown that the target thickness is 
not responsible for the experimental! width of the 
resonance. The two possibilities that remain are 
that the 10-kev width is due to the spread in 
the energy of the deuterons, or that the real 
width of the resonance level is 10 kev. Since 
this is about the same width as we obtained for 
the resonances from F+H!, we are inclined to 
think that at least most of the 10-kev experi- 
mental width is due to the spread in the energy 
of the deuterons, and so the real width of the 
level might be considerably less. 

We have studied the excitation curve for 
y-rays in more detail than is shown in reference 
16. Figure 4 gives a sample set of such data 
which were taken in the energy interval 1.24 to 
1.36 Mev. This curve shows that the level at 
1.30 Mev has a half-width of about 40 kev. 
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It is difficult to be certain how the y-rays 
originate without information about their energy. 
Consequently a series of absorption experiments 
on the y-rays were made at each of the reso- 
nances. We found that the absorption coefficient 
in Pb was near the theoretical minimum value 
and that the y-rays from each of the resonances 
had the same absorption coefficient. We have 
also measured the energy of the y-rays at each 
of the resonances by determining the range of 
the Compton electrons produced by the y-rays. 
We used a similar technique to that described 
in the following paper.** Curve (A) of Fig. 5 
gives the average results of the coincident 
Geiger-counter measurements on the y-rays from 
an ordinary carbon target. The results for each 
resonance agreed within the probable error with 
this average curve. The absorber thickness at 
which the number of coincident counts falls to 
half-intensity indicates a y-ray energy of roughly 
3 Mev, but the actual end point of the curve 
indicates an energy of about 5.1 Mev. This 
suggests that there is a strong y-ray with an 


%W.E. Bennett, T. W. Bonner and B. E. Watt, Phys. 
Rev. 59, 793 (1941) (this issue). 


energy of 3.0 Mev and a weak y-ray with an 
energy of 5.1 Mev. Similar results were obtained 
from a cloud-chamber study of these y-rays by 
Bonner, Becker, Rubin and Streib™ in Pasadena. 
When a potential of 1.0 Mev was used they found 
a y-ray with an energy of 3.0+0.2 Mev and 
another y-ray with approximately 5 percent 
relative intensity near 5.5 Mev. Using separated 
isotopes they found that the 5.5-Mev y-ray 
comes from C'. We have also carried out some 
experiments with 23-percent C™ which has been 
kindly supplied to us by Professor H. C. Urey. 
Curve (B) of Fig. 5 shows the absorption in 
aluminum of the secondary electrons produced 
by the y-rays. A thick carbon target was used 
for this experiment and the bombarding deu- 
terons had an energy of 1.5 Mev. It is apparent 
from curves (A) and (B) that the average 
quantum energy of the y-rays is considerably 
higher for the 23-percent C™ target. This con- 
firms the results quoted above that the higher 
energy component of the y-radiation is from C", 
and the curve gives a more accurate value of 


*T. W. Bonner, R. A. Becker, S. Rubin and J. F. 
Streib, Phys. Rev. 59, 215A (1941). 
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the maximum range of the Compton electrons 
than we obtained before. This maximum range 
is 10.3+0.3 mm of aluminum. From the range- 
energy relation for such electrons,” we find that 
this corresponds to a quantum energy of 5.5 
+0.2 Mev. This is in good agreement with the 
cloud-chamber results of Bonner, Becker, Rubin 
and Streib.* 

Since the y-ray excitation curves of Figs. 3, 4, 
and of reference 16 were obtained using ordinary 
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Fic. 7. Energy distribution of the neutrons from C+H? 
as inferred from recoil protons in a methane-filled cloud 
chamber. A target containing 23-percent C™ was bom- 
barded with 1.4-Mev deuterons. 


carbon targets, it is apparent that about 95 
percent of the y-rays recorded there were from 
the bombardment of C™. We have carried out a 
separate experiment to get the excitation curve 
for the y-rays from C. We were able to count 
y-rays from C® by using coincident Geiger 
counters with sufficient absorber between them 
to prevent 3-Mev rays from producing coincident 
counts. A target of 2-percent C' was used in 
these experiments. It was prepared by one of 
us (B.E.W.) from methane in which C" had 
been concentrated by thermal diffusion. Its 
thickness was found to be about 30 kev from 


the observation that the resonance at 1.43 Mev 
showed a width of 30 kev. Figure 6 shows the 
excitation curve for y-rays with an energy of 
5.5 Mev. The curve is smooth except for a 
resonance at 1.55 Mev which will be discussed 
in detail later in this paper. 

Since it seemed possible that the 1.43-Mev 
resonance for y-rays was of a different character 
from the other broader levels, we thought that 
the y-rays from that level might have a different 
quantum energy. We carefully searched for a 
sharp resonance in the yield of high energy 
y-rays at 1.43 Mev but found no effect. To 
further investigate the energy of the y-rays 
from the 1.43-Mev resonance, we put sufficient 
absorber between two coincidence counters to 
prevent 2.5-Mev y-rays from producing coin- 
cident counts. We then investigated the reso- 
nance at 1.43 Mev and obtained the same 
resonance curve as before. This showed that the 
y-rays from the 1.43-Mev resonance have an 
energy greater than 2.5 Mev and the previous 
experiment had already shown that the y-rays 
from the resonance were less than 4 Mev. All 
these results indicated that the radiation from 
the 1.43-Mev level is from C™ and that the 
y-rays have the same energy as those from the 
other broader levels. 

The intensity of the y-rays from a_ thick 
graphite target was compared with that obtained 
from the bombardment under similar conditions 
of a thick crystal of CaF, by 1.5-Mev protons. 
The y-rays from C+H? were found to be 3.0 
times as intense as those from CaF, when an 
electroscope shielded by one cm of iron was 
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Fic. 8. Integral range curve for recoil protons 
from the neutrons of C%+H?, (Fig. 7). 
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used to detect the y-rays. The absolute yield of 
the y-rays from carbon was determined at 1.0 
Mev by the use of a Geiger counter, assuming 
that it counted y-rays with an efficiency of two 
percent. For 1.0-Mev deuterons the yield of 
y-rays from carbon was 17X10° quanta per 
micro-coulomb of deuterons. This shows that 
the y-rays have about 4 times the intensity of 
those from CaF:;+H! at this same bombarding 
voltage. 


THE NEUTRONS FROM C+H? 


We have investigated the excitation curve for 
the neutrons from carbon bombarded by deu- 
terons. To detect the neutrons we have used a 
Wulf type electroscope filled with hydrogen at a 
pressure of 7 atmospheres. We have also studied 
the radioactive N"“ which is formed whenever a 
neutron is emitted according to reaction (2). 
The electroscope was placed in the forward 
direction to the incident deuteron beam and it 
gave, therefore, the excitation curve for neutrons 
in the forward direction. The amount of radio- 
active N™ formed in a thin paraffin target was 
measured by means of a thin-walled Geiger 
counter placed outside a thin window on the 
target tube. The procedure was to bombard the 
target for ten minutes, then shut off the high 
voltage and follow the activity of the target. 
The half-life of the N™ activity was found to be 
ten minutes in agreement with the known half- 
life of 9.93 minutes.” 

Figure 1 of reference 17 shows the excitation 


% A. G. Ward, Proc. Camb. Phil. Soc. 35, 523 (1939). 
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curves obtained from the ionization currents in 
hydrogen and from the radioactivity of the N*. 
Since some of the ionization in hydrogen is due 
to y-rays, a correction had to be made to com- 
pensate for this effect. By using a radium source 
filtered by 1.0 cm of lead, we found that the 
y-ray intensity measured by the electroscope 
filled with hydrogen was only 4.75 percent that 
of the argon-filled electroscope which was used 
in the y-ray experiments. The total intensity 
from a carbon target bombarded by 1.525-Mev 
deuterons was measured by using the electro- 
scope filled with argon and then with hydrogen. 
From the ratio of the two ionization effects it 
was possible to calculate that the y-rays were 
responsible for 19 percent of the ionization in 
hydrogen at 1.525 Mev. A corresponding correc- 
tion was made at all other voltages. Curve (3) 
gives this corrected ionization due to the neu- 
trons alone. Both curves (1) and (3) show 
resonances at 0.92, 1.16 and 1.30 Mev. These are 
the same resonances as those found for the 
emission of y-rays. However the y-ray resonance 
at 1.43 Mev does not show up as a resonance for 
neutron emission. Above 1.43 Mev the hydrogen 
ionization data show a resonance at 1.74 Mev 
and 1.82 Mev. The resonance at 1.74 Mev 
partially disappears after the y-ray correction is 
made, and it is uncertain whether the resonance 
at 1.74 Mev is truly a resonance for neutron 
emission. 

Although curves (1) and (3) of reference 17 
agree in showing the same resonances, the rela- 
tive heights of the resonances are different. The 
differences might be explained by supposing that 
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Fic. 10. Short range (1.1 cm) protons from C"+H?. The 
intense group of particles at 1.6 cm are deuterons from the 
incident beam which have been scattered at 90° into the 
cloud chamber. The closed circles have ordinates 6 times 
the value shown on the graph. 


NUMBER OF TRACKS 


neutrons from the different excited levels of the 
intermediate nucleus have different angular dis- 
tributions; whereas curve (1) gives the total 
yield of N™ (and neutrons), curve (3) gives the 
yield of neutrons only in the forward direction. 
Experiments to study these angular distributions 
would be of interest. 

Since it seemed likely that the 5.5-Mev y-rays 
from C™ were produced in reaction (3), we 
looked for low energy neutrons when a thick 
target of 23-percent C™ was bombarded by 
1.4-Mev deuterons. The neutron energies were 
measured by the method of finding the energy 
distribution of the recoil protons photographed 
in a cloud chamber containing methane. The 
neutrons were observed at an angle of 90+10° 
to the direction to the deuteron beam. Two 
thousand stereoscopic photographs were taken 
on which there were numerous recoil protons. 
One thousand and ten of these protons were in the 
forward direction (0—10°) and their track lengths 
were measured. The energy distribution of these 
recoil protons is given in the curve of Fig. 7. 
The group of neutrons at 0.9 Mev is the group 
from reaction (2) which was previously measured 
by Bonner and Brubaker. The Q value obtained 
in the present experiment was —0.19+0.05 
Mev. This agrees with the value of Q.= —0.25 
Mev which was obtained by Bonner and Bru- 
baker. The group of neutrons at 1.5 Mev was 
not observed by Bonner and Brubaker in their 
original experiments. For this reason and be- 
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cause of energy considerations this group must 
be from C*. An accurate calculation of the Q 
value for this group of neutrons was calculated 
by the method of Livingston and Bethe.*® The 
integral number vs. range curve of the recoil 
protons is given in Fig. 8. The extrapolated 
range of the recoil protons is 4.4+0.2 cm of 
standard air. The energy of the neutrons coming 
from the surface of the target at 90° to the 
incident beam was found to be 1.49+0.03 Mev. 
The disintegration Q value was then calculated 
from the relation Q3=(15/14)E,—(6/7)Ep. 
was found to be 0.40+0.05 Mev. The best value 
of the energy evolved in reaction (3), when the 
products are left unexcited, is that calculated 
from the masses of the atoms concerned and these 
give Q;=5.5+0.2 Mev. The neutron group at 
1.5 Mev must then leave the N™ nucleus in an 
excited state at 5.1+0.3 Mev. This excited N™ 
nucleus would give a y-ray of this energy when 
it returns directly to the ground state. This 
group of neutrons from C" then seems to explain 
the 5.5+0.2 Mev y-rays which are also produced 
in the disintegration of C™. 


PROTONS FROM C+H? 


Since the y-rays and neutrons from the bom- 
bardment of carbon by deuterons show many of 
the same resonances, we might expect the 15-cm 
protons from reaction (1) also to show the same 
resonances. These resonances correspond to ex- 
cited states in the intermediate N'* nucleus, 
and we might expect the emission of neutrons, 
protons and y-rays to be competing processes. 

The excitation curve for the emission of 15-cm 
protons was studied with an ionization chamber 
2.5 cm in depth connected to a linear amplifier. 
Such a deep chamber was used so that the 
protons would give large pulses and would be 
counted far from the end of their range. This 
was useful because the range of the protons from 
reaction (1) changes with bombarding energy. 
As an additional precaution to insure the count- 
ing of all protons, absorbers were added as the 
bombarding voltage was increased so that the 
protons were always counted at the same dis- 
tance from the end of their range. Figure 9 


% M.S. Livingston and H. A. Bethe, Rev. Mod. Phys. 9, 
281 (1937). 
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gives the excitation curve for the 15-cm protons 
from reaction (1). The curve shows the reso- 
nances at 0.92 Mev and 1.16 Mev that were 
found in the emission of y-rays and neutrons. 
The proton curve shows a resonance at 1.23 Mev 
which did not appear for y-rays or neutrons. 
The resonance at 1.30 Mev which was shown by 
both neutrons and y-rays appears to be missing 
in the proton excitation curve. Above 1.35 Mev 
the proton curve shows a rise of about 25 percent 
and then stays nearly horizontal up to 1.7 Mev, 
where it begins to decrease rapidly. 

The relative yield of 15-cm protons and neu- 
trons at 0.92-Mev bombarding energy was ob- 
tained by counting the number of N™ atoms 
formed and comparing it with the number of 
protons from the same thin target. At 90° to the 
bombarding direction, the ratio of protons to 
positrons from N"™ is 1.9. The number of y-rays 
in comparison to the number of protons was 
also obtained at 1.60 Mev. The number of 
counts produced by y-rays in a single Geiger 
counter was obtained at the same time that the 
15-cm protons were counted by the linear 
amplifier. After allowing for the differing solid 
angles and for the efficiency of the Geiger 
counter for counting 3-Mev y-rays (2 percent), 
we obtained 1.0 as the ratio of the number of 
y-rays to the number of protons. This shows 
that all three processes, the emission of neutrons, 
of 15-cm protons, and of 3.0-Mev y-rays are 
about equally probable. 

The apparent explanation of the 3-Mev y-rays 
from C” is that they come from a C® nucleus 
which is left excited to 3 Mev after the emission 
of a short range proton group. Consequently, we 
have looked for such a group of protons. The Q 
value for this group of protons should be equal 
to Q,—E, which is numerically equal to 2.71 
—3.0=—0.29+0.25 Mev. Since the expected 
proton group would have a range less than that 
of the scattered deuterons, we used an experi- 
mental arrangement which could detect particles 
with shorter ranges than the scattered deuterons. 
The deuteron beam was allowed to pass through 
an aluminum foil with a stopping power of 6.5 
mm of air into a gas target chamber which was 
filled with methane at a pressure of 2.5 cm. The 
protons and scattered deuterons were observed 
in a cloud-chamber at an angle of 90° to the 


direction of the deuteron beam. The disintegra- 
tion particles entered the cloud chamber through 
a Nu-Skin foil which had a measured stopping 
power of 1.3 mm of air. The cloud chamber was 
operated at a pressure, when expanded, of 0.5 
atmosphere of helium and water vapor. The high 
voltage apparatus was adjusted to give a beam 
of 1.7-Mev deuterons during this experiment, 
and the deuterons, after passing through the 
aluminum foil and 3.9 cm of methane, had an 
energy of 1.47 Mev at the center of the gaseous 
target. The part of the gaseous target from which 
it was possible to get disintegration particles 
entering the cloud chamber was 1 cm deep in 
the direction of the deuteron beam; thus the 
CH, target had a thickness equivalent to ap- 
proximately 0.03 cm of air. 

At this bombarding voltage we obtained about 
one track per expansion coming from the gaseous 
target. The scattered deuterons had an actual 
track length of 10.2 cm in the cloud chamber. 
We took a total of 4000 cloud-chamber pictures 
with this experimental arrangement and on these 
photographs we measured 2030 tracks. The 
length of the tracks in helium was converted to 
range in standard air by taking the stopping 
power of He equal to 0.20. To the range in the 
cloud chamber was added 0.39 cm, which was 
the air equivalent of the Nu-Skin foil and the 
thickness of methane absorber between the 
deuteron beam and the foil. The resulting range 
distribution of the tracks in terms of standard 
air is given in Fig. 10. The 15-cm protons from 
C® were observed in the cloud chamber, but 
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Fic. 11. Long range protons from a 23-percent C™ target 
bombarded by 1.00-Mev deuterons. The rise at 18-cm 
range is from the long range (15 cm) protons from C"+H?. 
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& Fic. 12. The excitation curve for 
> the long range (50 cm) protons from 
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DEUTERON'- ENERGY IN MV 


these went entirely across the chamber, so they 
are not indicated on the curve of Fig. 10. The 
group of particles with a range of 1.61 cm are 
the scattered deuterons. The expected range of 
the scattered deuterons was computed using the 
1938 Cornell range-energy relation and was 
found to be 1.50 cm which agrees to within 7 
percent with the experimental range of this 
group. The group of scattered deuterons appears 
to be symmetrical and has a half-width of 0.18 
cm of standard air. Besides the strong group of 
scattered deuterons there is a weaker group of 
particles with a range of 1.10+0.10 cm. The 
ratio of the number of 1.1-cm particles to the 
number of 15-cm protons was found to be 0.44. 
Since we previously found that the number of 
y-rays was 1.0 times the number of 15-cm 
protons, we should have expected that there 
would have been 1.0 times as many short range 
protons as 15-cm protons, if the short range 
proton groups had a spherical angular distribu- 
tion. The fact that the experimental ratio is 0.44 
indicated that the number of short range protons 
must be a factor of 2 to 3 greater in the forward 
direction. Another possible explanation of the 


'1.1-cm group is that it arises from a small 


amount of molecular hydrogen ions in the beam 
of mass 2. We found that our ratio of current in 
the mass 2 beam was usually about 2 to 3 times 
as large as the current in the mass 4 beam. We 
also had 7 percent as many protons (mass 1 
beam) as deuterons (mass 2). Consequently we 


would expect that about 1 percent or less of the 
mass 2 beam was composed of molecular hydro- 
gen atoms each of which had half the bombarding 
energy. The energy of these molecular hydrogen 
ions at the scattering target was calculated, and 
the range of these scattered protons was com- 
puted to be 0.70 cm of air. This eliminates the 
possibility that the group at 1.10 cm could be 
caused by scattered molecular hydrogen ions. 

The Q value corresponding to this 1.10-cm 
group of protons was calculated from the rela- 
tion: Qj = (14/13)E4—(11/13)Ep and a value of 
—0.52+0.07 Mev was obtained. From these 
data the energy of the y-ray should be 2.71+0.53 
= 3.24+0.08 Mev which agrees with the experi- 
mental value of 3.0+0.2 Mev. Since each low 
energy proton is accompanied by a quantum of 
y-radiation, the excitation curve for the produc- 
tion of low energy protons is the same as for the 
3-Mev y-rays (reference 16). 

It is interesting to note that the ratio of dis- 
integration protons from carbon to the number 
of scattered deuterons was 0.37, which indicates 
that at this angle a disintegration is nearly as 
probable as an elastic scattering. 

We have also investigated the protons that 
come from C" according to reaction (4). In order 
to be able to count weak groups of particles in 
the presence of neutrons we have used a deep 
ionization chamber (2.5 cm) connected to a 
linear amplifier. When the bias is set so that only 
the largest pulses are counted, the background 
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due to neutron recoils is found to be very small. 
We bombarded a thick target of 23-percent C™ 
with 1.00-Mev deuterons and obtained the curve 
shown in Fig. 11. The group of protons from C" 
has an extrapolated range of 53.7+2.6 cm. The 
value of Qs was found from this value of the 
extrapolated range to be 6.09+0.2 Mev. This 
agrees with the value of Q, calculated from the 
masses which is 6.13++0.19 Mev. 

No other group of particles from C™ was ob- 
tained with ranges greater than 20 cm although 
a group with an intensity of 75 that of the 50-cm 
group could have been detected. We also looked 
for protons from C' below the 15-cm group of 
protons from C™. No groups were observed, but 
only groups considerably more intense than the 
50-cm group could have been detected in this 
region. The experiment indicates that there are 
probably no excitation levels in C™ below 2.8 
Mev. The 3-cm group of alpha-particles from 
carbon was observed both with an ordinary 
carbon target and with a 23-percent C™ target. 
Since the number of alpha-particles was much 
greater from the C" target this definitely shows 
that this group had been properly attributed*’ to 
the reaction : 


Qs. (5) 


The ratio of the intensity of the 50-cm group 
of protons to the 15-cm protons from C” was 
determined with a target in which the C™ had 
been concentrated. With this target and 1.22- 
Mev deuterons, the ratio of the intensity of the 
55-cm group to the intensity of the 15-cm group 
of protons was 0.11 percent. A target in which 
the C was depleted by approximately the same 
factor was then bombarded and the ratio was 
only 0.029 percent. This showed that the C™ 
had been concentrated in the first target relative 
to the second by a factor of 3.8, whereas the 
separation factor had been calculated from the 
dimensions of the thermal diffusion apparatus as 
3.7. The average of these two ratios (0.07 per- 
cent) is consequently the ratio of the yields from 
ordinary carbon at 1.22-Mev bombarding energy. 
This ratio would vary considerably with different 
bombarding voltages since the yields for both 
groups of protons show resonances. 


*7J. D. Cockcroft and W. B. Lewis, Proc. Roy. Soc. 
A154, 261 (1936). 


The excitation curve has also been obtained 
for the 50-cm protons. A target of 2-percent C™ 
prepared from methane and with a thickness of 
30 kev, was used in this experiment. The deep 
ionization chamber was again used to count 
protons, and the bias on the recording circuit 
was made low enough to count all protons 
passing through the ionization chamber. Ab- 
sorbers were also added as the bombarding 
voltage was increased to insure that the protons 
were always counted at the same distance from 
the end of their range. The excitation curve is 
given in Fig. 12. The indicated error is the 
statistical error increased by a factor of two, 
since it was necessary to make observations in 
two separate runs, which were joined at 1.45 
Mev. The curve shows a maximum at 1.52 Mev. 
This may be interpreted as a resonance in the 
yield which would indicate an excited state in the 
intermediate nucleus (N™*) at 17.5 Mev. 


DISCUSSION OF THE RESULTS 
cr 


The sharp resonances obtained from the dis- 
integration of C” by deuterons are surprising 
because of the fact that the intermediate nucleus 
(N") is excited to 11 Mev and might be expected 
to emit a neutron or a proton in such a short 
time that the level in N'* would be wide. The 
width of the level is given by the uncertainty 
relation AEAT ~h. The level at 1.30 Mev had a 
width of 40 kevand consequently the intermediate 
N** nucleus must exist for a time 


1.04 10-27 
1,60 100.040 


= 1.6X10~*° sec. 


This shows that it takes more than 1.610-*° 
sec. to emit a neutron with an energy of 0.8 Mev 
from the excited N“* nucleus. The lifetime of the 
intermediate N™* nucleus can be estimated from 
the relation 


AT=R/2, 


where v is the velocity of the neutron in the N“* 
and R is the nuclear radius. If we take R=4 X10" 
cm and a neutron velocity of 1/20 that of light, 
we get a value of AE=1.2 Mev. This is more 
than 30 times the observed width and so it may 
be necessary to invoke a selection rule to explain 
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this result. The emission of neutrons and 15-cm 
protons from the 1.43-Mev resonance level in 
N'** must have a strong selection rule forbidding 
such modes of disintegration. 

The experiments with C” show what the 
three modes of disintegration, the emission of 
(1) a neutron, (2) a 15-cm proton, and (3) a 
short range proton, are all about equally probable 
from the levels in N'* corresponding to the 
resonances at 0.92 and 1.16 Mev. However the 
15-cm protons show a resonance at 1.23 Mev 
which is not shown by either of the two com- 
peting processes. It seems possible, however, 
that this level may be the same as that shown for 
neutron and y-ray emission at 1.30 Mev. A shift 
in the position of the maximum of 0.07 Mev 
might be produced through a different inter- 
ference with the background, or it might be 
caused by the theoretical considerations dis- 
cussed by 

The 15-cm proton group seems to show the 
resonance at 1.74 Mev as shown by the sharp 
decrease in the curve above 1.74 Mev, and 
another broad resonance at from 1.5 to 1.6 Mev. 
A broad resonance in this region is probably 
needed to explain the plateau which is experi- 
mentally observed from 1.45 to 1.75 Mev. 


C1 


In the disintegration of C by deuterons the 
competing processes are reaction (3), (4) and 
reaction (5). Unlike the case of C” the three 
competing processes for C™ have very unequal 
probability. The 5.5-Mev y-rays (consequently 
1.5-Mev neutrons) are approximately 70 times 
as numerous as are the 50-cm protons. The 
alpha-particles from reaction (5) were found to 
be roughly 50 times as probable as the 50-cm 
protons. Thus it seems that the emission of low 


28 G. Breit, Phys. Rev. 58, 1068 (1940). 


energy neutrons and alpha-particles from C™ are 
about equally probable but that the emission of 
50-cm protons is very unlikely. 

The excitation curves for 5.5-Mev y-rays and 
50-cm protons both show the same resonance at 
1.55-Mev bombarding voltage. This indicates 
that an excited state in the intermediate (N™*) 
nucleus at 17.5 Mev may break up either into 
low energy neutrons or 50-cm protons. It is 
surprising that an isolated resonance at this high 
excitation of the intermediate N' nucleus is 
observed. 

The principal difference in the excitation 
curves for 50-cm protons and 5-Mev y-rays is 
that the emission of y-rays appears to come 
mainly from a continuum, while the emission of 
50-cm protons comes mainly from the resonance 
at 1.55 Mev. The fact that very few of the 50-cm 
protons come from the continuum may explain 
why this reaction is so much more improbable 
than the emission of 5-Mev y-rays. 

The width of the level in the N'* at 17.5 Mev 
may be calculated as was done for the level in 
N"™* by the relation AEAT=h and AT=R/2. 
The width of this level is fixed by the neutron 
emission which is the most probable reaction. 
A calculation gives the value of AE=1.7 Mev for 
the width of this level. It may not be necessary 
to invoke the use of selection rules forbidding 
disintegration into neutrons and protons from 
this level to explain the experimental half-width 
of this level as 0.25 Mev. It is reasonable that 
the actual width of the level is smaller than the 
estimated width because the estimate does not 
take into account the internal rearrangements 
necessary for concentrating the necessary energy 
on either a neutron or a proton. 

It would be interesting to determine the 
excitation curve for the alpha-particles to see if 
they show the same resonance at 1.55-Mev bom- 
barding energy since they come from the same 
N* nucleus. 
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; and to have a maximum energy of 6.7+0.3 Mev. This indicates that the highly excited states 
‘Ss of Ne®® at 9.0 and 10.1 Mev, which were found in the reaction F!*(d, n) Ne®*, break up into O" 


are The Energy of the y-Rays from the Disintegration of Fluorine by 
1 of Protons and Deuterons 
d W. E. Bennett, T. W. BONNER AND Bos E. Watt , 
an The Rice Institute, Houston, Texas ' 
> at (Received March 31, 1941) ' 
ites 
aE. The maximum energy of the y-rays from the disintegration of a thick target of fluorine by 
nto protons has been determined at 0.90 Mev and 1.7 Mev. The energy of the y-rays was found : 
tis * to be the same at both bombarding voltages. The y-rays from the disintegration of fluorine by i 
‘igh deuterons have been studied with 1.3-Mev deuterons. The y-rays were found to be complex | 
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and He‘ instead of emitting y-rays. 


HE y-rays from the disintegration of 
fluorine by deuterons were discovered by 
McMillan! in 1934. From absorption of these 
y-rays in different elements, he obtained a 
value of their quantum energy of 5.4 Mev. 
Hafstad and Tuve? first found resonance maxima 
for these y-rays at a proton bombarding energy 
of 0.328, 0.892, and 0.942 Mev. Later Bernet, 
Herb and Parkinson* found many additional 
resonances from 0.33 up to 2.2 Mev. Above 
about 1.1 Mev the resonances were not com- 
pletely resolved, which would indicate that the 
resonance levels were becoming either wider or 
closer together. The energy of the y-rays has 
been carefully studied by Lauritsen, Lauritsen 
and Fowler‘ using the cloud-chamber technique. 
They found that the quantum energy of the 
y-rays was the same from all the resonances 
below 1.4 Mev and was equal to 6.2+0.1 Mev. 
The y-rays from the disintegration of fluorine 
by protons have been shown by McLean, Becker, 
Fowler and Lauritsen> to come from the reaction: 


F9+H!—(Ne2*) + He!+@! (1) 


They discovered the short range alpha-particles 
that are produced in reaction (1) and showed 
that these alpha-particles exhibited the 330-kev 
resonance. 


1E. McMillan, Phys. Rev. 46, 868 (1934). 


*L. R. Hafstad, N. P. Heydenburg and M. A. Tuve, 
Phys. Rev. 49, 866 (1936). 

+E. J. Bernet, R. G. Herb and D. B. Parkinson, Phys. 
Rev. 54, 398 (1938). 

‘ T. Lauritsen, C. C. Lauritsen and W. A. Fowler, Phys. 
Rev. 59, 241 (1941). 

*W. B. McLean, R. A. Becker, W. A. Fowler and C. C. 
Lauritsen, Phys. Rev. 55, 797 (1939). 


Because the character of the resonances 
above 1.1 Mev appears to change, we thought 
that some of these higher resonances might be 
caused by a different reaction and give y-rays 
of different energy. Consequently we have 
investigated the energy of the y-rays at a 
bombarding energy of 1.7 Mev. 

When fluorine is bombarded by deuterons, 
three different reactions are known to take 
place. Lewis, Livingston and Lawrence® in 1933 
reported the emission of a group of alpha- 
particles of 3.8-cm range. Later Burcham and 
Smith’ carefully studied the reaction: 


H?-(Ne?!*) + He! (2) 


and found five groups of alpha-particles with 
Q values equal to 9.84, 9.01, 6.89, 6.07 and 5.35 
Mev. The lower energy alpha-particles are 
emitted when the O" is left in excited states at 
0.83, 2.95, 3.77 and 4.49 Mev. 

Bower and Burcham*® found 5 groups of 
protons with ranges of 13.0, 18.8, 21.0, 24.2 and 
31.0 cm from the bombardment of fluorine by 
deuterons. These protons come from the reaction: 


(3) 


The lower energy proton groups correspond to 
excitation energies of 0.7, 1.0, 1.35 and 1.9 Mev 
in the nucleus. 

Bonner® has investigated the neutrons from 


*G. N. Lewis, M. S. Livingston and E. O. Lawrence, 
Phys. Rev. 44, 55 (1933). 

7™W.E. Burcham and C. L. Smith, Proc. Roy. Soc. A168, 
176 (1938). 

8 J. C. Bower and W. E. Burcham, Proc. Roy. Soc. A173, 
379 (1939). 

* T. W. Bonner, Proc. Roy. Soc. A174, 339 (1940). 
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MM OF ALUMINUM ABSORBER 


Fic. 1. The absorption curves for the secondary elec- 
trons produced by y-rays from F+H!. The closed circles 
refer to the measurements taken at a bombarding voltage 
of 0.9 Mev, while the open circles refer to the data taken 
at 1.7 Mev. The ordinates of the upper curve have been 
multiplied by 5 in order to magnify the curve near the 
end point. 


the reaction: 
(4) 


and has found seven groups of neutrons. The 
neutrons with energy less than the maximum 
indicate excited states in Ne*® at 1.5, 4.2, 5.4, 
7.3, 9.0 and 10.1 Mev. The excited states in Ne”® 
at 5.4, 7.3, 9.0 and 10.1 Mev are energetically 
unstable against alpha-particle emission and so 
they may either break up into O" and an alpha- 
particle or radiate a photon: We have investi- 
gated the energy of the y-rays from fluorine 
bombarded by deuterons to try to see how these 
highly excited states in Ne®® lose their excess 


energy. 
EXPERIMENTAL ARRANGEMENT 


The arrangement for bombarding targets with 
high energy protons and deuterons was the same 


as that described in the preceding paper.'® The 
energy of the y-rays was determined by observ- 
ing the maximum range of the Compton elec- 
trons in aluminum placed between coincidence 
Geiger counters. This method was first used by 
Bothe" and later was improved by Dee, Curran, 
and PetrZilka.” 

In the present work we used two alcohol-filled 
Geiger counters and recorded the passage of the 
electrons through both counters. They were 1.6 
cm in diameter and 6 cm long and were placed 
far enough apart that there was space for 2 cm 
of absorber between the counters. Each wall of 
the counters had a thickness equivalent to 1.25 
mm of aluminum. The coincidence circuit had a 
resolving time of 2.0X10-* sec. and would 
count coincidences properly when there were 
20,000 single counts in each counter per minute. 


RESULTS 


The y-rays from ThC" were first studied in 
order to check the apparatus and in order to get 
a calibration value on the range of the electrons 
from 2.6-Mev y-rays. The maximum range of the 
electrons was found to be 4.3 mm of aluminum. 
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ENERGY OF GAMMA-RADIATION IN MV 


Fic. 2. This curve gives the range in aluminum of second- 
ary electrons as a function of the energy of the y-rays which 
roduced them. It differs slightly from the curve determined 
ot Curran, Dee and Petrzilka taken under different ex- 
perimental conditions (see reference 12). 


10W. E. Bennett, T. W. Bonner, E. Hudspeth, H. T. 
Richards and B. E. Watt, Phys. Rev. 59, 781 (1941), 
(this issue). 

1 W. Bothe and J. Becker, Zeits. f. Physik 76, 421 
(1932); F. Rasetti, Zeits. f. Physik 97, 64 (1935); G. 
Bernardini and S. Franchetti, Ricerca Scient. 8, 406 (1937); 
S. Kikuchi, K. Husimi and H. Aoki, Proc. Phys. Math. Soc. 
Japan 18, 35 (1936). 

2S. C. Curran, P. I. Dee and V. PetrZilka, Proc. Roy. 
Soc. A169, 269 (1939). 
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Fic. 3. The absorption curve for the secondary electrons 
produced by 7-rays from F+H?. 


Next the y-rays from a thick target of CaF» 
which was bombarded by 0.9-Mev protons were 
investigated. The solid circles of Fig. 1 show the 
absorption curve for the secondary electrons 
from these y-rays. The maximum range of the 
electrons is 11.7+0.5 mm of aluminum, Above 
the maximum range of 11.7 mm the counting 
rate is greater than should be expected from 
chance coincidences between singles in the two 
counters. We made a test to see whether the 
counts with thick absorbers were chance coin- 
cidences by decreasing the current on the target 
by a factor of 2. If the counts beyond 11.7 mm 
of aluminum were chance coincidences they 
should have decreased by a factor of 4. Actually 
the counting rate decreased about a factor of 2, 
indicating that most of this background was 
caused by real coincidences. It seems probable 
that these coincidences are produced by a 
Compton electron in the second counter and its 
recoil photon which is scattered backward into 
the first counter. The intensity of the background 
on the absorption curve seems of the right order 
of magnitude to be explained by this process. 


Also this background usually drops by a factor 
of two when 14 mm of aluminum absorbers are 
added. This is about what is expected for the 
absorption of recoil photons since they would 
have an energy of 0.25 Mev. 

The range of the electrons produced by these 
fluorine y-rays gives us a calibration point on 
the range-energy curve at 6.2 Mev. Since Cur- 
ran, Dee and Petrzilka have shown that the 
range-energy relation is linear above 2.6 Mev, 
the two calibration points are sufficient to 
define a range-energy curve for our experimental 
arrangement. This curve is given in Fig. 2. 

The energy of the y-rays from a thick CaF, 
target which was bombarded with 1.7-Mev 
protons was next investigated. The absorption 
curve of the Compton electrons in aluminum is 
given by the open circles of Fig. 1. It is immedi- 
ately apparent from Fig. 1 that there is very 
little if any difference in the energy of the y-rays 
at these two different bombarding voltages. In 
view of the fact that a 5.5-Mev y-ray, which had 
an intensity of only 5 percent that of a 3-Mev 
y-ray which accompanied it, was easily detected 
by the same method (see preceding paper), we 
believe this experiment shows there is less than 
1 percent of radiation of 10 Mev or higher energy 
from a thick target of CaF, bombarded by 1.7- 
Mev protons. 

The y-rays from CaF; bombarded by deu- 
terons were studied at a bombarding voltage of 
1.3 Mev. The absorption curve for the secondary 
electrons from these y-rays is given in Fig. 3. 
The maximum range of the electrons is 12.7 mm 
of aluminum which gives a quantum energy of 
6.7+0.3 Mev. The thickness of aluminum where 
the counting rate falls to half-value is 1.9 mm 
for these y-rays, while the half-value for the 
homogeneous y-rays from F+H! was 1.9 mm 
of aluminum. This shows that the y-radiation 
from F+H? although no doubt complex, con- 
sists largely of the penetrating components. 

It seems natural to associate the 6.7+0.3-Mev 
y-rays found in this experiment with the excited 
level in Ne®® at 7.34+0.3 Mev. which was found 
from a study of reaction (4). The fact that no 
y-rays of energy higher than 6.7 Mev were ob- 
served indicates that the Ne? nucleus when 
excited to the 9.0- and 10.1-Mev levels breaks up 
into O'+ He! before it can emit a y-ray. 
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A Photographic Plate Spectrum of the Neutrons from the Disintegration of 
Lithium by Deuterons 
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The energy spectrum of the neutrons from Li+H? has been studied by the photographic 
emulsion technique. Confirmation is given of the broad 10.8-Mev neutron group reported by 
Stephens from He recoil data. Reasons are given for believing that not all the neutrons from 
this group correspond to the 2.8-Mev excitation level in the residual Be® nucleus. Evidence is 
also given for at least two higher (7.5 Mev and 10 Mev) levels in Be*. A discussion of the 
photographic emulsion technique is included, and a correction is calculated for the differing 
probability that long tracks be acceptable for measurement. The stopping power of the emul- 
sion for three different proton energies was found. 


INTRODUCTION 


HE energy distribution of neutrons from 

the disintegration of lithium by deuterons 
was first measured by Bonner and Brubaker.! 
A rather inhomogeneous distribution was found 
except for a neutron group at about 14 Mev 
which was attributed to the formation of Be® 
in its ground state: 


Mev. (1) 


Bonner and Brubaker’s distribution above about 
8 Mev was somewhat uncertain because the 
long range of the recoil protons (~200 cm) re- 
quired the use of mica sheets to stop them in 
their high pressure cloud chamber. The cor- 
rection (to compensate for the unequal prob- 
ability of observing different length tracks) was 
large and uncertain when the mica was in the 


chamber. 


By using He instead of H recoils, Stephens* 
was able to obviate the use of mica in the cloud 
chamber since the range of the recoil alpha- 
particles is much less than that of the recoil 
protons. Besides the homogeneous group from 
(1), he observed another group at about 10.8 
Mev which he ascribed to a wide three-Mev 
excitation level in Be*: 


Li’+H?—*Be®+n+12 Mev (2) 
*Be®—2 He'+3 Mev. 


1T. W. Bonner and W. M. Brubaker, Phys. Rev. 48, 


742 (1935). 
2 W. E. Stephens, Phys. Rev. 53, 223 (1938). 


In later work Staub and Stephens’ reported a 
plateau in the low energy neutron distribution 
which they attributed to the two-stage dis- 
integration, 


Li?’ Mev (3) 
He®’—He'+2+0.8 Mev. 


The data of Staub and Stephens’ above 8 Mev 
were uncertain since the neutron-helium scatter- 
ing cross section was not known for this region. 
At a lower energy (1.1 Mev) an anomaly in this 
cross section has been discovered.‘ This em- 
phasized the possibility that Stephens’ neutron 
group at 10.8 Mev might be spurious since 
higher levels in He® might give resonance neutron 
scattering which could be mistaken for true 
neutron groups. 

The object of the present experiment was to 
make a careful study of this neutron distribu- 
tion using, however, recoil protons instead of 
recoil alpha-particles. In this way information 
could be obtained concerning the reality of the 
10.8 Mev group found by Stephens. 

Recent work by Powell and others® has shown 
the value of the photographic emulsion tech- 
nique for studying nuclear processes. New 


a 935} Staub and W. E. Stephens, Phys. Rev. 55, 845 

19 H. Staub and W. E. Stephens, Phys. Rev. 55, 131 
5C. F. Powell and G. E. G. Fertel, Nature 144, 115 

(1939); C. F. Powell, Nature 145, 155 (1940); C. F. Powell, 

A. N. May, J. Chadwick, and T. G. Pickavance, Nature 

145, 893 (1940); T. R. Wilkins, J. App. Phys. 11, 35 (1940). 

For a general discussion, history, and bibliography on the 
hotographic emulsion technique see M. M. Shapiro, 
ev. Mod. Phys. 13, 58 (1941). 


796 


t 
\ 


ay 
| 
| 
| 
| 


NEUTRONS FROM DISINTEGRATION OF LITHIUM 797 


emulsions are now available which record indi- 
vidual proton tracks. Neutron energies can 
therefore be studied by measuring the ranges of 
the recoil protons which are produced when a 
neutron collides with a hydrogen nucleus of the 
emulsion. 

For neutron energy studies the emulsion 
technique has several advantages over the cloud- 
chamber method. Scattering material can be 
kept to a much lower minimum, and as the 
photographic plate is continuously sensitive, 
the recording of the neutrons requires only a few 
hours use of the disintegration equipment. 
Comparable experiments with cloud chambers 
would tie up the equipment for weeks or months. 
This advantage in recording the data is partially 
offset by the longer and more tedious measuring 
of the tracks with a microscope. For high energy 
neutrons the emulsion method has the advantage 
that it is equivalent to a cloud chamber of 
infinite diameter. Hence it is possible to avoid 
the uncertainty in the cloud-chamber data 
which is introduced by joining runs taken with 
different stopping powers and with different 
recoil nuclei. 


NOTES ON THE PHOTOGRAPHIC 
EMULSION TECHNIQUE 


(1) Emulsion used 


Although certain commercial process plates 
are suitable for recording individual alpha- 
particles, the only emulsion which was found to 
give really satisfactory recoil proton tracks was 
the special halftone plate which Ilford, Ltd. 
(London) coats to order with emulsions up to 
100 microns thick. More recently Ilford has been 
supplying 300 microns thick emulsions without 
the glass plates. Eastman Kodak Company also 
reports that they have developed a plate suitable 
for recording protons. 


(2) Exposure arrangement 

The x-ray intensity in the observation room 
of the Rice pressure Van de Graaff generator was 
low enough that the plates needed to be wrapped 
only in paper to exclude light. However, a thin 
aluminum sheet was used between the paper and 
the emulsion to prevent recording recoil protons 
which originated from the paper wrapping. This 


was done to save time in analyzing the plates, 
for, otherwise, many more tracks have to be 
examined before one is found which both begins 
and ends within the emulsion. For exposure the 
plates were placed 10 cm from the target with 
the target in the plane of the emulsion so that 
the neutrons entered the plate tangentially. A 
thick target of metallic lithium was bombarded 
by magnetically analyzed 1200-kev deuterons. 
A bombardment equivalent to 0.6 microampere 
hours gave a satisfactory number of tracks. 
Later a LiCl target was used for a second set of 
plates. This was done to determine the effect 
of a possible nitrogen contamination of the 
metallic lithium (see RESULTs). 


(3) Processing technique 


More than ordinary care should be taken in the 
processing of the emulsions, for it is important 
that background fog be kept to a minimum. The 
100-micron thickness of the special Ilford 
emulsions creates special processing problems. 
A strongly caustic developer (Eastman’s D-9, 
hydroquinone-NaOH) was used to soften the 
gelatine quickly and thus allow the developer 
to reach the bottom parts of the thick emulsion. 
A thorough rinse and immersion in a hardening 
bath helps condition the gelatine for the pro- 
longed fixation and washing which is required. 
A 100-micron thick emulsion takes about an 
hour to clear in a fresh fixing bath. 

If the thick emulsions are exposed to dry air 
after the processing, they tend to crack or to 
curl up from the edge of the plate and even pull 
glass lamina off with them. Therefore the de- 
veloped plates must be kept in a humidor. 
Application of a film of collodion aids somewhat in 
preventing this dehydration and cracking. The 
effect of dehydration on these thick emulsions 
makes them unsuited for long use in a vacuum. 


(4) Viewing and measuring technique 


A compound microscope with a well-corrected 
objective of large numerical aperture is needed 
for satisfactory examination of the tracks. 
For short tracks a calibrated eyepiece could be 
used for measurements. However, a more con- 
venient and accurate measuring arrangement 
was required since the track lengths of the recoil 
protons extended up to 10 times the field of view. 
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CASE A (S?d/SIN®,) 


CASE B, (S*d/siINe, 


| 


CASE B; (S*CASE B,) 


Fic. 1. Diagrams for the calculation of a correction for the differing probability that long tracks end within the 
emulsion. See text for explanation. 


A satisfactory arrangement was made by mount- 
ing the focusing equipment and drawtube of a 
compound microscope upon the bed of a measur- 
ing microscope. 

A fluorite oil immersion objective (N.A. 1.3; 
100 X) was used to give a small depth of focus® 
so that a decision could be made as to whether a 
given track was wholly contained in the emul- 
sion. To give the needed flatness of field two 
paired 7.5X hyperplane eyepieces were used in 
a binocular tube attachment. The binocular 
arrangement is almost a necessity where many 
measurements are to be made; although the 
stereoscopic effect is negligible, the saving of 
eyestrain is quite important. 

A mechanical stage reading to 0.1 mm was 
used to manipulate the plate in the search for 
measurable tracks. The actual length measure- 
ment was accomplished by motion of the measur- 
ing microscope bed parallel to the incident 
neutron direction. As a precaution against 
remeasuring the same track (important when 
several are in the field of view) settings of the 

* Measured to be about 4 microns by finding the length 


in focus at one time for tracks which went completely 
across the emulsion and hence whose angle could be found. 


mechanical stage were used as track coordinates 
and recorded along with the track length. 


CRITERIA FOR MEASURABLE TRACKS 


The most accurate method of measurement of 
neutron energies by proton recoils is to measure 
only those tracks which make an angle with the 
incident neutron beam so small that the recoil 
proton may be considered as having received 
the total energy (E,) of the colliding neutron. 
Let 6 be the angle between the recoil proton and 
the incident neutron; then the energy (E,) of 
the recoil proton is 


E,=E, 08°. 


To insure that the recoil proton gets essentially 
the neutron’s energy, @ was limited to about 12°. 
Instead of measuring a recoil track only if it 
were in a cone of half-angle 6=12° about the 
neutron direction, it was found easier to apply 
the angle criterion for a square pyramid whose 
half-angle at the vertex was 10°. Thus only 
tracks scattered into the solid angle subtended 
by ABCD of Fig. 1 satisfy the angle criterion. 
Crosshairs in the eyepiece were used to test the 
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azimuthal angle in the emulsion plane. Since the 
microscope’s depth of focus was known, the 
latitude angle to the emulsion plane was easily 
tested by noting whether a certain minimum 
track length was in focus at one time. 

The other criterion of acceptability of a track 
for measurement was that the total track length 
must be wholly contained within the emulsion. 
This was determined by focusing above and 
below the endpoints of the track and noting 
whether there were still background grains in 
focus. This is where an objective of large nu- 
merical aperature is needed to provide a very 
small depth of focus. 


CORRECTION FOR THE DIFFERING PROBABILITY 
THat LonGc Tracks BE ACCEPTABLE 
FOR MEASUREMENT 


Because of the finite thickness of the emulsion, 
it is evident (Fig. 1) that, for a given latitude 
angle, many more of the short than the long 
tracks will end within the emulsion. Hence in 
order to get the correct relative number of short 
and long tracks it becomes necessary to apply a 
correction for this differing probability that a 
long track be acceptable for measurement. 

Consider the protons scattered from some 
point P at a distance x from the bottom of the 
emulsion. Let des =number of protons scattered 
per solid angle dw, 6.=maximum azimuthal and 
maximum latitude angle which a recoil proton 
can make and still be acceptable for measure- 
ment, S=track length of the recoil proton, and 
d=emulsion thickness. Also let o and o’ be the 
number of protons of given S which are scattered 
from P into the solid angle subtended by ABCD 
and afSyé, respectively (Fig. 1). Consider a 
single recoil proton of ¢. Then the probability 
(W) that it will also end within the emulsion 
(and hence be acceptable for measurement) is 


W=o'/c. 


¢ is calculated upon the assumption that the 
scattering is spherically symmetric in the center 
of gravity coordinates. In laboratory coordinates 
this becomes 


do = 2k sin@ cosédéd¢ 


and as 6 is small 
6” 
c= sf f 2k sin®@ = 8k sin*4o, 
0 0 


since sin@’ sec¢. 
For the calculation of o’ two different ranges 
of track lengths must be considered: 


Case A: 

S=d/sin@>. For this length S there are no 
tracks with latitude angle greater than @) which 
are wholly recorded in the emulsion. o’ then 
turns out to be independent of the position of 
the scattering center P. From Fig. 1 it is seen 


that 
or” 0 $1 0 0 0 
Xsiné 
=4k(d/S) sino, 
since 
sind; = (x/S) seco 
sin6;; = secd 
=[(x—d)/S] 
tand; = (S/x) sin®o 
tangy =[.S/(x—d) ] sino. 
Hence 
Wa=0'/c=d/(2S sin6o) 
when 
S =d/ 
Case B: 


Here is dependent upon the 
position of the scattering center P. However, 
upon averaging over all positions the final 
result is 

Wp=1—(S sin6o)/2d 
when 
S=d/sin@o. 


In obtaining the above expression two sub- 
cases (i) and (ii) are considered: 

(i) SSd/sin®)> but of sufficient length that 
o’<o for all possible positions of P. Let g 
and h (Fig. 1, case B;) be at distance x;=d—S 
sin@) from the surfaces of the emulsion. Hence 
when P is interior to gh, i.e., xxSxS(d—x,), 
the acceptance probability W,, is the same as 
in case A, i.e., W,,= Wa. The acceptance proba- 
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Fic. 2. Variation of the correction 


3 4 5 6 7 8 9 10 
ENERGY IN MEV 


bility W,; (for x<x,;) then remains to be 
calculated and averaged with W,,. From Fig. 1, 
case B;, it is seen that for x<x;, o’ =(¢/2)+0” 
where o”’ is the number of protons scattered from 
P into the solid angle subtended by yée¢. Hence 


w/2 
= 2( f f +f f ) sin@ cosédédq@ 
or 0 0 0 


=4k(x/S) sin 0. 
Thus the average value (W),;)« for O=x=x; is 
= a’/o = 3+2x;/(4S sino). 
Hence the final average for all possible positions 
of P: 0=x3(d/2) is 
(d/2—xi) Wa] 
 =1—(S sin@o)/2d. 
(ii) All tracks shorter than case B;. From Fig. 1, 
case B;;, x;;=S sind) and hence for x;;=x=(d/2) 
then o’=o and therefore W,,,=1. For x<x;; 
then o’ <o and hence the correction is that cal- 
culated in (i) for (W);)4 except that the expres- 
sion is averaged between 0 and x;; instead of 0 
and x;. This gives 


(Wap w (3) (4S sino) 


and thus the final average over 0Sx3(d/2) 
yields 
Wap w+ Winn 

=1—(S sinOo)/2d, 


which is the same average value as found in (i). 


factor (1/W) with energy of the 
recoil proton. W is prob- 
ability that a recoil proton of given 
energy and scattered from P into 
solid angle subtended by ABCD 
(Fig. 1) will end within the emul- 
sion. 


For convenience in applying this correction to 
the observed distribution, 1/W was _ plotted 
against the recoil proton energy instead of range 
(S). The effective emulsion thickness (d) was 
measured to be 80 microns. Figure 2 shows the 
correction curve for this value of d and for 
=0.2. 


RESULTS 


About 3500 acceptable tracks have been 
measured and are recorded in 0.4-Mev energy 
intervals (curve J, Fig. 3). 900 of those tracks 
(plotted separately in Fig. 4) were from the 
plates exposed to the LiCl target instead of the 
metallic lithium one. This was done to decide 
whether the neutron group at 5 Mev (curve JJ 
Fig. 3) could be from the disintegration’ 


Mev, 


since a black layer of Li;N appears to form on 
the surface of a freshly prepared metallic lithium 
target. That this cannot account for this neu- 
tron group is shown by the appearance of the 
same group on the data taken with the LiCl 
target. (The number of neutrons from Cl at 
this bombarding voltage is negligible.) 


7™W. E. Stephens, K. Djanab, and T. W. Bonner, Phys. 
Rev. 52, 1079 (1937). 

8 W. G. Shepherd, R. O. Haxby and E. L. Hill, Phys. 
Rev. 52, 675 (1937). 


8 

| 
| 
| 


ection 
of the 
-_prob- 
given 
P into 
4 BCD 
emul- 


on to 
otted 
ange 

was 
s the 
1 for 


n on 
ium 
neu- 
the 
LiCl 
1 at 


Phys. 
Phys. 


NEUTRONS FROM DISINTEGRATION OF LITHIUM 801 


T 


Fic. 3. Total data from plates 
exposed at 90° to lithium targets 


which were bombarded by 1200-kv 
deuterons. Curve J gives the meas- 
ured distribution of recoil protons. 
Curve JJ gives the inferred neutron 
distribution obtained after the cor- 


rection of Fig. 2 and the correction 
for scattering cross section had been 
made to the recoil proton data. 


CORRECTED NUMBER OF TRACKS (e) 


3 
ACTUAL NUMBER OF TRACKS (.) 


T\ 


. 


CALIBRATION OF PLATES 


The stopping power of the emulsion for two 
different proton energies was found by the use of 
the monochromatic neutrons from the d-d 
reaction* (Fig. 5). A 20-kv thick heavy paraffin 
target was bombarded by 700-kv deuterons. 
Plates were placed 12 cm from the target at 0° 
and 90° to the incident beam. From about 200 
tracks measured on each plate, a mean range of 
135 microns and 72 microns was found for the 
plates placed at 0° and 90°, respectively. Since 
recoils up to 12° were measured, this mean 
recoil proton range corresponds to an energy 
3E,00° less® than the calculated neutron energy 
from Bonner’s precise determination of Q=3.31 
Mev.'® Thus the corresponding recoil proton 
energies on the plates at 0° and 90° are calcu- 
lated to be 3.73 Mev and 2.60 Mev, respectively. 
This gives a stopping power of 6.6 microns per 
cm of air for the 3.73-Mev protons and 6.5 
microns per cm of air for the 2.60-Mev protons. 
This indicated that the stopping power relative 
to air might decrease slightly with increasing 


*Dr. E. Hudspeth aided in measuring the plates ex- 
posed to the d-d neutrons. A preliminary report has 
already been published: H. T. Richards, and E. Hudspeth, 
Phys. Rev. 58, 382 (1940). 

*M. S. Livingston and H. A. Bethe, Rev. Mod. Phys. 
9, 290 (1937). 

10 T, W. Bonner, Phys. Rev. 59, 237 (1941). 
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proton energies. Hence a calibration point for a 
high proton energy was desired. This was ob- 
tained from the mean range (1415 microns) of 
the longest range recoil protons produced by the 
neutrons from the disintegration of lithium by 
deuterons. The reaction producing these neutrons 
is known to have a Q=15.0 Mev.?:* Hence the 
mean range of this recoil proton group should 
correspond to an energy of 14.1—3E,6.?=13.8 
Mev since recoils with an angle up to 6).=12° 
were measured. This gives a stopping power of 
6.9 microns per cm of air for the 13.8-Mev protons. 
A smooth curve through these three values of 
the stopping power was taken to represent the 
variation of stopping power with energy. 

The range number curves were transformed 
into energy number curves by the use of the 
1938 Cornell range-energy relations. Curve /, 
Figs. 3 and 4, gives the measured recoil proton 
distribution. Curve JI of the same figures 
represents the inferred neutron distribution 
obtained by correcting the recoil proton distribu- 
tion for geometry (different probability of long 
tracks being acceptable) and for the variation 
of the neutron-proton scattering cross section 
with energy." 

"Wigner’s formula was used: 


+ where «,=0.12 Mev and e=2.18 Mev. 
E is the neutron energy. [Rev. Mod. Phys. 8, 117 (1936). ] 
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Fic. 4. Similar to Fig. 3 except 
that these data are only from the 


plates exposed to the LiCl target. 
This shows that the neutron group 
at 5 Mev of Fig. 3 was not due to 
the disintegration of the nitrogen 
of the Li;N which appeared to form 


NUMBER OF TRACKS 


on the surface of the metallic lithium 
pe et _ for the bulk of the data 
ig 


ENERGY IN MEV 


DIsCUSSION OF RESULTS 


Neutrons of 3.5 Mev are the highest energy 
ones that could be produced from the disin- 
tegration of the Li® isotope by 1200-kev deu- 
terons :” 


Mev (4) 
Be?+¢e,—Li?+7+0.87 Mev. 


Furthermore, reaction (3) can only give neu- 
trons of less than 4-Mev energy. Hence all the 
neutrons above 4 Mev must be from reaction (1) 
and from the following modes of disintegration 
of the Li’ isotope: 


Mev, (5) 
(6) 
*Be’—2 He‘+(Q;. (7) 


Reaction (5) is certainly not the favored mode 
of disintegration since a three-particle disintegra- 
tion cannot give the observed homogeneous 
neutron groups at 5 Mev, 7.5 Mev, 10.8 Mev, 
and 14 Mev (Figs. 3 and 4). Therefore these 
neutron groups are assigned to reaction (1) and 
to reaction (6) where Be® may be formed in 
various excited states. 

The rather sharp 14-Mev neutron group 
(first reported by Bonner and Brubaker’) then 
corresponds to the formation of Be® in the ground 


2 This Q value is calculated from the mass of Be’ given 


18 6. Shoupp, Stephens and Wells, Phys. Rev. 58, 


state which (according to Wheeler’s recent 
analysis of other experiments) is a 100-ev 
wide 'S») level about 125 kev unstable against 
alpha-particle emission. Hence the observed 
0.5-Mev half-width of the neutron group may be 
considered as experimental" rather than real 
width. 

The wide neutron group at about 10.8 Mev 
is similar to that reported by Stephens? from his 
study of recoil alpha-particles in a cloud chamber. 
This equivalence of the alpha-particle and proton 
recoil data shows that there are no excited levels 
of He' in this region which give resonance neu- 
tron scattering. 

This neutron group at 10.8 Mev has previously® 
been identified with the 2.8-Mev level in Be® 
which has been observed in boron disintegra- 
tions.'® However, the observed half-width of this 
neutron group (Figs. 3 and 4) is at least 2 Mev. 
On the assumption that the half-width of the 
14-Mev group is experimental, then the natural 
width of the Be® level associated with the 10.8- 
Mev neutron group would be at least'® 


(9/8)(2.0—0.5) =1.7 Mev. 


18 J. A. Wheeler, Phys. Rev. 59, 27 (1941). 

4A recent determination in this laboratory of the ex- 
citation function for this disintegration indicates that the 
use of a thick target should in itself produce a 0.3-Mev 
experimental width to the neutron group. 

% Oliphant, Kempton and Rutherford, Proc. Roy. Soc. 
150, 241 (1935); P. I. Dee and C. W. Gilbert, Proc. Roy. 
Soc. 154, 279 (1936); J. D. Cockcroft and W. B. Lewis, Proc. 
Roy. Soc. 154, 246 (1936 ). 

6 50= =(M/M,)8E:’, from (7706), page 277 of reference 9. 
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Fic. 5. The d-d neutron 
spectrum used for the calibra- 
tion of the stopping power of the 
plates. The plates were placed 


at 90° and 0° to a 700 kv bom- 60 
barding beam. Since a_ thin 

target was used, these data also § 
give an indication of the re- £ 
solving power of the emulsion 4 ,, 


method. (For interpretation of 
the high energy tail on the 90° 
data and the group at 2.5 Mev 
on the 0° data see Richards and 


Hudspeth, Phys. Rev. 58, 382 20 
(1940).) 


This is over twice the 0.8-Mev width of the 
2.8-Mev level of Be® which is observed in boron 
disintegrations.'® Furthermore, the peak of the 
recoil protons from neutrons associated with a 
2.8-Mev level in Be® is calculated to come at 
11.6—(4)E,0°=11.4 Mev.® This is about 0.6 
Mev higher than the 10.8-Mev observed peak 
(Figs. 3 and 4). Thus it would seem that not all 
of the neutrons in this group correspond to the 
2.8-Mev level in Be*. The existence of another 
level in Be® at about 4 to 5 Mev would account 
for the remaining neutrons in this group and also 
explain its apparent width. Independent evidence 
for such a level is given by the recent discovery 
in this laboratory" of a fairly intense 4.9-Mev 
gamma-ray from deuteron bombardment of 
Li’. This gamma-ray follows the same excitation 
function as the neutrons and so it may be sup- 
posed to result from a 4.9-Mev excitation level 
which has odd angular momentum and/or odd 
parity so that the ordinary breakup into two 
alpha-particles is prohibited. A selection rule 
depending on this odd angular momentum 
and/or parity might then explain why this level 
is not observed in the boron disintegrations 
when the intermediate C" nucleus breaks up into 
Be® and He‘. 

The next lowest neutron group (~7.5 Mev) 
has a real half-width of probably about 1.5 Mev. 


Hence it requires a 7.5-Mev excitation level in 


* Bennett, Bonner, Richards and Watt, Letter to the 
Editor, Phys. Rev., to be published. 
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Be® which would be about 1.7 Mev wide. 
From the alpha-particles of the B'°+H? dis- 
integration Smith and Murrell'® infer a similar 
level in Be® although perhaps not quite so wide. 

The neutron group at 5 Mev (only partially 
resolved) must be associated with a wide 10-Mev 
level in Be*. Estimates of its width are quite 
uncertain since it is only partially resolved. 
However, it appears to be about as wide as the 
7.5-Mev group. This group at 5 Mev is not 
inconsistent with Staub and Stephens’ neutron 
distribution® from alpha-particle recoils although 
they drew their published curve through the 
proton-recoil points because of the uncertainty 
in the neutron alpha-particle scattering cross 
section. Since the present data remove this 
uncertainty it might be expected that their data 
with the He recoils would be the more reliable 
since their geometrical correction is quite small 
for the He recoils in this region. 

The present data show no plateau in the region 
0.8-3.9 Mev which could be ascribed to neutrons 
from the two-stage reaction (3) though part of 
the sharp rise at about 4 Mev may be from this 
source. However Rumbaugh, Roberts, and 
Hafstad'® found that for 800-kev deuterons 
about 8 percent of all neutrons are from the 
Li® isotope. Since these are all confined to the 


18 C, L. Smith and E. B. M. Murrell, Proc. Camb. Phil. 
Soc. 35, 298 (1939). 
a — Roberts, and Hafstad, Phys. Rev. 54, 657 
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lower energy portion of the curve, we may ex- 
pect that about a fourth of the neutrons with 
energy less than 3.5 Mev are from the Li® 
isotope. Hence, unless reaction (3) is the pre- 
dominantly favored mode of disintegration in 
this region, the neutrons from Li® and reaction 
(6) would be sufficient to obscure such a plateau. 
This would certainly be the case if the relative 
frequency of reaction (3) to reactions (5)+(7) 
were about 1 : 100 as Williams, Shepherd and 
Haxby”® report from their studies of the alpha- 
particles from the above reactions. 

It is uncertain whether the peak at 3 Mev is 
from the formation of a still more highly excited 
Be® or whether it is from the Li® disintegra- 
tion in which Be’ is formed, reaction (4). 
Other modes of disintegration of the Li® isotope 


20 Williams, Shepherd and Haxby, Phys. Rev. 52, 390 
(1937). 
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are: 

Li®+ H*—He'+ Mev 
Li®+H*?—He*+ He®+0.8 Mev 

He®—He!+n+0.8 Mev. 


For 1200-kev deuterons neither of these re- 
actions yields neutrons of over 2 Mev. Hence 
curve JI, Fig. 4, indicates that these are not 
prolific reactions. 

Nevertheless, the interpretation of the spec- 
trum below 3.5 Mev remains ambiguous until 
a careful study is made of the neutrons from a 
target of the separated Li® isotope.t 
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7 Rumbaugh, Roberts and Hafstad, reference 19, have 
made a rough measurement of this neutron distribution. 
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Deuteron Disintegration by Electrons 


B. Peters AND C. RICHMAN 
Department of Physics, University of California, Berkeley, California 
(Received March 14, 1941) 


The cross section for the disintegration of deuterons by electrons has been calculated. Near 
the threshofd the ‘‘magnetic’’ cross section is larger, increasing with the square while the 
“electric’’ cross section increases with the cube of the energy. They become equal at about 
600 kev above the threshold when each has a value of 3.5X10-*! cm*. The disintegration 
neutrons should be observable from 50 kev up. An experiment of this kind might therefore be 
used to check the theoretically predicted magnetic disintegration. This seems desirable since a 
parallel effect for photo-disintegration has so far not been observed. 


HE electron currents and voltages which 
have recently been obtained! with electro- 
static generators make it appear probable that 
the disintegration of deuterons by electrons can 
be observed. For the present such experiments 
will have to be confined to electron energies only 
slightly in excess of the deuteron binding energy. 
As in photo-disintegration, the only appreci- 
able transition probabilities are those to the *P 
and 'S states of the disintegrated deuteron. The 
corresponding cross sections are usually referred 
to as “electric’’ and “magnetic,” respectively. 
1D. L. Northrup, L. C. van Atta, R. J. Van de Graaf, 


J. S. Clark and C. M. van Atta, Phys. Rev. 57, 563 (1940) 
and Phys. Rev. 58, 199 (1940). 


The theory predicts, both for photo- and elec- 
tron-disintegration, a predominance of the mag- 
netic transition and corresponding spherically 
symmetric distribution of the emitted particles, 
for energies just above the threshold. Protons and 
neutrons ejected in the direction of the incident 
gamma-ray have been looked for but have not 
been found. Measurements were made with 2.62- 
Mev gamma-rays by J. Chadwick, N. Feather 
and E. Bretscher? and also by H. v. Halban.* No 
evidence for magnetic disintegration was found. 
V. Halban puts the ratio of magnetic to electric 

2 J. Chadwick, N. Feather and E. Bretscher, Proc. Roy. 


Soc. A163, 366 (1937). 
3H. v. Halban, Nature 141, 644 (1938). 
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cross section at this energy at less than 6 percent. 
It seems, therefore, worth while to investigate 
whether a magnetic effect in agreement with 
theory can be observed if electrons instead of 
gamma-rays are used for the disintegration. 

In calculating the magnetic cross section we 
employed the same method as that used by 
Bethe and Peierls‘ in their calculation of the 
electric cross section. 

The perturbing term in the Hamiltonian is 
given by 

V=—(eh/2Mc)(uet+y'o’)-curl A, 
where .V/ is the mass and uy, y’, @, o’ are the mag- 
netic moments and spin operators of the neutron 


and proton respectively. 
For A we use the Mller potential® 


i 
4rh*ea’* aa o— 


hen 
Ap*—(AE/c)? 
Here a is the Dirac current operator, a’ and a are 
the Dirac amplitudes of the wave functions 
describing the incident and scattered electron. 
Ap=po—p’ and AE=E,—E’ are the change in 
momentum and energy suffered by the electrons 
in the collision process. 
In the approximation of zero range for the 
neutron-proton force, the ground state of the 


deuteron can be described by 


where 7 is the distance between neutron and 
proton, x3 the spin function of the triplet state, 
1/a=(h*?/ Me)! the “radius” of the deuteron and 
« the binding energy of the ground state. 

For energies up to 25 Mev the wave-lengths 
of the electrons are much larger than the size 
of the deuteron and therefore Ap/h<a. It will 
then be sufficient to use the first term in the 
power series expansion of the exponential in A. 
The final state of the deuteron will then also be 
an S state. We shall describe it by 


ve= [sin(p/hr+ 5)x}/r, 


where p=(ME)! is the momentum of the 
disintegrated deuteron, E= Ey—«—E’, its energy 
and x, the spin function of the singlet state. 
For final states, whose energy is small compared 
to the depth of the potential well of the neutron- 
proton force, the logarithmic derivative at the 
edge of the well will not depend on energy. Its 
value at zero energy is known from scattering 
experiments. In the approximation of zero range 
we obtain therefore coté=(y,£)'. The experi- 
mental value for vy is 80 kev. 

The cross section for the transition is given 
by the well-known formula 


do, =(22/hvo) | Vor! *pp'dE'dQ’ 


where v)=c*po/Eo is the velocity of the incident electron, p’ = p'E’/82°*h'c? is the density of final 
electron states in the energy range dE’ and the solid angle dQ’. p=(1/42*h)(M/E)} is the density of 
final deuteron states in the energy range dE=d(E,)—e—E’). The evaluation of the matrix element 


gives 
Von |2=32 (-) al(p/h)yi+akE!} | }x:) 
(a?+ p?/h®)?(y+E) [Ap?—(AE/c)?}? 


We now average over the three orientations of the magnetic moment of the deuteron in its ground 
state and the two spin orientations of the incident electron, and sum over the two spin orientations 
of the scattered electron. Expressing the energy in units of mc? and the momenta in units of mc, we 
obtain 


2 

3a a? 1 
po +p” AE* 
pol4 2Ap?—AE? 


*H. A. Bethe and R. Peierls, Proc. Roy. Soc. A148, 146 (1935). 


*C. Moller, Zeit. f. Physik 70, 786 (1931), 
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Integration over the scattering angles gives: 


2 dE'(E,—«—E’)! E.E’+pop’-1 
3 he a? 1 E,—E’ 


The electric cross section is obtained by a parallel calculation. Here the perturbing potential is 
given by V=eo—(e/c)A- (h/2iM)grad,. The final deuteron state which is a P state can in zero-range 
approximation be described by the P component of a free wave ~z=exp[(i/hk)(p-r) ]. The cross 
section so obtained is given by 


This result is not identical with the formula given by Bethe and Peierls.‘ Their expression, which we 
believe to be in error, is twice as large as ours and its last term is 3p’/ po instead of 2p’/po. 


(2) 


RANGE CORRECTION 


If the finite range of the neutron-proton force 
is taken into account, the electric cross section 
has to be multiplied by a factor which varies 
very slowly with energy and depends somewhat 
on the type of interaction assumed. These factors 
have been tabulated by G. Breit and E. U. 
Condon.* For a range a=2.8X10-" cm the 
factor decreases approximately linearly from 
1.95 at the threshold to 1.85 at 2.5 Mev above 
the threshold, if one assumes Majorana inter- 
action; and increases from 2 to 2.25, if one 
assumes ordinary forces. We write therefore for 
the range correction : 


B(Majorana) = 1.95 —0.02(W,— 
B(ordinary) = 2+0.05(Wy—e), 


where Wo=Eo—1 is the kinetic energy of the 
incident electron in units of mc’. ; 
The range correction is very much smaller in 
the magnetic case. Since the final wave function 
is large near the origin, a decrease of the ground 
state wave function yp» at small distances will 


- reduce the value of the matrix element. This 


decrease will largely offset the increase in the 
normalization constant of Wo. For a square hole 
potential the factor 8, is 1.04 near the threshold. 


vir se?\2 1 (e+2) 
6 a? e(e+2)! 


For high energies both formulae take the form 
o=« In (E¢?/se), 


where «x is a constant and s is a factor of order 
unity. They are identical with the cross sections 
which one can derive by the method of virtual 
quanta, if one uses c£/v for the upper and h/mcé 
for the lower impact parameter (= (1—v?/c*)-}). 
The formulae which we calculated can therefore 
no longer be applied if h/mct=hc/Ey<1/a, that 
is when the electron passes through the deuteron. 
If the condition E»<hca is not satisfied, then 
impacts are possible for which the momentum 
transfer Ap is larger than ha. As we have men- 
tioned before, the expansion of the potentials 
in a power series and therefore the neglect of 
the variation of the field across the deuteron is 
then no longer justified. Comparison with the 
formula derived by the method of virtual quanta 
shows that the improper inclusion of these im- 
pacts of large momentum transfer introduces an 
error into our formula, although the energy trans- 
fer in the disintegrations is never large and 
always of order e. 

At low energies, up to about 250 kev above 
the threshold the formulae can be written 


Wo-e Wo-«\} 


€ 2y 


6G. Breit and E. U. Condon, Phys. Rev. 49, 904 (1936) and Phys. Rev. 51, 56 (1937). 
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TABLE I. Cross section for disintegration of deuterons 


by electrons. 
ENERGY IN 
EXCESS OF 
BINDING 
ENERGY 50 KEV 0.5 MEV 0.75 Mev 2.5 MEV 


(cm*) 3.5X10-" 2.0X10-" 5.2x10-" 
cm? (2.2 10-8) (5.6X10-%) 7.2X10-% 


Total cross section 
(cm*) 8.3X10-" 4.5X10-"% 
(4.7K (9.7 107%) (8.5 x 107%) 


Putting in the numerical values for the con- 
stants, we obtain: 


T 
Tota! CS 
% / 
7 
z / SMagnetic CS 
23 
‘4 a 
A 
bo 
& Electric CS (Mojorana) 
ENERGY THE THRESHOLD 


Fic. 1. Electric, magnetic, and total cross section for the 
disintegration of deuterons by electrons. 


= 2.8 X 10-*°[ Mev cm’, (3) 
Mev—{1+12.5(Wo—e) Mev}*] cm*. (4) 


We have used (3) and (4) to calculate the 
cross sections at 50 kev and have evaluated (1) 
and (2) numerically for 0.5 Mev, 0.75 Mev and 
2.5 Mev. The results are given in Table I (the 
figures in parenthesis refer to ordinary, the others 
to Majorana-type interaction) and in the graph 
of Fig. 1. 

With an electron current density of 10 micro- 
amperes, a cross section of 10-* cm? corresponds 


to the production of about 10‘ neutrons per 
minute per milligram of deuterium per cm?. Such 
a cross section should therefore still be observ- 
able. It corresponds to an energy of about 50 
kev. 

We want to express our thanks to Dr. J. 
Schwinger and to Professor J. R. Oppenheimer 
for having suggested the problem and for many 
helpful discussions. 
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1941 PHYSICAL REVIEW 
Zeeman Effect of the Forbidden Lines of Pb I 


I. Pure Quadrupole and Magnetic Dipole Lines 


F. A. JENKINs AND S. MrozowskI1 
Department of Physics, University of California, Berkeley, California 
(Received March 17, 1941) 


A very intense source of the forbidden lines, consisting of a high frequency discharge through 
helium containing a small amount of lead vapor, is used to photograph the Zeeman patterns in 
fields up to 4000 gauss. Reproductions are given of the transverse effect for the main hyperfine 
component, due to the even isotopes, in a magnetic dipole line, 44618, and in two quadrupole 
lines, \A4659 and 5313. The reversal of polarization in 44618 is clearly shown. The splittings 
of the weaker hyperfine components emitted by Pb?’ are also observed for \A5313 and 4618, 


_ including, in the latter case, fields strong enough to show the Back-Goudsmit effect. Splittings 


and intensities are in satisfactory agreement with theory throughout. These are the first 
observations on the Zeeman effect of the hyperfine structure of forbidden lines, and they 
establish the applicability of the electric dipole intensity formulas (with interchange of polariza- 
tion) to magnetic dipole lines and of the Rubinowicz quadrupole formulas (with suitable 
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substitutions of quantum numbers) to quadrupole lines. 


N a previous paper,! one of us has described a 

very intense source of the forbidden lines of 
lead that are emitted in transitions between the 
levels of the normal (6*) configuration. It was 
possible to investigate the hyperfine structure 
due to Pb?’ in the case of three of these lines, 
and to measure the relative intensities of the 
five lines observed. According to the selection 
rules, the strongest line, (4618 ('S)—*P,) is 
permitted only for magnetic dipole radiation, two 
others, 45313 and 44659 
only for electric quadrupole, while the remaining 
ones, \7330 (1D.—*P;) and 9250 (‘D.—*P») in- 
volve both types of radiation. Since the existing 
data on the Zeeman effect of forbidden lines are 
rather meager, and in no case include the hyper- 
fine structure, it seemed worth while to in- 
vestigate the effect in these lines. 

The Zeeman effect of 44618 has been studied 
by Niewodniczanski? and this constitutes the 
only existing example of the effect in a pure 
magnetic dipole line. With our improved source 
we obtained the line with an intensity at least 
500 times greater than that attained by Niewod- 
niczanski. This made it possible to observe the 
Zeeman effect not only for the strong central 
component due to the even isotopes of masses 


1S. Mrozowski, Phys. Rev. 50, 1086 (1940). 
2H. Niewodniczanski, Acta Phys. Polonica 3, 285 


(1934). 


204, 206 and 208, but also for the much weaker 
hyperfine components due to the odd isotope, 
207. The same was true for the pure quadrupole 
line 45313, even though this is some five times 
weaker than \4618. For the other, still weaker 
quadrupole line, \4659, only the central com- 
ponent could be studied. Of the two mixed lines, 
the stronger one at \7330 has been investigated 
in detail, including quantitative intensity meas- 
urements of the Zeeman components. In Part I 
of this paper we shall give the results on the 
pure quadrupole and magnetic dipole lines, 
while the case of the mixed line, in which a new 
interference effect between the two types of 
radiation has been discovered, will be considered 
later in Part II. 


EXPERIMENTAL 


The lines were excited by a high frequency 
discharge in helium containing a small amount of 
lead vapor. The characteristics of this source 
have already been described in some detail,' 
so it will only be necessary to give here the 
modifications which became necessary for its 
operation in a magnetic field. Figure 1 shows, 
to scale, the cylindrical quartz discharge tube T 
and the asbestos-insulated Nichrome oven F 
situated between the poles P, P of the magnet. 
For large fields it is necessary to have the dis- 
charge current in the direction of the field, so in 
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this case the high frequency potential was ap- 
plied to external electrodes E, E consisting of 
strips of platinum 5 mm wide placed above and 
below the tube. This disposition of the elec- 
trodes has the advantage that the deposit of 
carbon which gradually forms near the elec- 
trodes! does not occur on the window. The 
intensity with which the forbidden lines can be 
produced is, however, much less than with the 
previous arrangement using ring-shaped elec- 
trodes near the ends of the tube. The loss of 
intensity is due to the fact that when one ap- 
proaches the best conditions for the forbidden 
lines (high helium pressure, or high vapor pres- 
sure of lead), the bright part of the discharge 
splits in a strong magnetic field into two regions 
close to the electrodes, leaving the middle with 
fairly low luminosity. This is no doubt caused 
by the spiraling of the electrons around the 
magnetic lines of force, which leads to an increase 
of ionization near the electrodes, and a more 
intense emission of the allowed lines in these 
regions. The forbidden lines, however, are 
emitted with nearly uniform intensity across 
the space between the electrodes, as a result of 
the diffusion of metastable atoms. These lines 
could be easily distinguished visually from the 
allowed lines by this difference in intensity 
distribution. The best conditions for the lines 
investigated here were a helium pressure of 4 
to 5 mm and a temperature somewhat below 
800°C, such that the PbII lines are still very 
bright and the He lines have not yet disappeared 
completely. 

With the smaller magnetic fields used in 
investigating the Zeeman effect of the hyperfine 
structure, the ring-type electrodes could be used, 
and gave an intensity of the forbidden lines even 
greater than that obtained without magnetic 
field. At a field greater than about 1000 gauss, 
the discharge concentrates at each end toward 
certain points inside the electrodes. The result- 
ing increase in current density destroys the 
metastable atoms, and although the intensity 
of all allowed lines increases with further increase 
of the magnetic field, the forbidden lines become 
weaker. 

The magnet used’ could give a field of 4000 


7A. V. Hershey, Phys. Rev. 56, 908 (1939). Thesis, 


University of California, 1939. 


gauss constant to within one percent over the 
space occupied by the tube. Since the oven had to 
be maintained at a high temperature over long 
periods of time,’ water-cooled copper plates 
were placed over the pole faces to prevent over- 
heating of the magnet when it was not in opera- 
tion. The optical apparatus, consisting of a 
one-prism glass spectrograph and silvered quartz 
Fabry-Perot plates, were the same as those 
used in the previous work. 


10 cm 


Fic. 1. Arrangement of the source. 


Measurements of the Zeeman splittings were 
made on a comparator, including several orders 
of each Fabry-Perot pattern. By interpolation on 
a parabolic curve, linear distances were ex- 
pressed as fractions of the separation of orders, 
and for each component the average fraction 
from the different orders was then converted to 
wave numbers. 

For the calibration of the magnetic field, the 
Zeeman splitting of the singlet line \4916 
(61P,:—8 Sio) of mercury was carefully meas- 
ured at six different values of the field. The 
mercury spectrum was excited with the same 
apparatus by merely introducing mercury vapor 
into the tube. With Av=4.67X10-> Hg, and 
with g=1.010 for this line,* values of H were 
found which were accurately proportional to the 
current in the magnet. 


RESULTS 


Zeeman effect of the main components 


From both the experimental and theoretical 
sides, the investigation of the Zeeman effect 
of the main component due to the even isotopes 
is much simpler. As was shown by Mrozowski,! 


* P. Gmelin, Ann. d. Physik 28, 1079 (1909). 
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this component contains 79 percent of the 
intensity of each line, in agreement with the 
abundance of the isotopes. Its Zeeman pattern 
could be photographed with relatively short 
exposures (30 sec. in the case of \4618) while the 
pattern for the hyperfine components emitted 
by the odd isotope required more than ten times 
this exposure, and was partly overlapped by the 
main pattern. Because of this overlapping, no 
quantitative intensity measurements were at- 
tempted for the patterns reported here, but a 
qualitative comparison with theory gives an 
important conclusion mentioned below. 

The fourth column of Table I gives the results 
of measurements on the main patterns. The 
magnetic dipole line (4618 gives a }$-normal 
triplet, but with the center component polarized 
perpendicular to the field, the outer ones parallel. 
This reversal of the normal polarization is shown 
very clearly in Fig. 2(b), where the line can at 
once be distinguished from electric dipole lines 
by this characteristic. This photograph was 
obtained by consecutive illumination of the two 
halves of the slit of the spectrograph with light 
of different polarization. The exposure time in 
the x case was made twice as long as in the 
other in order to equalize the blackenings in the 
line 44618. This explains why the 7 components 
of the neighboring electric dipole lines appear so 
strong in comparison to the « components. The 
« component of 44618 nearly overlaps with a 
= component of a neighboring electric dipole line. 
The reversal of the normal polarization has 
already been observed by Niewodniczanski,’? but 
the reproduction of the line in his article was not 
satisfactory. According to the theory® the mag- 
netic dipole radiation field is obtained from the 
electric dipole case by substituting H for E and 
—E for H and by replacing the electric moment 
by the magnetic moment. Hence the selection 
and intensity rules for the Zeeman effect of 


TABLE I. Zeeman patterns of the main components. 


TRANSITION RADIATION Oss. CALc. 


4659.18 4So—'P, 0 (1.48) 0 (1.501) 
H 


pole 
5313.00 (1.23) 2.53 (1.269) 2.538 
4659.02 'D:—*P, Quadrupole (1.20) 2.45 (1.23) 2.46 


5A. Sommerfeld, Atombau und Spektrallinien (1939), 
Vol. 2, p. 735. 
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magnetic dipole lines are the same as for electric 
dipole lines, but with the polarizations inter- 
changed. Thus our result gives convincing 
evidence for the magnetic dipole character of 
this line. 

The two quadrupole lines give a characteristic 
pattern of four lines of equal intensity (Fig. 2(d) 
and (f)); corresponding to AM=+1(mr) and 
+2(c). The transition A/=0 is also allowed, 
but has zero intensity for either transverse or 
longitudinal observation, appearing only as an 
elliptically polarized component in diagonal ob- 
servation.* The calculated splittings given in 
the last column of Table I were obtained by using 
the g values determined by Back’ from the 
allowed lines involving these levels, namely for 
and 0/0, for *P;, 1.501, for *P»2, 1.269, 
and for ‘Dz, 1.230. It will be seen that our ob- 
served splittings are consistently a little lower 
than the computed ones. This may possibly be 
due to the fact that the g value we adopted for 
the calibrating Hg line \4916 was a little too 
small. Our measurements were quite accurate 
for the lines \4618 and \5313, but somewhat 
less reliable for \4659, because here the com- 
ponents were not as sharp and were partly over- 
lapped by neighboring lines. 


Zeeman effect of the hyperfine structure 


We have been able for the first time to ob- 
serve the effect of a magnetic field on the hyper- 
fine structure of forbidden lines. The two cases 
4618 (magnetic dipole) and 45313 (quadrupole) 
were investigated. Reproductions of the spectro- 
grams appear in Fig. 2(c) and (e), and the re- 
sults of the measurements are presented graph- 
ically in Fig. 3. For \4618, because of its greater 
intensity and smaller hyperfine splitting, ob- 
servations could be carried to large enough fields 
to show the Back-Goudsmit (‘strong field’) 
effect, while for 45313 only the ‘‘weak field” 
effect could be studied. 

In obtaining the theoretical patterns in Fig. 3, 
the splitting factors of the levels for weak fields 


as 3 Segré and C. J. Bakker, Zeits. f. Physik 72, 724 
7 E. Back, Zeits. f. Physik 37, 193 (1926). 
8On the Hg spectrograms the singlet line 5790 
(6!P,;—7'D:2) was also measured, and showed a splitting 
8 percent greater than \4916. This does not agree with 
Gmelin,‘ who concluded that \5790 gave a normal triplet. 
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ZEEMAN EFFECT IN LEAD sil 


A 4618 


Fic. 2. Zeeman patterns of the forbidden 
lines. (a) 4618. J7=2875 gauss; 5-mm = 
spacer; exposures 50 sec. (b) 44618 and ad- @e 
jacent electric dipole lines. H = 2890 gauss; 
3.5-mm spacer; exposures 40 min. and 20 e 
min. (c) 44618, components only, showing 
hyperfine structure. 7=0,670, 1090 and 
1590 gauss; 15-mm spacer; exposure 4+ min. = 
for H=0, 15 min. for the others. (d) 44659, — 
x components. The arrow indicates the 
pattern for this weak line, which is over- 
lapped by strong adjacent lines. Exposure 
45 min. (ec) 45313, ¢ components showing 
hyperfine structure H=0, 225, 500, 765 
gauss; 9-mm_ spacer; exposures 1 hour. 
30 min. (f) A5313, 2 and @ components. 


40 min. 


A 4659 


were obtained from the Goudsmit formula 
+1) 
2F(F+1) 


For strong fields, the splitting due to the nu- 
clear moment is given by the coup- 
ling constant A being determined from the 
hyperfine splittings of the levels without field, 
namely ]. 
(The over-all separation of the hyperfine struc- 
ture, Av= Finox( t+ 1) Foin( Fmin+ 1) ] 
=AI(2J+1), if J>J.) The splitting in inter- 
mediate fields could have been calculated ex- 
actly from the theoretical equations. However, 
since only the weak field effect was observed for 
\5313, and since for \4618 the hyperfine splitting 
is small and not known very precisely, this cal- 
culation was not made. Instead, a graphical 
interpolation between the weak and strong field 
levels was used, since this should give the transi- 
tion effect within the accuracy justified by our 
measurements. 

The upper diagram of Fig. 3 shows the pre- 
dicted and observed splittings for the magnetic 
dipole line, 44618, including both polarizations. 
The upper level, 'So, does not split, so that the 
diagram gives directly the splitting of the lower 


H=2890 gauss; 3.5-mm spacer; exposures (d) ea (e) 


A 4618 


(c) 
= 


= 
§ 
— 
—_ 
= 
re A533 


state,*P;. In the theoretical weak field patterns 
of the hyperfine structure components, for which 
the separations and intensities are indicated at 
the left of the figure, it will be seen that all 
components except the outer # components in 
the high frequency member of the hyperfine 
doublet should appear in both polarizations. 
The agreement between experiment and theory 
is complete, both as to separations and intensi- 
ties. Since the theoretical intensities are those 
calculated from the ordinary electric dipole 
formulas, we see that one can extend to the 
hyperfine structure the rule that everything is 
the same as in the electric dipole case, but with 
reversed polarizations. At the highest fields 
reached, the Back-Goudsmit effect is almost 
complete, and here also the observations agree 
satisfactorily, though they are rather incomplete. 

The lower diagrams of Fig. 3 show the results 
for the quadrupole line 45313. In this case also 
(\Sy—*P2), the splittings are those of the lower 
state. Because the pattern is more complicated, 
the and components are plotted separately 
in parts (a) and (0b) of the figure. The theoretical 
intensities for weak fields were obtained from the 
Rubinowicz formulas* by substituting F for J. 


*A. Rubinowicz and J. Blaton, Ergeb. d. exakt, 
Naturwiss. 11, 176 (1932). 
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Fic. 3. Predicted and 
observed splittings for the 
magnetic dipole line \4618 


Since the hyperfine structure is considerably 
wider in this case, we could not reach fields large 
enough to show the Back-Goudsmit effect. This 
line is much weaker and in order to get the hyper- 
fine structure components on not too long 
exposed plates the ring-type electrodes had to be 
used. Within the range of magnetic fields used, 
the observed pattern agrees well with that pre- 
dicted for weak fields. The applicability of the 
modified Rubinowicz formulas to the Zeeman 
effect in the hyperfine structure of a quadrupole 
line is thus verified for the first time. Not only do 


(upper diagram) and for 
the quadrupole line \5313 
(lower diagram). In both 
cases since the upper 
level 1So does not split, 
the diagrams give di- 
rectly the splitting of 
the lower states. The 
heavy lines refer to the 
main component; the 
lighter ones to the hyper- 
fine structure. Broken 
lines are the parallel com- 
ponents; the solid lines 
the perpendicular ones, 
Squares represent com- 
ponents observed in 7 
polarization; circles those 
in oa. Open squares or 
circles are for lines ob- 
served as completely re- 
solved from the main 
pattern while filled ones 
are for those which are 
in contact with or partly 
overlapped by the main 
components. 


we observe the lines AM = +2 in the perpendicu- 
lar effect, and AM= +1 in parallel, but also the 
intensities are qualitatively right for the quad- 
rupole case. For example, the highest frequency 
component is considerably weaker than the next 
one (as indicated by the smaller open square in 
Fig. 3(a)), corresponding to the predicted ratio 
1 to 3. The intensity formulas for electric dipole 
give just the reverse ratio, 3 to 1, for these two 
lines. 

The only feature of our pattern which does not 
agree with the theory for weak fields is the in- 
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tensity of the lowest frequency x component 
(smaller filled square in Fig. 3(a)), which is 
apparently much weaker than it should be. 
This may possibly be explained as an incipient 
Back-Goudsmit effect, since the line goes to zero 
intensity in strong fields and the effect is ob- 
served at the highest field we used for this line. 


The apparent weakness of this component is 
probably enhanced by the fact that the next 
component lies on the wing of a main line, and 
therefore looks relatively too strong. In any 
case, there is no reason here to doubt the validity 
of the existing theory of the Zeeman effect of 
quadrupole lines. 
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The Dynamics of a Particle in General Relativity 


C. Lanczos 
Purdue. University, Lafayette, Indiana 
(Received September 12, 1940) 


Some applications of the Gaussian integral transforma- 
tion to the conservation law of the matter tensor yield the 
complete dynamical law for a material particle in a Rieman- 
nian space, in the form of two relations which correspond 
completely to Newton’s two fundamental laws. The 
“moving force’ of the motion law comes out as a volume 
integral which can be transformed into a surface integral. 
The inertial mass and the coordinates of the center of mass, 
on the other hand, are genuine volume integrals which 
cannot be transformed into boundary integrals. For weak 
and static matter the principle of the geodesic line can be 
established, irrespective of spherical symmetry. But even 
for static and infinitesimal matter an interesting exception 
occurs if the total mass of the particle happens to drop 


I. INTRODUCTION 


VER since the discovery of general rela- 
tivity the question of the motion of a 
particle in a Riemannian field has been raised 
and discussed. The original paper of Einstein! 


_went a long way toward the establishment of a 


field dynamics, by showing that the matter 
tensor of general relativity satisfies the conserva- 
tion laws of energy and momentum as a mathe- 
matical necessity. He also succeeded in bringing 
that conservation law rigorously into the form 
of a classical divergence equation, by adding to 
the matter tensor certain auxiliary quantities 
which had no invariant significance but were 
interpretable as the momentum-stress-tensor of 
the gravitational field.” 

However, the conservation law of momentum 


1 A. Einstein, Ann. d. Physik 49, 769 (1916). 


? Cf. reference 1, p. 809. 


down to a quantity of second order, as is indeed the case 
if the scalar condition T =0 is satisfied inside of the particle. 
Then the acceleration is smaller than it should be accord- 
ing to the geodesic principle which has the consequence 
that the gravitational mass of the sun, determined in the 
customary manner by the acceleration of its planets, may 
be underestimated. This gives a new clue for the under- 
standing of the anomalous value of the light deflection 
on the limb of the sun, discovered by Freundlich and his 
collaborators. However, the theoretical prediction of an 
increase of 83 percent, deduced for a static and spherically 
symmetric particle, is too large in comparison with 
the 25-40 percent increase, indicated by Freundlich’s 
measurements. 


and energy in itself is not able to establish a 
motion law for the material particle. The 
mechanics of discrete points is based on the 
two fundamental announcements of Newton: 
(a) momentum equals mass times velocity; (b) 
the time rate of change of momentum equals 
moving force. If we have to reinterpret these 
laws in terms of field quantities, the second law 
offers no difficulties. The conservation law of 
momentum can be considered as the direct 
counterpart of the second law of Newton. The 
real difficulty lies in the first law, viz. in finding 
a kinematic interpretation for the momentum 
of the field. 

In an earlier attempt*® the author was able to 
solve the motion problem of a particle for the 
case of special relativity. The method is based 
on the application of the Gaussian integral 


3C. Lanczos, Zeits. f. Physik 59, 514 (1930). 
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transformation. The momentum of the particle 
can be proved to be equal to the product of the 
total mass times the velocity of the center of 
mass. H. P. Robertson‘ extended the method 
later under certain restricting conditions, to the 
case of general relativity. 

Einstein, Infeld and Hoffman® attacked the 
motion problem from a different angle. Based 
on the consideration that the moving force 
acting on the particle must have its origin in the 
external field, and cannot depend on the matter 
distribution inside of the particle, they derive 
the moving force by discarding the central field 
entirely. They apply the method of successive 
approximations, combined with expanding the 
outside field in negative powers of the distance. 

The difficulty of the author’s own method lies 
in the fact that the moving force comes out as a 
volume integral, extended over the matter- 
occupied central field of the particle. The un- 
satisfactory feature of this result has been 
pointed out different times by Einstein in 
private correspondence with the author. The 
author is greatly indebted to him for his clarifying 
remarks. Indeed, let us assume a purely static 
particle. We may modify the metrics inside of 
that particle in any arbitrary static manner and 
still the field remains static and no motion occurs. 
Hence at least the time independent part of the 
moving force must be determined by the external 
field alone. 

At last the author succeeded in transforming 
the volume integral of the moving force into a 
boundary integral, extended over the border of 
the particle, or any closed surface of the external 
field which includes the particle. This has 
removed the last difficulty. The motion law is 
now rigorously established, without taking 
recourse to successive approximations. 

An interesting feature of the new motion law 
is that the right side of the law appears as a 
boundary integral, while on the left side the mass 
of inertia and the coordinates of the center of 
mass of the particle are determined by genuine 
volume integrals which cannot be transformed 
into surface integrals. The different nature of the 
two types of integrals leads to a startling new and 


a 937} P. Robertson, Proc. Edinburgh Math. Soc. 5, 63 
5A. Einstein, L. Infeld and B. Hoffman, Ann. Math. 39, 
65 (1938). 


unexpected consequence, even for a static and 
spherically symmetric particle, for which the 
principle of the geodesic line seemed to be firmly 
established by all previous considerations. 


2. THE Two NEwTONIAN Laws oF Motion, 
OBTAINED BY THE GAUSSIAN INTEGRAL 
TRANSFORMATION 


The detailed description of the method em- 
ployed by the author may be found in his earlier 
paper? and shall not be repeated here. Moreover, 
the same method has been utilized by the author 
in a very recent publication® in which he was 
able to give a theoretical explanation for the 
necessarily positive character of the mass. 
An abbreviation utilized in that paper shall be 
employed again for the process of integrating at 
a given time moment over the entire volume of 
the particle. We shall put: 


(2.1) 


On the boundary of the region of integration and 
outside of it the matter tensor 7 is supposed 
to be zero. 

We start out with the divergence equation of 
the matter tensor: 


—TyT”. (2.2) 


It is important to employ the purely contra- 
variant 7“ and not the currently used mixed 
T*, because for the kinematic interpretation 
of the momentum the symmetry property 
T*=T* of the matter tensor is absolutely 
essential. 

The following notation will be introduced. 
In these formulae the index 7 shall assume only 
the space values 1, 2, 3, and not the value 4: 


=p’, (2.3) 
(TH) = M, (2.4) 
tim THY, (2.5) 


Since. is the momentum density, 7% the 
energy density, we have to interpret p‘ as the 
total momentum, and WM as the total mass 
(energy) of the particle. The 3 quantities & 
represent the coordinates of the center of mass 


®C. Lanczos, Phys. Rev. 59, 708 (1941). 
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of the particle. In distinction to the mass of 
inertia Wf we will denote by m the gravitating 
mass of the particle, as it appears in the gq, of 
the external field. (In the previous paper® A 
was denoted by m.) 

Finally we put, without any restriction of the 
indices: 


i 


=F. (2.6) 


The application of the Gaussian integral 
transformation to the conservation law (2.2) 
yields at once: . 

dp'/dx4= F', (2.7) 
dM /dx4= (2.8) 


These equations show that the 4 volume integrals 
(2.6) have to be interpreted as the 4 contra- 
variant components of the moving force. The 
equation (2.7) expresses the second law of 
Newton: “The time rate of change of the 
momentum is equal to the moving force.” 
The equation (2.8) expresses the energy theorem. 

We now multiply the 4th equation of (2.2) 
by x; where again the index 7 shall not include 
the value 4. We may write 


and the Gaussian integral transformation yields: 
=(d/dxs) +{T'T, xi}. (2.10) 


In view of the notations (2.3), (2.4), and (2.5) 
we get: 


xi} (2.11) 
= 


Here we have now the kinematic interpreta- 
tion of the momentum. The momentum p‘ of 
the particle comes out as the product of ‘“‘mass 
times velocity,” in accordance with 
law. The ‘“‘velocity”’ has to be understood as the 
velocity of the center of mass. An additional 
term is present: 


pi= (2.12) 


but it is not difficult to see that this term must 
be exceedingly small, owing to the small values 
of For spherically symmetric particle 
p' vanishes by reasons of symmetry. Even for 
nonspherical particles p'‘ will vanish except if the 
particle has a stationary spin. But even then jp‘ 
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can only be a small correction term which can 
contribute only negligible amounts to the law of 
motion. Hence it seems justified for most 
purposes to neglect p', and the important 
kinematic relation results: 


pi'= Mdé'/dx,. (2.13) 


The kinematic condition (2.13), together with 
the motion law (2.7), establishes the entire dynamics 
of the particle, in full harmony with Newton's two 
laws of motion. In deducing these laws, no refer- 
ence was made to any specific structure of the 
material particle. The two fundamental rela- 
tions (2.11) and (2.7) are entirely general. 


3. THE MovinG Force As A SURFACE 
INTEGRAL 
According to (2.6) the moving force F' is 
determined by a volume integral, extended over 
the entire particle: 


(3.1) 


It shall now be shown that this integral is trans- 
formable into a surface integral because the 
quantity 


Ki=-T’T,,, (3.2) 
permits the following representation : 
(3.3) 


For this purpose we have to return to some 
ingenious relations which exist between the 
calculus of variation and general relativity. 
They have been discovered right at the start of 
the theory by Einstein,’ Hilbert® and Klein® 
and further developed by Weyl,'® Pauli," and 
others. These beautiful relations have escaped the 
attention of the author until recently because 
they operate with quantities which have no 
covariant significance and thus seem to con- 
tradict the basic spirit of relativity. Actually 
these relations are of fundamental importance 
and permit at once, with a slight modification 
of the procedure of Einstein, the transformation 
of the moving force into a boundary integral. 

We consider the function G which is the 
Lagrangean function of the Einsteinian field 


7 A. Einstein, Berliner Berichte, p. 1111 (1916). 
8D. Hilbert, Goettinger Nachrichten, p. 395 (1915). 
* F. Klein, Goettinger Nachrichten, p. 235 (1918). 
1H. Weyl, Ann. d. Physik 54, 117 (917). 

"W. Pauli, Physik. Zeits. 20, 25 (1919). 
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equations 
8 a a 8B 
(3.4) 


The invariant scalar curvature R is in the follow- 
ing relation to G: 


(3.5) 
where we have put: 


log 


OX 


(3.6) 


The variation of the invariant action integral 
yields: 


OX. 


-f ite) 


0 
veve) far (3.7) 
OX. 


with the notations: 
Lik, dr=dx}: -dx4. (3.8) 


The second and third terms are transformable 
into surface integrals but for our present pur- 
pose we leave them as volume integrals. 

We now specify the variation of the gx to the 
following particular type: 


Xk Ox; OXa 


The significance of this particular variation is 
that the points of the continuum are displaced 
by the infinitesimal vector ed‘ and thus distort 
the reference system, but without changing the 
geometry. A variation of the form (3.9) repre- 
sents merely an infinitesimal transformation of 
the coordinates. The vector \‘ is an arbitrary 
function of the 4 variables x4. 

A variation of this nature varies the action 
integral in a definite manner. We get: 


Hence we have two different expressions for the 
same integral: the right side of (3.7) and the 


(3.10) 


right side of (3.10). Equating these two expres- 
sions we get a relation which must hold inde- 
pendently of the region of integration. This is 
only possible if the integrands themselves are 
equal. We thus get: 


0 
OX, 
Xa 


This important relation holds for any choice of 
the vectorial function \‘. The previous literature 
dropped the terms with the V‘ and restricted the 
\‘-quantities to constants." But in the form 
(3.11) the relation holds for arbitrary X*. The 
explicit evaluation of the relation (3.11) is 
somewhat cumbersome under general conditions 
but a great simplification can be achieved if we 
make use of Einstein’s procedure! and consider 
a reference system in which 


g=1. (3.12) 


The fact that the right side is a perfect diver- 
gence, holds, of course, quite generally. The 
restriction of the reference system due to the 
condition (3.12) is merely a matter of con- 
venience because the general expressions do 
exist in any reference system, only their form is 
somewhat more complicated. 

The variation (3.9) of the gi, contains the first 
derivatives of the gx. In the vectorial quantity 
V‘ the gx are differentiated. The variation 6V‘ 
will thus contain the second derivatives of the 
gi, and the formula (3.11) shows that even 6V‘ 
is differentiated. Hence even third derivatives 
are present at the right side of (3.11). However, 
by combining the terms properly, the third 
derivatives all cancel out and the right side of 
(3.11) comes out finally for any arbitrary \‘ as 
follows: 

(a) Contribution of the G terms: 


| (- an}... (3.13) 


with the notation 
p'=0g'*/0Xa. 


(3.14) 


1 Cf. W. Pauli: ‘‘Relativitaetstheorie,"’ Ency. d. math. 


Wiss. (Teubner, Leipzig, 1921), Vol. V, 19, 621; H. Weyl; 
Space-Time- Matter (Methuen & Co., London), p. 269. 
13 Cf. reference 1, p. 801. 
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(b) contribution of the V* terms: 


Od aw 
p* — }. (3.15) 
OX. OXa OXy 


The corresponding expression derived by the 
previous literature reduces to the first two terms 
of (3.13) since all the terms with the derivatives 
of the \‘ drop out for the special choice \‘= const. 

The relation (3.11) permits us to decide 
whether a given expression which contains the 
matter tensor as a factor, integrated over a 
given volume, can or cannot be transformed into 
a surface integral. This depends on whether or 
not the integrand can be written in the form 
T’5gy» with a 6g, which must have the form 
(3.9), with an arbitrary choice of X*. 

If we test, e.g., the volume integrals (2.4) and 
(2.5) which define the total mass and the center 
of mass of the particle, we find that they cannot 
be reduced to surface integrals. They are genuine 
volume integrals. The same holds for the mo- 
mentum quantity j‘ defined by (2.12), which had 
to be added to the mechanical momentum of the 
particle. But the moving force F‘ behaves 
differently. 

Let us make the choice 


(3.16) 
where 7 is a fixed index. Then (3.9) yields: 
5g — (3.17) 


and the general form of the relation (3.11) shows 
that we are able to write: 


(3.18) 


The quantities ¢* assume in the reference system 
(3.12) the following values, after substituting in 
(3.13) and (3.15) the special value (3.16) for X*: 
ik aig k 
=| (3.19) 


OXq OXg OXe 


— p'p*+p* 


These ¢** are symmetric in 7 and k. They differ 
from the ¢;* used by Einstein,’ even for in- 
finitesimal fields, although here one can see 
directly that in this approximation the moving 
force may be calculated from either set of values. 


* Cf. reference 1, p. 806. 
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But the difference appears in the higher ap- 
proximations and our development shows un- 
equivocally that for the evaluation of the moving 
force F‘ the new ¢* have to be employed, and 
not the older ¢;*. 

In consequence of (3.1) and (3.18) the moving 
force can be written as follows: 


(3.20) 


The application of the Gaussian theorem to the 
indices ¢=1, 2, 3 transforms the volume integral 
into a surface integral, extended over the bound- 
ary of the particle, or any closed surface outside 
of the particle, (but including the particle,) since 
in the outer field 7“ remains zero. 

The term with o=4 remains a genuine volume 
integral, contributing to the moving force the 
amount 

Fi=(d/dx,) {t*}. (3.21) 
For static conditions F‘ vanishes. But even for 
nonstatic conditions F‘ remains a time derivative 
and thus must have an oscillatory character and 
cannot contribute anything to the stationary 
part of the force. The stationary part of the 
moving force is determined by the external 
field alone. 


4. THe Motion LAw or A STATIC PARTICLE 
IN FIRST AND SECOND APPROXIMATION 


We now assume a Static particle which shall 
be originally at rest. The acceleration of that 
particle shall be determined, caused by the 
adiabatic superposition of an infinitesimal ex- 
ternal field. For our present discussion the 
original form (3.1) of the moving force shall be 
utilized, without transforming it to a surface 
integral. 

We notice that the force is by its very nature 
a second-order quantity. This is only natural 
if we realize that in first approximation the 
field equations are linear and the superposition 
principle holds. No interference of given field 
and superposed field exists in first approxima- 
tion, and thus the moving force must vanish. 

We now come to the second-order approxima- 
tion. The quantity 


(4.1) 


is a product of two infinitesimal factors. Hence 
both factors have to be known with first-order 
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accuracy only. We can put: 
‘= (4.2) 


abandoning the distinction between covariant 
and contravariant components. 

Now we know that originally the particle is 
at rest, and thus F* vanishes. We only have to 
take the variation of F’, caused by the introduc- 
tion of the infinitesimal external field. We thus 
get: 

However, since the matter tensor 7, is not 

disturbed by the external field in first approxi- 

mation, owing to the non-interference of the two 
fields, we have: 

=0. (4.4) 

Moreover, the variation of the I’,,,; reduces 


itself to the I,,,; quantity of the external field 


which shall be denoted by T,,. ; 

This quantity is practically a constant along 
the particle and can be put before the integral 


sign. We thus get: 
Fis (4.6) 
The previous paper® has shown that all the 
expressions {7} vanish for a static particle, 
with the possible exception of {744}. We know 
further that in sufficiently large distance from 
the center—quite independently of any assump- 
tion of spherical symmetry—g4, may be written 
as follows: 
gu=1—m/r, (4.7) 


where m, the gravitating mass of the particle, is 


determined by: 
m= sa} (4.8) 


Hence | Ts4} =40m. (4.9) 


But, according to (2.4), the same _ integral 
defines also the mass of inertia, neglecting 
quantities of higher order: 


M=4nrm. (4.10) 


The motion law thus takes the form: 
(4.11) 
The mass cancels out on both sides and we get: 
(e) 


at'/dxe= — i- (4.12) 


This equation is equivalent to the “law of the 
geodesic line’’ which has always been considered 
as the natural dynamical law of general rela- 
tivity. We have established this law without any 
reference to a special distribution of matter, 
except the static condition. In particular the 
proof has not involved at all the condition of 
spherical symmetry, nor is it necessary that the 
gia Shall be zero. Hence even rotating particles 
with a stationary spin are included in the law 
(4.12). 


5. THE EXCEPTIONAL CASE 7'=0 


It has been mentioned before that the volume 
integrals which appear on the two sides of the 
dynamical equation, are of two different types. 
The integral on the right side is transformable 
into a surface integral and is thus independent 
of the inner field. Whereas the integrals on the 
left side which define the total mass and the 
center of mass of the particle are genuine volume 
integrals which cannot be transformed into 
surface integrals. In view of this situation it 
seems by no means self-evident that the mass 
shall drop out completely from the law of 
motion. The moving force represents merely 
an interaction of the external field originally 
present with the superimposed external field. 
The mass of inertia, on the other hand, is 
determined by the inner field. It is obviously 
possible to make static changes of the inner gx 
distribution, without modifying the external 
field. If these changes are sufficiently strong, 
they will modify the value of the mass M, 
although the gravitating mass m remains un- 
changed. This would necessarily modify the 
motion law of the particle. 

Such a situation actually arises, even without 
taking recourse to strong fields. In a previous 
investigation® the author has pointed out that 
there is considerable interest connected with a 
matter distribution which is restricted by the 
condition that the scalar 7, the so-called ‘“‘scalar 
of Laue,” shall vanish. A particle of this nature, 
without further specifications but the static 
condition, shows the following interesting fea- 
tures: (a) The total mass m of the particle 
vanishes in first approximation. This is a plaus- 
ible result, since in view of the enormously small 
value of the gravitational forces in comparison 
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with the electrical forces it seems reasonable 
indeed to consider the mass as a second-order 
phenomenon. (b) In second approximation the 
mass comes Out as a necessarily positive quan- 
tity, which again is in perfect agreement with 
the a priori by no means self-evident fact that all 
known masses are positive. 

Now for such a particle gravitating mass 
and inertial mass are no longer reducible to 
each other, without specific conditions. For the 
spherically symmetric case, however, the pro- 
portionality of gravitational and inertial mass 
can still be established. But the numerical factor 
which is ordinarily 47, has to be modified to: 


M=222m/3, (5.1) 


i.e., the customary factor is almost doubled. 

Let us see what happens with such a particle 
if we place it in an external field. We have seen 
that the moving force F‘ depends only on the 
external field of the particle. That external field 
vanishes in first approximation while the field 
of second order has the same form as ordinarily 
a spherical field of first order, except that m 
has dropped down to a second-order quantity. 
Thus the expression 

i, (5.2) 
which was on the right side of the motion law 
(4.11), will remain unchanged. On the left side, 
however, according to (5.1), a new numerical 
factor appears, with the consequence that now 
the acceleration becomes: 


d?ti/dx.2= ,/11. (5.3) 


The law of the geodesic line is no longer fulfilled. 

Let us draw the consequences of the motion 
law (5.3) if adopted for the planetary system. 
The gravitational mass of the sun is not de- 
termined directly by measuring the gq, as it is, 
but indirectly, by the acceleration produced in 
the planets. Hence the motion law of the planets 
has to be known. If the motion occurs actually 
according to the law (5.3), while we assume the 
law of the geodesic line, we underestimate the 
gravitational mass of the sun by the factor 
11/6; i.e., the actual mass of the sun is in the 
ratio 11/6=1.83 larger than the mass ordinarily 
attributed to it. But then all those gravitational 
effects which are not connected with the motion 
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of the planets, should be in fact 83 percent 
larger than customarily expected. 

Now the ingenious and exceedingly careful 
measurements of E. Freundlich and his col- 
laborators® have indeed shown that the light 
deflection on the limb of the sun is larger than 
the theoretically predicted value. The increase 
seems to be between 25 and 40 percent, while 
an increase of about 80 percent seems to be out- 
side of the empirical scope. We must not forget, 
however, that the factor of the relation (5.1) 
has been derived under certain simplifying 
conditions. The assumption 7'=0 has been made, 
which seems to the author, for different reasons, 
a plausible hypothesis. (For one thing, it gives 
a natural explanation why the total mass 
should be a quantity of second order.) But be- 
sides that condition, a static and spherically 
symmetric structure of the metrical field has 
been assumed. Both of the latter assumptions 
are probably generally not satisfied. 

Even so the author believes that the posses- 
sion of the exact law of motion, which dispenses 
with the cumbersome and unreliable successive 
approximations, opens the way to new possi- 
bilities, and may eventually give the clue to the 
proper interpretation of the results of Freundlich 
which need no longer be construed as a con- 
tradiction to the Einsteinian field equations. 

The present paper does not draw all the pos- 
sible consequences of the now firmly established 
dynamical law of a particle. Only static condi- 
tions have been discussed so far. Entirely new 
perspectives open up if we permit inner fields 
to be present of an oscillatory character, and 
study their interaction with a superimposed 
constant external field. Considerations based 
on this possibility led the author to the con- 
viction that the scope of the natural field equa- 
tions 7, =0 for the empty space is actually much 
wider than customarily assumed. One is able to ac- 
count even for the electric interaction of particles, 
without taking refuge to a specifically constructed 
“stress-energy tensor of the electromagnetic 
field,” without abandoning the Riemannian basis 
of geometry, and without introducing a strange 
“fifth dimension.’”’ The author hopes to report 
about these developments at some later occasion. 


4% E. Freundlich, H. v. Klueber, A. v. Brunn, Zeits. f. 
Astrophys. 3, 171 (1931). 


| 
the 
ered 
rela- 

any 

tter, 

the 
n of 

the 
icles 

law | 
ume 

the 

pes. 
able 
lent 

the 

the 
ime 
into 
1 it 
lass 

of 
rely | 
ally 7 
eld. 
is 
isly 

nal 
M, 
un- 
the | 
put 
ous 
nat 
la 
the | 
lar | 
re, 
tic 
q 

on | 


MAY 15, 1941 


PHYSICAL REVIEW 


VOLUME 59 


Electric Breakdown of Glasses and Crystals as a Function of Temperature* 


A. von Hipret Ann R. J. MAURER 
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(Received March 25, 1941) 


The effect of order and disorder upon the breakdown 
characteristics of ionic materials has been investigated. 
While the dielectric strength of the alkali-halide crystals 
increases steeply with temperature in the lower tempera- 
ture range, the trend in glasses is found to be the opposite. 
A comparison of the characteristics of silica glass with 
quartz confirms the fact that the positive temperature 


coefficient of the dielectric strength is produced by the 
periodic structure of the crystals. Also the absolute value 
of the breakdown strength is appreciably increased in the 
disordered amorphous state. The situation is explained 
partly by the influence of temperature and structure upon 
the scattering of electrons. The increase in resistivity of 
metals with temperature is an analogous phenomenon. 


OME time ago it was reported! that the 

dielectric strength of ionic single crystals 
ike KBr shows an unexpected temperature 
dependence (Fig. 1): The breakdown strength 
rises sharply in the lower temperature range, 
reaches a maximum just above room: tempera- 
ture, and falls steeply again. A reasonable ex- 
planation for the increasing strength of the 
material seemed to be offered by the fact that 
with rising temperature the periodicity of the 
structure becomes more and more disordered by 
lattice vibrations. In consequence, electrons 
passing through the lattice will be scattered with 
increasing probability and slowed down; a 
higher field strength must be applied to produce 
breakdown by impact ionization. The increase in 
resistivity of metals with temperature is an 
example of an analogous situation. 

If the periodic order is a decisive factor for the 
breakdown characteristic of a material, the 
increase in breakdown strength should be absent 
in glasses, because the structure is disordered 
from the start. Also their strength should be 
higher than that of comparable crystalline 
materials. Our paper investigates this situation. 


APPARATUS AND PROCEDURE 


Simultaneously and independently of our 
research, Austen and Hackett? and Austen and 
Whitehead* have studied the temperature de- 
pendence of the breakdown strength on several 
materials, including KBr. The specific goal of 

* Publication No. 6 of the Solar Energy Conversion 
Research Project. 

1R.C. Buehland A. von Hippel, Phys. Rev. 56, 941 (1939). 

2 A.E,W. Austen and W. Hackett, Nature 143, 637 (1939). 


3A. E. W. Austen and S. Whitehead, Proc. Roy. Soc. 
A176, 33 (1940). 


this investigation in the British Electrical and 
Allied Industrial Research Laboratory was a 
test of Froehlich’s breakdown theory. Though 
the results obtained on the same material in 
England and in our laboratory (see Fig. 1) 
show a similar trend, they are not in too good 
agreement. It seems, therefore, necessary to give 
an accurate and critical account of our method 
for comparison with that used by the English 
authors. 

Figure 2 shows the latest edition of our break- 
down vessel. As in earlier constructions,! flat 
test specimens are used and clamped by spring 
action between the two highly polished copper 
electrodes E,; and E, of the tubular specimen 
holder shown at the right. The space between 
the slightly rounded electrode surfaces and the 
edge of the sample is filled under vacuum 
with an insulating material possessing a very 


TLL 
4 


TEMPERATURE (°C) 


Fic. 1. The electrical breakdown strength of potassium 
bromide. -O- Buehl and von Hippel, -¢- earlier data, — 
Austen and Whitehead. 


*H. Froelich, Proc. Roy. Soc. A160, 230 (1937); 172, 
94 (1939); Phys. Rev. 56, 349 (1939). 
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ELECTRIC BREAKDOWN 821 


high viscosity at the test temperature. There- 
upon the holder is inserted in the pressure bomb 
at the left and exposed to nitrogen pressure of 
between 50 and 100 atmospheres, while the 
embedding material is plastic. Consequently it 
forms a tight-fitting insulating ring, filling the 
edges between electrodes and sample and pre- 
venting edge breakdown. In addition, the gas 
strengthens the insulation and suppresses corona 
discharges on exposed high voltage leads. 

The ground electrode E,2 contains an insu- 
lated center section of 3 mm diameter for current 
measurements, inserted with one to two thou- 
sandths of an inch clearance into a copper block 
serving for shielding and heat exchange. A 
spring contact connects this current electrode 
through a spark plug L with a feedback current 
amplifier similar to one previously described.° 


S | 


| 
SPECIMEN HOLDER DETAIL 


Fic. 2. High pressure breakdown vessel. 


5S. Robert, Rev. Sci. Inst. 10, 181 (1939). 


(a) (b) 


Fic. 3. Electrode arrangements. EF electrodes, C crystal, 
W wax, G graphite, T tinfoil. Curvature of electrodes and 
thickness of samples exaggerated. 


An Esterline-Angus recorder in the output covers 
full scale the ranges from 1X10-" to 1107 
ampere; its clock-drive is coupled to a voltage 
recorder.! Our judgment about the reliability of a 
breakdown point is based upon the current- 
voltage-time records obtained, and upon a 
microscopic examination of sample and elec- 
trodes after breakdown. Secondary destruction 
of the test piece is prevented by use of a pro- 
tecting resistor (usually 100 megohm) and of an 
automatic switch short-circuiting the electrodes 
immediately after breakdown. 

Electrode E; is inserted in a Pyrex block P 
fitting into the copper wall of the specimen 
holder. By spring contact, connection is es- 
tablished with the high voltage bushing B, a 
Pyrex insulator inserted in a stainless steel cup 
and bolted down upon a cylinder S of the same 
metal. The specimen holder fits tightly into a 
copper block C screwed into the bottom of this 
cylinder and carrying a copper extension rod R. 
This rod, shaped for good heat transfer, pro- 
vides the main temperature regulation; it is 
heated electrically or immersed in a cooling 
liquid. An additional temperature bath sur- 
rounds the copper block and steel cylinder; 
heat insulation is provided by suspending the 
breakdown vessel on three Pyrex rods. The 
temperature of the sample is determined by two 
thermocouples 7 in the wall of the vessel which 
have been calibrated initially against a thermo- 
couple replacing the test sample. 

Our method of avoiding edge breakdown 
(Fig. 3a) was applied successfully between 
—80°C and +150°C by using capacitor oil, 
pyranol, beeswax, and Apiezon wax in over- 
lapping ranges, but for liquid-air temperature 
no solidifying material of good properties has 
been found. The highest strength resulted if 
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Fic. 4. The breakdown strength of soda lime glass. 


nitrogen was slowly liquefied in the breakdown 
vessel from its tank, and used as embedding 
medium under high pressure. This method was 
consequently adopted in the lowest temperature 
range. In rechecking our earlier data in the 
new apparatus, the liquid-air point was found 
lower than before (see Fig. 1), but was certainly 
appreciably higher than the value reported by 
the English authors. 

Austen and coworkers tried to avoid edge 
effects by grinding one side of their test specimen 
into lens shape, painting the curvature with 
graphite, and providing contact over tinfoil 
rings to cylindrical brass electrodes (Fig. 3b). 
If the crystal surface was polished, a low and 
inconsistent breakdown strength resulted, while 
finely ground surfaces gave higher, reproducible 
values. This observation apparently indicates 
contact trouble between graphite and _ the 
polished crystal, which is avoided on the rough 
finish. But grinding of single crystals, which 
split easily, is likely to produce small cracks 
which remain unnoticed under the matt surface 
and create voids and field distortion. Current 
measurements, which would detect such dis- 
turbances, have not been made by the English 
authors. It may be mentioned here that the 
voltage should be raised very slowly to avoid 
possible transient effects produced by polariza- 
tion or formation effects; in our tests about one 
hour elapses between application of the voltage 
and breakdown. 


BREAKDOWN OF GLASSES + 


Figure 4 presents a typical breakdown charac- 
teristic of soda lime glass. The material came 
from one and the same tube and was blown into 
bubble windows then silver-plated in vacuum 


Fic. 5. Thermal and disruptive breakdown of glass. 
Current and voltage vs. time. The dashed line indicates a 
change of the current scale. 


by evaporation, and tested in a thermostat. 
For comparison, samples were blown into thin 
sheets and measured in the breakdown vessel 
(Fig. 2) at room temperature. In the bubble 
windows edge effect is automatically avoided; 
the good agreement between the results obtained 
by both methods serves as a test of our apparatus. 
A comparison between Fig. 1 and Fig. 4 
shows that the steep positive temperature 
coefficient of the single crystal is in glasses 
replaced by a flat negative slope. Contrary to 
earlier observers, we found no range in which 
the breakdown strength was independent of 
temperature, and our current-voltage-time char- 
acteristics seem to indicate that there exist only 
two distinct mechanisms of dielectric break- 
down: At low temperatures disruptive break- 
down is characterized by an abrupt jump of the 
current to high magnitudes without warning, 
while the steep decrease of the characteristic, 
Fig. 4, corresponds to ‘‘thermal’”’ breakdown 
announced by a creeping increase of the current 
at constant voltage. Figure 5 gives examples of 
both situations: a disruptive breakdown of 
Pyrex at 24°C and of soda lime at —39°C, and 
thermal breakdown in soda lime at 40°C and 
84°C. The composition of the sample, its thick- 
ness, the heat transfer, and statistics decide 
when the first mechanism is replaced by the 
second. While in soda lime glass of good ionic 
conductivity under homogeneous stress, a heat 
breakdown can be clearly observed at — 30°C, 
a disruptive breakdown’ may take place at 
+400°C if a pointed cathode is used, enhancing 
electron emission and impact ionization. 


®L. Inge and A. Walther, Arch. f. Elektro. 19, 257 
(1928) ; 22, 410 (1929). P. Moon and A. S. Norcross, Trans. 
A.1.E.E. 49, 755 (1930). 

7 A, von Hippel, J. App. Phys. 8, 815 (1937). 
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Sitica GLAss Versus QUARTZ 


The opposing trends observed in the low 
temperature range of Fig. 1 and Fig. 4 lend 
support to our expectation that the periodic 
order is a decisive factor for the scattering 
probability of electrons in solids. For establishing 
this effect beyond a doubt, the same material, 
SiO2,; was investigated in the crystalline and the 
glass modification. Samples cut in different 
orientations were taken from the same crystal. 
The experiments were made difficult by the high 
breakdown strength of the quartz which necessi- 
tated the grinding and polishing of crystal 
sections down to about 0.005 cm. But the final 
results given in Fig. 6 confirm: The positive 
temperature coefficient is characteristic of the 
ordered state, it disappears and changes to a 
negative slope if the periodicity is destroyed. 

Figure 6 gives more details. At — 80°C, silica 
glass is about twice as strong as the quartz, 
while its density is about 20 percent less. The 
increase in breakdown strength, thanks to the 
amorphous state, is therefore very striking; on 
the other hand, it has to be expected that at some 
higher temperature the crystalline material will 
be stronger than the glass, as observed. That the 
crystal parallel to its axis has a slightly higher 
dielectric strength than when perpendicular to it, 
is somewhat surprising in view of the fact that 
the conductivity in the axis direction is a thou- 
sand times larger than in the other orientation. 
But this conducitivity is ionic and has therefore 
no direct bearing on the dielectric strength, 
while it may indirectly be responsible for the 
effect of producing field distortion by polariza- 
tion. We have no evidence so far, in the case of 
the alkali halides,’ that slow electrons are 
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Fic. 6. The electrical breakdown strength of silica. 
Quartz: -@- field parallel to axis, -o- field perpendicular 


to axis, silica glass: -,-. 


§ A. von Hippel and J. W. Davisson, Phys. Rev. 57, 156 
(1940). 
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Fic. 7. The electrical breakdown strength of silica glass. 
-@- 6X 107 cm, -o- 18 10~ cm, -A- 14 10™ cm. 


sensitive to the lattice direction; therefore we 
might expect no influence of the orientation on 
the strength in quartz crystals. On the other 
hand, we found that the optical axis in quartz 
is not a_ preferential breakdown direction; 
the observed path makes an angle of 41° with 
the axis. Therefore an existing effect may be small. 

It is known that the properties of glasses 
change in wide limits with their prehistory. The 
breakdown strength of quartz glass is no excep- 
tion to this rule. The characteristics of three 
samples of different thicknesses and blown at 
different times (Fig. 7) show the same trend, 
but the absolute value of the breakdown strength 
varies appreciably. 


CONCLUSIONS 


Our experiments have established that the 
periodic order in ionic materials has a decisive 
influence upon magnitude and temperature 
coefficient of their breakdown strength. A 
theoretical calculation of the temperature effect 
has been made by Froehlich‘ under the assump- 
tion that fast electrons are scattered by lattice 
vibrations. The dependence predicted is much 
smaller than the effect observed; in the case of 
silica it should be negligible, while in fact the 
change is of the same magnitude as in the alkali 
halides. We believe that Froehlich’s theory is 
untenable for other reasons, and that the scatter- 
ing effect should be calculated for slow electrons. 
But even then we would not expect too good an 
agreement, because two more physical effects 
enter into the picture, as is reported in the 
following paper. 
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Scattering, Trapping, and Release of Electrons in NaCl and in Mixed 
Crystals of NaCl and AgCl 


A. von HipPet Anp G. M. LEE 
Massachusetts Institute of Technology, Cambridge, Massachusetts 


(Received March 25, 1941) 


The dielectric strength of the alkali-halide crystals decreases steeply with temperature in 
the higher temperature range. This fact cannot be explained by thermal breakdown, but it 
points towards a new mechanism which has been investigated in rocksalt and mixed crystals 
of NaCl and AgCl. The facts observed can be explained by the assumption that the intensity 
of the electron avalanches is decisively influenced by trapping of secondary electrons in the 
lattice, and their release by heat vibrations. New spectroscopic evidence for the trapping 


mechanism is presented. 


N the preceding paper! an intimate connection 
was found between the periodic structure of 
an ionic crystal and the increase of its dielectric 
strength with temperature. Still unexplained 
remains the fact that the breakdown charac- 
teristic, as our measurements on KBr show,’ 
traverses a maximum and falls steeply again at 
higher temperatures. Thermal breakdown cannot 
be held responsible: our current-voltage charac- 
teristics indicate a discontinuous break without 
creeping ; the breakdown path follows preferential 
directions; and moreover, the heat conductivity 
of the alkali halides is comparable to that of 
iron. To learn more about the situation, we 
investigated the breakdown characteristic of 
rocksalt and the influence of increasing amounts 
of AgCl upon it. 


MIXED CrYsTALS OF NACL AND AGCL 


Sodium ions in rocksalt can be replaced to 
some extent by silver ions, effecting a remarkable 
increase in the breakdown strength of the crystal 
at room temperature. This had been previously 
established® but because the research had to be 
interrupted, a characteristic of dielectric strength 
versus composition had been given, based upon 
the concentration of the melt.* In our new 
investigation we have checked the silver con- 
centration of the crystal samples used: pieces of 


1A. von Hippel and R. J. Maurer, Phys. Rev. 59, 820 
1 


(1941). 
8 35) C. Buehl and A. von Hippel, Phys. Rev. 56, 941 
1939). 

3A. von Hippel, Zeits. f. Physik 88, 358 (1934). 

* The crystals used by us were grown by the Kyropoulos 
method (S. Kyropoulos, Zeits. f. anorg. Chemie 154, 308 
(1926)). The samples used in breakdown measurements 
were split from them to between 1 to 2 107 cm thickness. 


the same section were dissolved in a 0.01 normal 
solution of KCN containing an addition of KI 
and ammonia. The solution was back-titrated 
with a 0.01 normal AgNO; solution, and the 
turbidity point observed.‘ This method gave 
very consistent and reproducible results in the 
crystals weighing as little as 76 g and showed 
that 10 percent, 5 percent, and 1 percent molar 
concentration of AgCl in the melt corresponded 
to about 3 to 5 percent, 1.5 to 3 percent, and 0.3 
to 0.5 percent concentration in the crystal. The 
silver content of a mixed crystal showed local 
variations from section to section, according to 
the concentration of the remaining molten ma- 
terial. 

Figure 1 gives our new characteristic of break- 
down strength versus composition at 22°C. The 
values measured agree well with the old ones if 
the changed horizontal scale is taken into ac- 
count. Preliminary measurements indicate that 
the breakdown strength of the clean silver 
chloride is of the order of 1.7X10° v/cm. 


3.0 T 
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Fic. 1. Breakdown strength of mixed AgCl-NaCl 
crystals at 22°C. 


4 J. M. Kolthoff and J. B. Sandell, Textbook of Quantita- 
tive Inorganic Analysis (Macmillan, New York, 1936), 
p. 545. 
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SCATTERING OF ELECTRONS IN NaCl 825 


Mixed crystals with more than about 7 percent 
of silver content are unstable, and the samples in 
the limiting range are inclined to form multi- 
crystals and to crack easily. Consequently the 
characteristic can be measured only over a 
small range of concentrations, but its general 
shape and the known curves for mixed crystals 
between alkali halides* seem to imply that the 
increase in dielectric strength is an effect analo- 
gous to the increase in resistivity observed in 
metallic solutions. 


BREAKDOWN CHARACTERISTICS AS A FUNCTION 
OF TEMPERATURE 


Using the apparatus and procedure described 
in the preceding paper, we have measured the 
temperature characteristics of our mixed crystals 
(Fig. 2). The result is surprising and reveals a 
more complicated picture than was anticipated. 
At low temperatures the breakdown character- 
istics of pure NaCl and of the mixed crystals 
have about the same slope; increasing silver 
content displaces the curves to higher values. 
The higher the silver concentration the sooner 
the dielectric strength reaches a maximum and 
the steeper is the negative slope thereafter. 
Consequently the curves cross, and while at 
room temperature and below, the breakdown 
strength of rocksalt is increased by the addition 
of silver, the opposite is true above 100°C. 

It has been demonstrated in the preceding 
paper! that the motion of electrons in ionic 
crystals is impeded by temperature vibrations 
disordering the lattice structure. Nevertheless, 
the formation of electron avalanches producing 
breakdown is favored by heating of the material, 
as the falling part of our characteristics indicates. 
It can only be concluded, if the avalanche pic- 
ture is accepted, that the number of free electrons 
produced by impact per unit length of path 
increases rapidly with temperature. 

At a first glance this conclusion seems strange 
but it can be supported by other evidence: We 
are accustomed by Townsend’s theory of gas 
breakdown to assume that electrons multiply in 
avalanches in exponential fashion: 


n=MNo f (1) 
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Fic. 2. Temperature dependence of the breakdown 
strength of mixed AgCl-NaCl crystals. 


This formula assumes that the secondary elec- 
trons stay free and participate in the ionization 
process; the factor a describes the probability of 
ionization. For breakdown in electro-negative 
gases this is not true any more, because electrons 
may be trapped while they are slow. The intensity 
of ionization decreases with increasing proba- 
bility of capture, c: 


n=nNo esp] (2) 


A discharge is quenched by the addition of 
electro-negative gases. The negative ions formed 
may, on the other hand, release the trapped 
electrons again under the influence of violent 
collisions, as Loeb has first shown’ for the O2~ ion. 
Consequently a probability of reliberation, r, 
has to be introduced, and the increase of an 
avalanche, started by mo primary electrons and 
growing by impact ionization, has to be described 
by the general expression : 


exo] fa (3) 


This situation apparently exists in our ionic 
crystals. Electrons are trapped in the alkali 
halides, as the appearance of the ‘F”’ band indi- 
cates ;° the spectral position of this band corre- 
sponds to the energy needed for their photo- 
electric liberation. Electrons can be reliberated 


5 L. B. Loeb, Phys. Rev. 48, 684 (1935). 
a an for instance, R. W. Pohl, Proc. Phys. Soc. 49, 3 
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——— ABSORPTION BANDS 
PRODUCED BY X-RAYS 


30 MIN. OF A+ 435 mp 
CRYSTAL AT -180°C 
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Fic. 3. Absorption bands in NaCl. 


by heat vibrations, as Stasiw first found,’ be- 
cause the band of trapped electrons at high 
temperatures moves in an electric field like an 
electronic cloud. The probability of trapping is 
increased by the addition of AgCl to NaCl;° the 
migration of the color centers slows down and 
can be entirely suppressed if sufficient silver is 
added. If we could show that at moderate tem- 
peratures the probability of reliberation is suff- 
cient for the electrons to be reinstated quickly 
into the avalanches, our results of Fig. 2 would 
be understood. The increase in breakdown 
strength with temperature would be due. to the 
scattering and trapping of electrons; the de- 
crease would result from their release, which 
apparently is favored by mounting silver con- 
centration. 

Avalanches producing breakdown will have a 
critical magnitude. On the basis of the mecha- 
nism proposed above, one would expect that the 
length of these avalanches is dependent on tem- 
perature. For sufficiently thin samples (<10~ 
cm) the breakdown strength should therefore be 
a function of thickness varying with temperature. 


THE TRAPPING OF ELECTRONS 


If the ‘“F”’ band, as one of us believed until 
recently,? corresponds to the trapping in the 
free lattice, the chance of liberating electrons 
from the band indicated by this spectral position 
through thermal motion at the temperatures of 
our experiments would be small indeed. On the 


70. Stasiw, Goettingen Nach. No. 50 (1933). 
8 A. von Hippel, J. App. Phys. 8, 829 (1937). 
* A. von Hippel, J. Chem. Phys. 8, 605 (1940). 


other hand, de Boer’ and Mott and his co- 
workers" have proposed another model for the 
“F” band: an electron captured at a lattice 
defect where a negative ion is missing. If that is 
true, a band corresponding to the electron in the 
lattice should be found at a longer wave-length, 
because this state would be an excited state of 
the first one. To decide this question, J. P. 
Molnar has at our suggestion reinvestigated the 
absorption spectra of the alkali-halide crystals 
between 10,000 and 2000A. If the band in ques- 
tion exists, it should be found by exposing the 
crystals to x-rays at low temperatures. Electrons 
ejected into the lattice under this condition could 
be trapped with reasonable probability in the 
flatter but much more frequent potential hole of 
the free lattice? and stay there, while at higher 
temperatures they would be shaken out fre- 
quently enough so that they would all be found 
in the deeper hole corresponding to the ‘F” 
absorption. 

This expectation has been fulfilled. Molnar has 
found in all the alkali halides a band accompany- 
ing the ‘‘F”’ band at longer wave-lengths. It is of 
the small width to be expected for atomic centers, 
can be bleached, and disappears by heating.* 
Figure 3 shows the situation for NaCl, indicating 
the position of the band at 1.75 ev. A theoretical 
calculation for the position of the band, made by 
a cycle process before this result was known,’ 
gives 1.63 ev; the agreement, while encouraging, 
is better than the theoretical data justify. 

That electrons are trapped high above the “F” 
band in our mixed crystals, and then released by 
heat vibrations, can be demonstrated in a differ- 
ent way. Under x-ray illumination the crystals 
are strongly phosphorescent at room tempera- 
ture. The light emitted is bluish, and preliminary 
measurements indicate that the end state of the 
electrons lies in the near ultraviolet shortly below 
the ‘‘F’’ band. Measurements of the initial state 
and its sensitivity to temperature and _ field 
strength and of the light emitted will give us, we 
hope, the final link between the temperature 
dependence of the breakdown strength and the 
spectroscopic properties of our crystals. We are 
engaged in this study. 


( 107, H. de Boer, Rec. Trav. Chim. Pays-Bas 56, 301 
1937). 

tN. F. Mott, Trans. Faraday Soc. 34, 822 (1938). 

* Doctorate thesis of J. P. Molnar, Massachusetts 
Institute of Technology (1940), to be published shortly. 
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When a disk cut from a cubic crystal of a ferromagnetic material is placed in a magnetic field 
parallel to its plane, the torque exerted on the disk by the field generally increases continually 
with the strength of the field, approaching a finite limit (saturation) in very high fields. An 
“anomaly” for one orientation of a crystal of iron-silicon alloy has recently been observed by 
Tarasov. We have confirmed his findings and observed a similar effect for a series of orienta- 
tions of disks cut parallel to the (110) and (100) planes. As the field strength is increased, the 
torque passes through a well-defined maximum before beginning the final approach to satura- 
tion. This anomalous peak is now explained in terms of the domain theory; no new concepts 
are necessary. A graphical method is described for calculating torque and magnetization curves 
for any direction of the applied field. Torque curves depend markedly on the demagnetizing 
factor of the disk and, in many cases, a high demagnetizing factor prevents observation of a 
maximum in the curve of torque vs. field strength. 


HEN a disk cut from a single crystal of a 
ferromagnetic material is placed in a 
magnetic field, it is generally subject to a torque 
on account of the anisotropic character of the 
crystal, even when the crystal is cubic. The 
anisotropy constant of the crystal is probably 
best determined by measurement of this torque. 
In most experiments hitherto performed the 
torque increases continually with field strength 
and approaches a finite value simply related to a 
single anisotropy constant, and the observed 
relation between torque and field strength agrees 
with calculation. 

Recently, however, Tarasov! has found that 
under certain conditions the torque passes 
through a maximum before declining to its 
limiting value in high fields. The purpose of this 
note is to show that this supposedly anomalous 
behavior is in accord with theory, and to give a 
method of calculating torque curves for any 
orientation of the crystal with respect to the 
field. Theory predicts, and our experiments 
show, that a maximum occurs in the torque vs. 
field strength curve for most crystal orientations 
provided the demagnetizing factor of the speci- 
men is small. In order to have a small demag- 
netizing factor our experiments were performed 
with a single crystal disk 22 mm in diameter and 
0.22 mm thick. 


1L. P. Tarasov, Phys. Rev. 56, 1224 (1939). 


SIMPLE TORQUE CURVES 


When the specimen of magnetic material, 
disk or oblate ellipsoid, is cut with its principal 
plane parallel to the (100) plane of a single 
crystal, the equations for calculating the torque 
per unit volume L, and the magnetization J, in 
the presence of a field of strength //, are derived 
from the expression for the energy density E. 
The energy density is composed of two parts, 
Ex due to the magnetic anisotropy of the crystal 
and Ey, due directly to the presence of the field. 


E=Ext+En 
= — HI, cosé. (1) 


Here K;, is the anisotropy constant of the cubic 
crystal, the s’s are the direction cosines of the 
local magnetization, with respect to the three 
crystal axes and @ is the angle between the direc- 
tion of the field and that of the local magnetiza- 
tion. In accordance with the domain theory it is 
assumed that each domain is magnetized to 
saturation (J=J/,); the torque acting on any 
domain is then 


L=—dEx/dé (2) 


ergs/cm* and the torque on the whole specimen 
is the resultant of the torques on all of the 
domains. The strength of the field is obtained by 
minimizing E and is 
dEx/dé 
I, sind 


(3) 
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When the specimen is demagnetized the result- 
ant is zero. When the field is increased in strength 
the directions of magnetization in most of the 
domains are changed by displacement of the 
boundaries between domains, until they are 
directed along that direction of easy magnetiza- 
tion that is nearest to the direction of the field; 
when this process is near completion the torque 
first attains appreciable magnitude, and the 
calculations of this article begin at this point. 


y 
i \ Is 
A; \ / 
\ 


Fic. 1. Single crystal disk in a magnetic field showing 
directions of easy magnetization in the (001) plane, and 
nonparallelism of H and J;. 


When all of the domains are similarly oriented 
Eqs. (1) to (3) are easily applied. In the (001) 
plane (see Fig. 1) Eqs. (2) and (3) reduce to 


L=K, (sin4a)/2, 
H=K;,sin4a/2I, 


in which @ is the angle between J, and the 
nearest crystal axis [100]. . 
The torque in this plane is greatest when a is 
22.5°. H and L have been calculated using the 
above equations with a+@=22.5°, the field 
being applied always 22.5° from a crystal axis. 
In any experiment the strength of the true field 
H is not equal to that of the applied field, Hz, 
because of the demagnetizing field, AH, which 
is determined by the magnetization and dimen- 
sions of the specimen. For comparison with 
experiment 160 oersteds (the demagnetizing 
field for saturation) were added to H to obtain 
H,, and L was plotted against H, as shown in 
Fig. 2. The experimental curve shown for com- 
parison was obtained using a very thin specimen 
cut and ground parallel to a (001) plane of a 
crystal containing 3.6 percent silicon and the 
remainder iron, prepared as described by Wil- 
liams.2 The ratio, m, of diameter (22 mm) to 
thickness (0.22 mm) was 100, and the corre- 
sponding demagnetizing factor is N=0.10. 
2H. J. Williams, Phys. Rev. 52, 747-751 (1937). 


The demagnetizing field when the magnetization 
is near saturation is thus NJ,=160. Both experi- 
mental and theoretical curves are consistent with 
an anisotropy constant K,=350,000 ergs cm=-, 
in fair agreement with results of Tarasov.* 

The agreement between theory and experi- 
ment indicates that for this simplest case the 
theory is good. 


GENERAL CASE 


Lines of constant crystal energy, Ex, are 
plotted in stereographic projection in Fig. 3, 
using Eq. (1). The angles ¢ and y are latitude 
and longitude and the spherical triangle shown 
includes all points of the complete sphere that 
are nonequivalent crystallographically in a 
cubic crystal. Each point on the projection 
represents a direction determined by the corre- 
sponding point on the sphere, and the center of 
the sphere. 

This diagram can be used to find the torque 
on a crystal when placed in a field having the 
direction indicated by any point H on the 
diagram, as shown in Fig. 4. When the field is 
weak the point representing the magnetization 
of the domain J, will be at [100], the nearest 
direction of easy magnetization, and when the 
field strength is increased the rotation of J, will 
be represented by some line going from [100] to 
H. The path followed by J, will be straight if H 
lies on the line between [100] and [110], or on 
the line between [100] and [111], and the 
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Fic. 2. Torque/volume vs. strength of applied field. 
Solid line is observed relation, dashed line the relation 
calculated from simple domain theory using the previously 
known value of the anisotropy constant K,=350,000. 
The intercept on the H, axis is the demagnetizing field de- 
termined by the dimensional ratio of the disk (NJ, = 160). 


3L. P. Tarasov, Phys. Rev. 56, 1231-1240 (1939). 
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Fic. 3. Lines along which the density of the energy of magnetic anisotropy E is constant. [See Eq. (1).] 
Stereographic projection of a cubic crystal with scales of latitude ¢ and longitude y. 
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Fic. 4. The position of the magnetization vector J, of a 
domain is determined graphically as shown here for a given 
orientation of the field H and angle @ between H and J,. 


calculation of the torque can be made most 
easily in these cases by using the equations that 
have been used repeatedly for this purpose.‘ 
If H is applied in any other direction the path 
followed by J, will be curved and may be de- 
termined graphically in a way to be described. 
Whatever the path that the end of the /, 
vector follows, it is evident that a maximum in 
L may be expected when J, lies in the region in 
which the lines of constant crystal energy are 
closest together, because L is the derivative of 
Ex with respect to 6. In Fig. 3 the region in 
which the equipotential lines are closest together 
is marked with small arrows in the (001) and 
(011) planes. A maximum in the L vs. H curve 


*E.g. by R. M. Bozorth, Phys. Rev. 50, 1076-1081 
(1936). 
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Fic. 5. Curves showing torque vs. field strength for disk 
parallel to (001) plane. Theoretical curve is for infinitely 
thin disk; observed curves are for disks having ratio of 
diameter to thickness of 100 and 7.4. 


is thus to be expected when H and [100] lie 
on opposite sides of this region. For example, in 
the (001) plane the field must “‘pull’’ the J, vector 
over the place (y= 22.5°) where L is a maximum, 
and the same maximum value of Z occurs in the 
L vs. H curves for each position of H between 
y= 22.5° and y=45°. It may thus be determined 
that for over 80 percent of all of the possible 
orientations of HT with respect to the crystal 
axes, L passes through a maximum. As men- 
tioned below, the existence of too large a de- 
magnetizing factor may in some cases cause the 
maximum to disappear. 

To illustrate the determination of the path 
and of the L vs. H curves, two cases will be con- 
sidered for each of which the torque passes 
through a maximum and declines to its final 
value. In the first case H/ is applied in the (001) 
plane 40° from the [100] direction (¢=0°, 
y=40°), in the second it is applied in the (110) 
plane 75°, 70° or 65° from the [001] direction 
(@=15°, 20° or 25°, y=45°). In the first case 
the path is obviously a straight line, in the 
second case the way in which the path is to be 
determined may be outlined as follows. Draw 
circles of various diameters about the point H 
and note the points at which these circles are 
tangent to the equipotential lines. At these 
points the energy Ex+ Ey is a minimum and the 
locus of such points is the required path. Along 
this path Ex is plotted against 6, the angle 
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between J, and H, and the slope of this curve 
is —L as given by Eq. (2). His found by dividing 
L by I, sin@ in accordance with Eq. (3). Paths 
and curves so determined are shown in what 
follows: 


(001) PLANE 


In the (001) plane, with // directed 40° from 
[001], the appropriate equations are Eqs. (2) 
and (3) with 

6=40°—a. 


The curve so calculated with K,=350,000 is 
shown in Fig. 5. 

For comparison with theory data were taken 
on the thin crystal described above. The ob- 
served curve shows the expected maximum and a 
limiting value in high fields that agrees with the 
calculated limit. However, the height of the ob- 
served maximum lies considerably below the 
theoretical one. This may be attributed to the 
action of the demagnetizing field, and this point 
of view is substantiated by the fact that no 
maximum at all was observed when a much 
thicker disk was used having a ratio of diameter 
to thickness of m=7.4. 

The way in which the demagnetizing field 
lowers the value of LZ in intermediate fields is 
shown in Fig. 6. The applied field 77, has turned 
the magnetization vector 7, through the angle a 
from the nearest direction of easy magnetization 
[100]. In this position J, gives rise to an anti- 
parallel demagnetizing field, AZZ, which adds 
vectorially to H, to form the true field, J/. 
Now H will be nearer to the [110] direction than 
will HZ,; consequently, if AH or the angle be- 
tween HT, and [100] is large enough, the true 


Fic. 6. Orientations of vectors 7, and H=H,—AH in 
(001) plane, to illustrate effect of demagnetizing field AH 
on the direction of H and the magnitude of the torque. 


\ (001) 3 
/ \ 
| “ 
| \ \ 4s° | 
\ BH 
f 
/ 


rve 
ing 
ths 
hat 


om 


(2) 


\H 


TORQUE ON A SINGLE CRYSTAL 831 


x103 
200 
PARAMETER, ANGLE 
Ape BETWEEN Hg AND [100] 
OR [010] IN DEGREES 
160 + 
\ 
36.7 
120 
Fic. 7. Curves of torque 
vs. reciprocal of applied § 
field strength for various §& 
directions of the field ap- 
ied in the (001) plane a 
tion to Hi = ~. « 
48.7 
63.7 
40 
) 8 16 , 24 32 40 48 56 
Fig? RECIPROCAL OF FIELD STRENGTH x 105 


field will lie almost parallel to [110] and some of 
the domains will spring to another equally stable 
position lying between [110] and [010]. These 
domains will have a torque opposite in sign 
to those lying on the other side of [110], and 
the resultant of the torques of all the domains 
will be considerably lowered. The thicker the 
crystal or the smaller its diameter the greater 
will be AH and the smaller the torque in all 
except very intense fields. Since the crystal 
used had a ratio of diameter to thickness as high 
as 100, it is not surprising that such a maximium 
has never been observed before. 

Other data for this specimen for other posi- 
tions of the field with respect to the crystal 
axes are shown in Fig. 7. Here L is plotted vs. 
the reciprocal of H, in the way suggested by 
Schlechtweg> and used by him and Tarasov.' 
It may be noted that, in general, the portions of 
the curves for which H, is largest are not as free 
from curvature as those published by Tarasov, 
but this method of plotting is nevertheless a 
convenient one for extrapolating to infinite 
field strength. 


(110) PLANE 


When a disk is cut parallel to a (110) plane, 
and magnetized in a direction lying between 


5H. Schlechtweg, Ann. d. Physik [5] 27, 573 (1936). 


adjacent [100] and [111] directions, the calcu- 
lation of L and H can be made by the method 
often described.‘ But when J// lies between 
adjacent [110] and [111] directions the problem 
is complicated by two factors: The path followed 
by J, is curved, and there are two equidistant 
directions of easy magnetization. Heretofore the 
torque in the (110) plane has been calculated for 
all azimuths by assuming that the magnetization 
vector of each domain is confined to the plane 
by the high demagnetizing field in the direction 
perpendicular to the plane; but this assumption 
is invalid because the magnetization in half of 
the domains may have a considerable component 
perpendicular to the plane in one direction and 
the other half an equal component in the op- 
posite direction so that the over-all magnetiza- 
tion is confined to the plane. 

In the illustrative example to follow it is 
assumed that in the weaker fields the magnetiza- 
tion is equally divided between the two direc- 
tions [100] and [010]. Calculation of the path 
of the vector rotating from just one of these 
positions is, of course, sufficient because the 
path of the other is the mirror image of the first. 

Figure 4 shows the graphical construction of 
the path of the vector, J,, as it rotates from 
[100] to H which is located at ¢=15°, y=45°. 
The contour lines are the same lines of constant 
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crystal energy as those of Fig. 3. Arcs of circles 
are drawn, each point on which represents a 
vector making the same angle 6 with the field, 
H. The point of tangency of this circle with a 
line of constant crystal energy is a point on the 
required path followed by J,, for the crystal 
energy is constant along the one line and the 
energy due to the field, HJ, cos@, is constant 
along the circle and the sum of the two energies is 
thus a minimum at the point of tangency. Arcs 
for other values of @ determine points of minimum 
energy for other positions of J, corresponding 
to other field strengths, and the locus of such 
points is the required path of the vector, J,. 

Lines of constant @ are circles on the stereo- 
graphic projection as well as on the sphere, but 
the centers of the circles on the projection are 
not at the projection of the center of the circles. 
The center of the projection can be found simply 
in the following way. Draw a straight line from 
[100] through H. Determine the angular dis- 
tance, ao, from [100] to H by laying off from 
[100] along the projection of the (001) plane a 
linear distance equal to that from [100] to H 
and reading the angle on the y scale. On the 
straight line passing through [100] and H, 
mark off the angles ao+@ and ao—8, using the 
same scale. These two points are at opposite 
ends of a diameter of the required circle which 
may then be easily drawn. The distance / on the 
projection from [100] to any point y=y, 
¢=0, is given by 


l=ly tan(y,/2), 
where /, is the radius of the projected sphere. 
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The next procedure is to plot the crystal 
energy Ex as a function of 6, the angle between 
I, and H, for all points on the path of J,. In 
accordance with Eqs. (2) and (3) the slope of 
this function determines the torque Z,; and the 
field strength H: 


L,= —dEx/dé, 
dEx/d@ Ly 


I, sine 


for all domains originally parallel to [100]. 
The torque due to all of the domains, those 
originally parallel to [100] and to [010], is 
best determined by finding the magnitude and 
direction of the resultant J of all of the domains 
and multiplying this by the sine of the angle, 6’, 
between this and H. Thus, the torque on each 
cubic centimeter of specimen is 


L=HI, cos@ siné’, 


where 8 is the angle between J, and the (110) 
plane measured along the great circle perpendicu- 
lar to the (110) plane. Approximate values of the 
angles 8 and @’ are easily read from Fig. 7; 
for the case under consideration calculations 
have shown that the angle 8 is very nearly equal 
to the difference between 45° and the value of y 
for the position of J,, and 6’ is nearly equal to the 
difference in the values of @ for J, and for H. 
If necessary a more accurate evaluation of 8 and 
6’ can be made analytically, or graphically using 
a stereographic net. 


| 
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Finally, the component of the resultant J 
that is parallel to H is 


I=TI, cosé. 


This quantity is useful in determining the de- 
magnetizing fields that must be calculated in 
order to compare the theoretical torque curves 
with the results of experiment. 

The results of the calculation are shown in 
Fig. 8. Three different directions were chosen 
for the application of the field, those at 15°, 
20° and 25° from the [110] direction in the 
(110) plane (Y=45°, ¢=15°, 20° and 25°). The 
heavy dashed lines were calculated as described 
above; the light dashed lines show the relation 
expected on the former assumption according to 
which the magnetization in each domain is 
confined to the plane of the disk. For comparing 
the theory with experiment the applied field, 
H., is calculated by adding to the true field, H/, 
the demagnetizing field, AH, appropriate for a 
disk having a ratio of diameter to thickness of 
4.0. The experimental curve is shown in the 
same figure as a heavy solid line, and was ob- 
tained using a disk 0.8 cm in diameter and 0.2 
cm thick which had its edges cut away to make 
the demagnetizing field more nearly that of an 
oblate ellipsoid of the same dimensional ratio. 
The demagnetizing factor used in the calcula- 
tions was N=1.75 and the demagnetizing field 
strength at saturation is therefore 2800. The 
other constants used were J, = 1600, K = 350,000. 

The calculated curves bear a close resemblance 
to those observed. It is to be expected that the 
calculations will not be exact, for no account has 
been taken of the deviation of the direction of the 
true field from that of the applied field, a devia- 


tion similar to that illustrated in Fig. 6 for the 
(100) plane. 

It should be pointed out that the magnitude of 
the demagnetizing field, or the field at which the 
torque curve leaves the // axis, is in itself strong 
confirmation of the truth of the assumption that 
the domains are directed initially along [100] 
and [010], rather than along [110] as formerly 
assumed. In the latter case the calculated torque 
curve begins to rise at 2800 cos20° or 2600 
oersteds, in the former case at 2800 cos46° or 1900 
oersteds, when H, is applied at ¢=20°. The ex- 
perimental curve shows an initial rise at 1900 to 
2000 oersteds. 


x10? 
H= 
5000 


200 
PARAMETER, ANGLE BETWEEN 
Hag AND [100] IN DEGREES 


23.5 
31.9 


L, TORQUE IN ERGS PER CM? 
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a3 


41.9 
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Fig, RECIPROCAL OF FIELO STRENGTH x1o~4 


Fic. 9. Torque vs. reciprocal of applied field strength 
for various directions of the applied field in the (110) or 
(110) plane (see Fig. 4). 


Curves showing LZ as dependent on 1/H are 
given in Fig. 9 for other directions of the applied 
field in the (110) plane. 

It would appear that the simple linear rela- 
tion observed by Tarasov! to exist between L 
and 1/H, at high values of H., has no general 
experimental or theoretical validity, but refers 
only under limited conditions of crystal orienta- 
tion and specimen form. 
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Letters to the Editor 


ROMPT publication of brief reports of important dis- 
coveries in physics may be secured by addressing them 
to this department. Closing dates for this department are, for 
the first issue of the month, the eighteenth of the preceding 
month, for the second issue, the third of the month. Because 
of the late closing dates for the section no proof can be shown 
to authors. The Board of Editors does not hold itself responsible 
for the opinions expressed by the correspondents. 
Communications should not in general exceed 600 words 


in length. 


Effect of Neutron Energy on the Total Decay 
Curves of Fission Products 
S. K. HAYNES 


Brown University, Providence, Rhode Island 
April 7, 1941 


HE study of the total decay curve & the fission 

products of U collected by recoil was first made by 
Joliot' using chiefly Rn-Be neutrons. A marked increase 
in the ratio of the intensity of long periods (~1 hour) to 
that of short periods (~15 min.) was produced when the 
source and target were surrounded by paraffin. Later, 
Bjerge, et al.,? using neutrons from Li plus deuterons 
(0.75 Mev) found that within +10 percent the decay 
curves were identical for U+fast neutrons, U+C neutrons, 
and Th-+fast neutrons. This proved true for both irradia- 
tions of a few minutes and of 100 minutes. Frisch* showed 
that the results of Bjerge et al. (Table I) fitted reasonably 
well the curve 


Jr(t) =const. "—(t+T)-"), (1) 
where T is the time of irradiation, ¢ the time of decay 


and n is from the Sargent law \=kE”" (n=5). This result 
is calculated on the assumption that all energies have 
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equal statistical weight. Thus the results of Bjerge et al. 
would not indicate that the products of U2%5 (C neutrons) 
were necessarily identical with those of U5 (fast neutrons) 
but only that the number of periods present was so large 
as to smear out any differences. 

In the spring of 1940 at Paris, the author, at the sug- 
gestion of Professor Joliot, undertook to resolve, if possible, 
the discrepancies between the Paris and Copenhagen 
measurements and also to extend them to higher energies. 
Very incomplete measurements (cut off by the fall of 
France) with Rn-Be neutrons, seemed to confirm the 
difference previously observed by Joliot. Furthermore the 
results with C neutrons, in accord with Joliot’s previous 
work, failed to agree with the Frisch formula by a factor 
of 2. 

The experiments were continued for high energies 
during the summer of 1940 with neutrons from the Harvard 
cyclotron. The results are shown in curves A-E, Fig. 1, 
and in Table I in descending order of neutron energy.‘ 
The deuteron beam was controlled by hand and the mean 
current was estimated and held constant by taking readings 
every half-minute. The good agreement between the 
successive runs of curve A, and also of curve B, testifies 
to the reproducibility of the decay curves. Except possibly 
for curve C, deviations due to varying beam current were 
undoubtedly less than 5 percent. 

The counting system consisted of an alcohol argon 
counter with 0.010” brass walls. It was used with a scale 
of 32. The whole was checked for additivity and a resolving 
time of 2.3 10~ sec. was found. A blank fission-catcher- 
was placed behind the actual catcher on each run and its 
activity was subtracted. A constant radioactive source 
was used at intervals during each run to check the con- 
stancy of the counter. 

Curve E (thermal and low energy fast neutrons) agrees 
excellently with Frisch’s equation (Table I) and thus with 
Bjerge’s curves for both thermal and fast neutrons. 

The shape of the other curves does not vary rapidly 
with energy, although numerous additional products® are 
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TABLE I. Relative number of counts from decay products. 


MIN 
AFTER | FRISCH | FRISCH F 
IRRADI- 60 1 BJERGE 
ATION MIN MIN A B Cc D E ET AL 
4 28.1 20.4 36.6 |33.7 (31.4 30.4 |28.0 19.7 
10 17.6 13.5 23.5 |22.0 |20.4 {19.3 |18.6 13.2 
20 11.6 9.3 15.3 (14.4 |14.2 |15.5 |12.3 9.81 
40 6.97 5.98 8.06} 7.85| 7.55) 7.78) 6.98 6.05 
60 4.92 4.40 5.28] 4.95) 5.06) 4.99) 4.64 4.29 
100 3.15 2.93 2.94} 2.84) 2.80) 2.86) 2.79 2.66 
140 2.25 2.14 2.11) 2.08) 2.15) 2.06) 1.95 2.02 
180 1.74 1.68 1.58) 1.58) 1.61) 1.51} 1.56 1.62 
240 1.27 1.26 1.25] 1.19) 1.25) 1.20) 1.18 1.26 
300 1.00 1.00 1.00) 1.00) 1.00) 1.00) 1.00 _ 


known to result from fission at high energies. There does 
seem, however, with increasing energy, to be a definite 
tendency towards shorter periods (Table I) which amounts 
to nearly 30 percent in the case of bombardment with Li 
neutrons. Since these are composite curves of many 
parent and daughter products, a 20-30 percent change at 
4 to 10 min. might indicate a considerably larger increase 
at } to 1 min. It thus appears that not only are new 
products formed, but also the new parent substances have 
shorter periods, and their probabilities of formation are 
relatively large at neutron energies of 10-25 Mev. 

Circumstances have prevented the continuance of this 
work such as the investigation of very short periods and 
of the irregularity in curve D and a conclusive check on 
the Joliot-Bjerge discrepancy. Because of larger intensities, 
the Copenhagen results are probably more reliable. 

The author wishes to express his gratitude to Professor 
Joliot and to the Institute of International Education and 
Williams College, under whose auspices he was privileged 
to work in Paris, and to acknowledge with deep apprecia- 
tion the kindness of Harvard University, and the assistance 
of Drs. I. A. Getting, E. P. Clancy and L. Fussell. 

1 F. Joliot, J. de phys. et Rad. 10, 159 (1939). 

2T. Bjerge, K. J. Brostrém and J. Koch, Nature 143, 794 (1939). 

30. R. Frisch, Nature 143, 852 (1939). 

4The neutrons of curve C were considered to be of higher energy 
than those of curve D because the activity of the control relative to 


the fission was 15 percent less for D than for C, 20 percent less than 
for B,and <} of A. 

§ Y. Nishina, T. Yasaki, H. Ezoe, K. Kimura and M. Ikawa, Nature 
146, 24 (1940); Y. Nishina, K. Kimura, T. Yasaki and M. Ikawa, 
Phys. Rev. 58, 660 (1940); 59, 323 (1941); E. Segré and G. T. Seaborg, 
Phys. Rev. 59, 212 (1941). 


Lifetimes of Nuclear Levels with Respect to 
Electric Multipole Radiation 


IRVING S. LOWEN 
Washington Square College, New York University, New York, New York 
May 1, 1941 


HE liquid drop model has been used to estimate the 
lifetimes of nuclear excited levels with respect to 
electric multipole y-ray emission involving various changes 
in angular momentum and different excitation energies.' 
The purpose of this note is to give the results of a calcula- 
tion of these lifetimes by a method which avoids some of 
the approximations usually made, beyond those already 
embodied in the liquid drop model. The distortion of the 
nucleus is assumed to be given by 


TABLE I. Lifetimes of nuclear states for various excitation energies and 
angular momentum changes. 


ENERGY 
or Ex- CHANGE IN ANGULAR MOMENTUM 
CITATION 
(EV) l=3 l=4 
sec. 0.83 hr. 3.5 yr. 2.0 X10 yr. 
X104 sec. 0.20 sec. 33 days 2.0 X 108 yr. 
20 X104 |9.3 sec. 4.8X10™ sec. 43 sec. 1.9 yr. 
r\! 
m 
> Yim (8, (1) 
l=2,m a 


a being the nuclear radius, ¥;, a real spherical harmonic 
normalized to unity, while r’ represents the vector to a 
point in the distorted nucleus whose position in the undis- 
torted nucleus is given by r.? Since in an excited state, the 


a's will be given by ain =aim sinwt, the total energy of 
these vibrations is 
3a* 2 21+ 


2 
v 
A MX l ’ 


E= 


AM being the mass of the nucleus.* 
We make use of the classical expression for the rate of 
radiation from such a system of currents 


1 
where 


j= jo(r) exp(iwr-n/c)dr, 


provided the current density at rf is given by jo(r) coset. 
In the evaluation of this integral for j, it is convenient to 
use the expansion 


kr 


the amplitude of the current density being obtained from 
(1). For the /th mode of vibration, we have 
(0)? 

9Z7e2 1+ 1 

Se 
where it has been assumed that all the a’s belonging to a 
given / and different m values are equal. In obtaining this 
result, the integral for j has been evaluated by retaining 
only the first term in the power series expansion of the 
Bessel function.‘ 

For the mean lifetime 7: of a level with given / value, 
we have 


1 32Z%? 1+1 (aw/c)*! 


71 4 


We have calculated the lifetimes in Table I from this 
formula for the case of Hg, using a radius a= $A le?/me?. 
In conclusion, I should like to acknowledge my indebted- 
ness to my colleague E. Feenberg for helpful discussions. 


1H. A. Bethe, Rev. Mod. Phys. 9, 226 (1937). 

2? Terms in/ =0, 1 are absent because of the conditions of incompres- 
sibility and fixed center of gravity. 

3Strictly, the frequencies occurring here are determined by the 
surface and Coulomb energies, whereas we regard them as poe. 

4 In the derivation of Bethe mentioned above, an order of magnitude 
estimate of the integral is used. 
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Further Evidence for a Single Component in the 
Primary Radiation 
W. F. G. Swann 


Bartol Research Print tt of The Franklin Institute, 
Swarthmore, Pennsylvania 


April 29, 1941 


N former communications,'~“ the writer has given 
reasons to suppose that there is no electron component 
to the primary cosmic radiation, that the mesotrons are 
secondary to primaries, and that the only electrons in the 
atmosphere are secondary to the mesotrons.’ He has 
pointed out, moreover, that insofar as mean life con- 
siderations are more potent than energy loss in deter- 
mining variation of mesotron intensity with altitude, such 
variation in intensity should follow approximately an 
exponential law with a single coefficient when plotted 
against true distance. When intensity is plotted against 
water equivalent distance, the curve becomes compressed 
parallel to the distance axis in such a manner that the 
result is an apparent increase of absorption coefficient with 
altitude.’ 

The purpose of the present communication is to add 
further evidence in support of the foregoing views. Our 
interest concerns that region of the atmosphere which 
lies below the altitudes at which the approach to a maxi- 
mum is felt.® 

Figure 1 gives measurements made by W. F. G. Swann 
and W. E. Danforth.’ In curve A, logio J/Jo, the logarithm 
of the ratio of intensity to that at sea level is plotted against 
true altitude in kilometers. In curve B, logio I/Jo is plotted 
against water equivalent altitude. The constancy of the 
slope of A bears out the conclusions cited above, and the 
change of slope of B shows how, on a water equivalent 
basis, the existence of an increase of absorption coefficient 
with altitude, and a consequent existence of more than one 
component, is simulated. 

If x represents true altitude in kilometers, and 
I=Ioexp[—ux] we obtain, from the slope of A, 1=0.156 
X 2.3 =0.36. Now if v is the velocity of the mesotrons, and 
t the time to travel the distance x, we may write ux = yt 
=0.36X3 X 105 = 105 where we have assumed v equal to 
the velocity of light for purely kinematical purposes. Thus 
exp[— 10°]. 


4 1 4 
log, L log, 
10 
A B 
a 
os 
meters of 
Kilometers 


Fig. 1. A—intensity variation vs. true altitude; B—intensity variation 
vs. water equivalent altitude. 
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If the rays disappeared simply by mean life considera- 
tions, we should have J=J,)exp[—t/r] where + is the 
mean life. A comparison of these two expressions gives 
7=10- sec. If the mean life for a stationary mesotron is 
10-6 sec., this tells us that the mesotrons would have to 
have a mass 10 times the rest mass, which would give them 
an energy of 10° ev. If we should assume Yukawa’s value 
5X10? sec. for the mean life for a mesotron at rest, the 
corresponding energy would be 210° ev. This value is 
sufficiently near the probable average value of the energy 
of the mesotrons to provide support to the foregoing views. 
A distribution of energy introduces a distribution of 
apparent absorption coefficients, but high energies give 
small absorption coefficients, and the story provided by 
energy distribution for high energies concerns increase of . 
hardness for great depths, where, of course, disappearance 
of rays by energy absorption becomes more important. 

1W. F. G. Swann, Phys. Rev. 56, 209 (1939). 

2 W. F. G. Swann, Rev. Mod. Phys. 11, 242 (1939), (see in particular 
7 IWR Swann, Phys. Rev. 58, 200 (1940). 

4W. F. G. Swann, Phys. Rev. 59, 770 (1941). 

> Further evidence to support these views has also been supplied by 
M. Schein, W. P. Jesse and E. O. Wollan, Phys. Rev. 59, 615 (1941). 

® The maximum becomes accounted for on the old and new views, 


but different fashion. 
F. G. Swann and W. E. Danforth, J. Frank. Inst. 228, 43 (1939). 


The Stable Isotopes of Nickel 


G. E. VALLEY* 


Research Laboratory of Physics, Harvard University, 
ambridge, Massachusetts 


April 23, 1941 


HE relative abundances of the nickel isotopes have 
been measured by a large number of workers.'!~> 

In these results discrepancies exist either between the 
various values of the ratio Ni®™ : Ni™, or between the 
atomic weight calculated therefrom and the chemical value. 
The present results were obtained in August 1940, with 
the type of mass spectrometer in which an initially mono- 
energetic ion beam is passed through a semicircular 
magnetic momentum analyzer, the emergent ion current 
being proportional to the isotopic abundance. Ions were 
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formed by electron bombardment of a beam of Ni vapor 
obtained from a hot, nickel-plated platinum wire. The 
nickel was especially purified for the experiment. 

Figure 1 shows a typical mass spectrum of nickel and 
serves to demonstrate that no corrections were necessary 
due to incomplete resolution. The values obtained are: 

Ni 58 60 61 62 64 

67.4 26.7 1.2 3.8 0.88 
From the nickel masses of Okuda et al. and the oxygen 
ratio 1.000275,’ the value 58.71+0.02 is obtained for 
the chemical atomic weight. This agrees within the 
experimental error with the chemical value 58.69+0.01.° 
The ratio Ni™ : Ni® obtained here is in agreement with 
the recently published value of Straus,’ and in approximate 
agreement with the photographic results of Dempster.‘ 

The assistance of Dr. H. H. Anderson who prepared 
and purified the specimen is gratefully acknowledged. 
The advice of Professor K. T. Bainbridge was most helpful. 

* At present National Research Fellow at ecward University. 

1F. W. Aston, Proc. Roy. Soc. 149, 35). 

2 J. de Gier and P. aw Proc. K. Akad. Amst. 38.8, 810 (1935). 

ow. A. Lub, Proc. K. Akad. Amst. 423, 253 (1939). 

A. J. Dempster, Phys. Rev. 50, 98 (1936). 

: H. A. Straus, Phys. Rev. 59, 430 (1941). 

6 T. Okuda, K. Ogata, H. Kuroda, S. Shima and S. Shindo, Phys. 
Rev. 59, 104 (1941). 

7W. R. Smythe, Phys: Rev. 45, 299 (1934). 

8 B. F. Murphey, Phys. Rev. 59, 320 (1941). 


*G. P. Baxter, M. Guichard, O. Hénigschmid and R. Whytlaw-Gray, 
J. Am. Chem. Soc. 62, 669 (1940). 


Radiative Collision of Neutrons with Protons 


TOSHINOSUKE Muto 


Institule of Physics, Kyushu Imperial University, Fukuoka, Japan and 
The Institute for Advanced Study, Princeton, New Jersey 


April 28, 1941 
S is well known, capture with radiative emission and 
elastic scattering are the principal collision processes 
of neutrons with protons. Although the capture process is 
known to be predominant only in the narrow range of very 
slow neutron energies, elastic scattering will probably 
play an important role over a wide range in the primary 
neutron energy owing to the larger contribution to the 
cross section of the partial waves with relatively high 
angular momenta. For very high neutron energies, how- 
ever, we must consider the possibility of radiative emission 
in the collision of neutrons with protons, since the protons 
set in motion during the collision are capable of emitting 
radiation. In order to find the order of magnitude of the 
cross section, we have worked out such a radiative col- 
lision of heavy particles. The process is similar to that of 
the radiative scattering (bremsstrahlung) of charged par- 
ticles in an atomic field. The calculation, therefore, has 
been treated in a similar manner to that of the brems- 
strahlung of charged particles. 

For the nuclear interaction between neutron and proton 
we have assumed, as usual, a static potential with an 

exchange operator of short range, 


V=—Vo(e*"/r)O, (1) 
in which O is the exchange operator (Heisenberg’s or 


Majorana’s), r the mutual distance, \= mc/h the reciprocal 
Compton wave-length of a meson, and Vo the constant 
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dependent upon the depth of the potential. As a prelimi- 
nary, the motion of heavy particles has been treated non- 
relativistically,! the radiative emission induced by the 
magnetic dipole moment of a neutron being neglected, and, 
further, we have used Born’s approximation. 

For the case that the neutrons with a given momentum 
py® and an energy Ey® impinge upon protons initially at 
rest, we have obtained the following expression for the 
differential cross section: 


1 e? dhy sin’@ 


x pot)? (2) 


in which hvy=the emitted quantum, M is the mass of a 
heavy particle, cos? =(k,-pp)/kypp, cos6’ = (k,: 
-hk, and pp represent the momenta of the photon and 
recoil proton, respectively. dQ, and d&p are their elements 
of solid angle. According to the conservation equations of 
energy and momentum, pp is found to depend upon @, ¢ 
and 6’ if the primary energies of the colliding neutron and 
of the emitted quantum are known. For the special case 
that hy/Mc*<1 and py®/ Mc <1, the bracket expression in 
(2) may be shown to depend on the angles through the 
function pp only. 

For an estimate of the order of magnitude of the total 
cross section ¢, for the emission of a quantum between hy 
and hvy+dhy, it is sufficient, in the above approximation, 
to consider the average value [2M(Ey®—hv)]}-+-g of the 
reduced bracket expression in d¢, taking into account that 
pr varies from 0 to (py°?—2Mhyv)}. 

Then, integrating d@ over all angles, we get 


Vo\? 1 dhv g(Ey®, hv, me*) 
2M hv Ey® \ 


where g is a rather lengthy expression dependent upon 
Ey®, hv and mc*®, which reduces to (Ey®—hyv)/(m*c?/2M) 
for Ey’>>hyv. 

For neutrons of 10° electron volts energy, the above 
formula gives a cross section of less than 10-” cm? for the 
emission of y-rays of 10° electron volts. The corresponding 
cross section for the bremsstrahlung of an electron in the 
Coulomb field of a proton is found to be less than 10-* cm?, 
being very large compared with that of the former process. 
This is to be expected owing to the very short range of the 
nuclear force and the relatively large mass of a proton. 
Since the cross section for elastic scattering becomes* 
approximately 10~* to 10~*’ cm? for these neutron energies, 
the radiative collision of heavy particles may be reasonably 
neglected in the study of neutron scattering by protons. 
In connection with cosmic-ray phenomena, however, it 
seems desirable to calculate the cross section of this process 
for extremely high neutron energies.! 


i- (3) 


! The relativistic calculation, based upon the present meson theory, 
of the radiative collision of heavy particles is now in progress; a full 
a) will be published elsewhere, together with details of this work. 

A. Bethe and R. F. Bacher, Rev. Mod. Phys. 8, 122 (1936). 
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Comment on the Elastic Constants 
of Alpha-Quartz 


A. W. Lawson 


Randal Morgan Laboratory of Physics, University of Pennsylvania, 
Philadelphia, Pennsylvania 


March 13, 1941 


TANASOFF and Hart! have recently described some 
admirable experiments on the vibrations of variously 
oriented crystal quartz plates. Their results, however, 
for the elastic constants of quartz, calculated on the 
assumption that the plates were infinite, evince a small 
but unexplained discrepancy with the classical elastic 
theory for substances possessing the symmetry of quartz, 
namely crystal class D3. 

Now Cady? has shown that in the calculation of the 
eigenfrequencies of a piezoelectric crystal, the forces 
arising from the polarization produced by the vibration 
must be taken into consideration. This has not been done 
by Atanasoff and Hart, and the object of this letter is to 
show that when cognizance is taken of the converse 
piezoelectric effect, the discrepancy between elastic theory 
and experiment is removed. 

If, in an anisotropic substance, u is the displacement 
vector, @ the strain tensor, ¢ the stress tensor, ¢ the elastic 
tensor, F the volume rate of force developed by the stress 
tensor, E the electric field intensity, P the polarization, 
and e the piezoelectric tensor, we have: 


Ou; 
¥ij Cija89a3— ~ CijyEy, (2) 
(4) 
ap 


where &; is a principal electric susceptibility. Now by 
symmetry, for an infinite plate oriented perpendicular to 
the x; direction, we have: 


9a8 = (5) 
‘js (6) 


If cognizance is taken of Eqs. (5) and (6), the second of 
Maxwell’s equations yields: 
where K;, is a principal dielectric constant. 
The first and third of Maxwell's equations yield, to 
within the precision of the calculation: 


@E,/dx?=0 for (8) 


The foregoing relation is a consequence of the fact that 
the wave-length of the resulting electric wave is very long 
compared with the wave-length of the elastic wave. 

We now assume that the external field is zero. The 
differential equations (3), (7) and (8) must be solved 
subject to the boundary conditions: 
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a Y 
K,E,+= CaiyIay = 0. (9) 


This requires that £, =0 for i#~+. Our differential equation 
for the propagation of elastic waves then becomes: 


2 
= { |e. (10) 


If we assume 


uj= Aj cos Mi piwe 


the condition that Eq. (10) be satisfied is the secular 
system: 


Ag—wA;=0, (11) 
a 
where 
* 
Cijait 
and 
f=w/2r, (12) 


where s is the thickness of the plate and f is the observed 
resonant frequency. 

The use of Atanasoff and Hart’s data for the resonant 
frequencies of vibration of variously oriented plates, and 
the assumption that quartz is trigonal yields the elastic 
constants of quartz and several internal checks. These 
constants are shown in Table I and compared with those 
computed from Atanasoff and Hart’s incomplete theory 
and those observed statically by Voigt.* 

The internal consistency of the data verifies that 
quartz is trigonal. This may best be shown by calculating 
some resonant frequencies and comparing these with 
observed frequencies not used in the calculation of the 
elastic constants. Thus for a (30°, —38°12}’) plate: 


« calc. = 64.86 X 10"° dynes/cm?, 
« obs. = 64.92 x 10! dynes/cm?, 


and for a (30°, 0°) plate: 


« calc. = 49.36 X 10" dynes/cm*, 
« obs. = 49.39 X 10" dynes/cm?*. 


The differences between the observed and computed 
values are within the experimental errors of the quantities 
involved. 

The value used for e,; was 5.29X10* e.s.u. units as 
reported by Van Dyke,‘ and the values for K; were taken 
directly from Cady’s paper. 

One other point should be noted. Atanasoff and Hart 
used an air gap between the crystals and the plates 


TABLE I. Elastic constants of alpha-quartz. The units are 10 dynes/cm*. 


ATANASOFF ATANASOFF 
AND HART AND HART 
VoIGct UNCORRECTED CORRECTED 
cu 86.8 87.55 86.75 
7.1 6.07 6.87 
ca 58.2 57.19 57.86 
c14 17.2 17.25 17.96 
C13 14.4 13.3 11.3 
€33 107.5 106.8 106.8 


(9) 


ion 


10) 


lar 


11) 


BS 
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exciting the vibrations. Thus the external field is not zero, 
and this produces another shift in the fundamental 
frequency, of the same order of magnitude but opposite 
to that contemplated here. However, Cady has shown 
that this correction decreases as 1/n* and since m is very 
large in the experiments under review, this correction is 


negligible. 


1J. V. Atanasoff and P. J. Hart, Phys. Rev. 59, 85 (1941). 

2 W. G. Cady, Physics 7, 237 (1936). 
3W. Voigt, Lehrbuch der Kristall physik (Teubner, Leipzig, 1928). 
4K. S. Van Dyke, Phys. Rev. 42, 587 (1932). 


Note on the Structure of Liquid Helium 


L. I. Scuirr 


Randal Morgan Laboratory of Physics, University of Pennsylvania, 
Philadel phia, Pennsylvania 


April 15, 1941 


1. The existence of surface transfer films of liquid 
helium II of the order of 100 atoms in thickness! has 
sometimes seemed difficult to understand, although 
actually this thickness is in accord with the usual theory 
of interatomic forces. The difference between the 1/r° 
van der Waals interactions of two helium atoms and of a 
helium atom and an atom of the solid container yields a 
1/a* dependence on the distance from the surface. In the 
case of a metallic container, there is also a comparable 
1/a* energy arising from interaction with the electrons of 
the metal.? If a solid surface projects out of a reservoir 
of liquid helium II, the thickness a, of the film at a height 
h cm above the free surface of the liquid is obtained by 
equating the total interaction energy to the gravitational 
potential energy of a helium atom. This gives for copper: 
an=4.3X 10~°/h' cm, which is in excellent agreement with 
the rather uncertain experimental value of 3.5 cm 
for a sheet extending a few cm above the liquid surface.' 
The numerical value 4.3 for copper is replaced by 4.7 for 
silver, ~4 for glass, and 2.2 for rocksalt (100 plane ex- 
posed). Although such films probably exist in general, 
this simple calculation is valid only when the fluid under 
consideration is sufficiently nonviscous so that it can come 
into equilibrium under these forces before it evaporates; 
this explains the failure to observe such films with helium 
I. Accurate measurements of film thickness would give 
useful information about long range interatomic forces. 

2. The phase transition in liquid helium appears to be a 
continuous one in the sense that all properties of the state 
that do not involve differentiation with respect to the 
temperature are continuous.’ It is therefore plausible to 
assume that helium I consists of a single phase A and 
helium II of a mixture of two phases A and B, phase B 
being absent at the A-point and growing continuously at 
the expense of phase A as the temperature decreases. 
This point of view was proposed by Tisza* and developed 
by H. London ;' they showed that the experiments could 
be understood, qualitatively at least, by assuming that A 
is disordered and has a normal viscosity, while B is ordered 
and has little or no viscosity. However, there appears to 
be no need to assume, as has sometimes been done,®* 
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that the liquid in the surface layer is qualitatively different 
from that in the volume of the fluid, since the capillary 
experiments can be understood equally well in terms of 
the viscous drag exerted by a nearby wall on phase A. 
Moreover, a different structure of the surface layer could 
be expected only if the interaction of an atom with the 
wall were comparable with the interaction with its neigh- 
bors. The calculations of paragraph 1 indicate that this 
is true only for the first 2 or 3 atom layers, while the 
thermomechanical effect is appreciable for much larger 
channels.” 

3. The phase .1 referred to in paragraph 2 can be 
understood in terms of an Einstein-Bose gas model, in 
which each helium atom moves in the periodic field that 
approximately represents the effect of correlations between 
the positions of the atoms. This model explains the x-ray 
data, the position of the A-point® and the positive temper- 
ature coefficient of viscosity of helium I.*-* Its application 
to phase B does not, however, explain the observed 
specific heat anomaly,® and the failure of a more complete 
treatment!® to accomplish this makes it seem likely that 
no gas model will succeed here. It is natural then to 
assume a more typically liquid model for phase B, in 
which there exists a discrete energy level for each atom. 
Because of the weak attractive forces and large zero-point 
kinetic energy of helium atoms there will not be more 
than one such state per atom, if any, and it would probably 
lie ~10-* ev below the continuum of phase A states. 
These “bound” atoms would diffract over several atom 
distances, so that the aggregation would not possess the 
rigidity of a lattice. The energy of a given bound atom 
could depend on the number of its neighbors that are 
similarly bound; the inclusion of an “ordering energy” of 
this type would favor the transition from A to B and 
sharpen the specific heat curve. Further, volume-conserving 
deformations would not alter the energy of phase B, and 
the adiabatic nature of laboratory manipulations as 
compared with the period of the system (~10~ sec. for 
B states depressed by 10° ev) would make excitation to 
continuum states extremely improbable. There would thus 
be no resistance to deformation, no mechanism for dissi- 
pating flow energy, and hence no viscosity. While this 
model has some physical plausibility to recommend it, 
it must be regarded as tentative because of the great 
difficulty of calculating the several parameters that appear 
in it in terms of atomic properties. 

The author takes pleasure in thanking Professor F. 
Seitz for interesting discussion of these ideas. 

' J. G. Daunt and K. Mendelssohn, Proc. Roy. Soc. A170, 429 (1939). 

2 J. Bardeen, Phys. Rev. 58, 727 (1940). 

3W. H. Keesom and G. E. MacWood, Physica 5, 737 (1938), have 
obtained data that indicate a possible discontinuity in the viscosity 
at the A-point. The improbability of the coexistence of distinct thermo- 
dynamic states having the same energy and free energy makes a further 
experimental investigation of this point desirable. 

4L. Tisza, Nature 141, 913 (1938). 

5H. London, Proc. Roy. Soc. yh 484 (1939). 

6 F. London, Phys. Rev. 54, 947 (1938). 

7J. F. Allen and H. Jones, * wh 141, 243 (1938); J. F. Allen and 
J. Reekie, Nature 144, 475 (1939). 

* Cf. reference 6. From the present point of view the A-temperature 
depends essentially on the properties of phase A, and may be chan 
somewhat by introducing an effective mass for the helium atoms due 


to their motion in the periodic field (cp. electronic effective mass in 


the theory of metals). 
*W. H. Keesom and A. P. Keesom, ia 2, 557 (1935). 
0 L. I. Schiff, Phys. Rev. 59, 758 (1941) 
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On the Production of Mesotrons 


J. FRANKLIN CARLSON, Department of Physics, Purdue University, 
Lafayette, Indiana 


AND 


MARCEL SCHEIN, Department of Physics, University of Chicago, 
Chicago, Illinois 


April 22, 1941 


HE recent letter of Schein, Jesse and Wollan'! makes 
it seem doubtful that the generally accepted ex- 
planation of the atmospheric transition curve of the 
cosmic rays is valid. If, as the evidence suggests, the 
primary rays entering the atmosphere are largely protons 
and not electrons, one must explain the production of 
both the penetrating component, the mesotrons, and the 
soft component, electrons and y-rays, and their intensity 
curves in the atmosphere. An obvious assumption is that 
the protons create mesotrons in collisions with nuclear 
particles. The mesotrons in turn decay, transferring a 
part at least of their energy into an electronic component 
which gives rise to the soft component. 

Since the nature of the mechanism of these processes is 
not known, we are suggesting a simplified quantitative 
model of the process based on a few reasonable assumptions. 

The main assumptions we make are (a) the creation is 
an “explosion” in which the proton loses all of its energy 
and creates m mesotrons; (b) the mesotrons have roughly 
equal energy, and (c) the multiplicity m, the average 
energy of the mesotrons and the production cross section 
increases with the proton energy over the region from 
210° to 2X10 ev. The multiplicity and the cross 
section remain constant for higher energies. 

If we assume for protons of 10! ev a cross section of 
1.6 10-*° cm? per nuclear particle, and a multiplicity of 
ten, we find from the solution of the diffusion equation 
for the mesotron intensity as a function of the pressure, 
a curve which starts out at zero for zero pressure and 
rises sharply. In the first 1/12 of the atmosphere it reaches 
a broad maximum of about twice the initial proton count. 
It then drops off with a slight positive curvature. The 
total intensity of the protons and mesotrons has a slightly 
higher maximum and drops off more rapidly. Since the 
average mesotron intensity never greatly exceeds two, it 
follows that most of the mesotrons produced decay near 
the top of the atmosphere. 

Because of the extreme rarity of the atmosphere near 
the top, the time required to traverse a small fraction of 
a homogeneous atmosphere is ample to allow disintegra- 
tion. Thus the average number of mesotrons which have 
not decayed is very small compared to the multiplicity n. 
The greater part of the mesotrons which decay in the 
upper atmosphere will have lost little energy by ionization. 
There would, therefore, be between 8 and 9 electrons in 
the upper atmosphere per incident proton. These disinte- 
gration electrons with energies of about 5X 10® ev should 
give rise to some cascade processes and be responsible for 
the soft component. 

High energy mesotrons should produce a high east-west 
asymmetry since for high energies the angular divergence 
of the mesotrons produced should be small. At lower 
energies there should be some angular spread and conse- 


quently less asymmetry. Comparison of the difference 
curves for the mesotron intensities at Chicago, latitude 
51°, and Texas, latitude 38°, shows that the greater part 
of the mesotrons must be produced by protons of energy 
above 6X 10° ev. We would thus expect the cross section 
to be larger and the multiplicity greater for higher energy 
protons. 

Near the equator where the mean energy of the protons 
is higher, if the cross section increases with proton energy 
and the mean energy of the mesotrons increases, one finds 
a larger number of mesotrons at the maximum due to the 
lower rate of decay. For the extremely high energy protons 
which are responsible for the mesotrons below sea level, 
this mesotron intensity should follow very closely the 
energy spectrum of the protons if the multiplicity ap- 
proaches a fixed value. Comparison of these results with 
experiment is at least qualitatively satisfactory. 

Professor R. Serber? has kindly pointed out to us the 
following facts. From the area of Pfotzer’s curve the 
multiplication ratio (intensity at maximum /intensity at 
top) is about 10.7 whereas experimentally it is about 5. 
This suggests that only half of the energy appears in 
ionization, the remainder presumably going into neutrinos. 
Also the same conclusions follow from the absolute rate 
at the top determined from absolute sea level count times 
the experimental ratio of sea level to top, when compared 
with Bowen, Millikan and Neher’s calculation of the 
incoming numbers. From the cascade theory the multi- 
plication should be about 13.5. The comparison of the 
multiplication ratios shows that the energy gets into the 
shower producing component quickly. 

These results are only preliminary and suggestive. A 
more careful analysis is in progress which will be reported 
at an early date. 


1 Schein, Jesse and Wollan, Phys. Rev. 59, 615 (1941). 
2 R. Serber, private communication. 


Proton-Induced Fission 


G. DESSAVER AND E. M. HAFNER 
University of Rochester, Rochester, New York 
April 2, 1941 


UCLEAR fission induced by the capture of a charged 
projectile has been established by the work of 
Jacobsen and Lassen! after some preliminary experiments 
of Gant.? These authors used deuterons of energies between 
8 and 9.5 Mev on thorium and uranium targets. For 
6.0-Mev protons on uranium, Bohr and Wheeler® predict 
a value of 10-*8 for the fission cross section, which should 
be observable. 

We have bombarded thorium metal and uranium oxide 
with the 6.9-Mev proton beam of the Rochester cyclotron. 
Our previous observations with uranium targets had led 
us to the conclusion that a (p, m) reaction probably also 
takes place. For this reason, care had to be taken that 
fission caused by secondary neutrons from the target was 
not interpreted as being proton induced. 

The proton beam entered an evacuated cylindrical 
brass chamber (Fig. 1A) in which targets could be placed 
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NEUTRON FISSION 


BRASS - THORIUM METAL MIL 
LUCITE URANIUM OXIDE 


in any desired position. In a side tube (which was designed 
to accommodate a proportional counter) projecting at 
right angles from the body of the chamber, and well 
removed from the beam, foils were mounted so as to catch 
fission fragments emitted from the target surface. A 
0.1-mil aluminum foil in front of the catching foils was 
found sufficient to retain all ordinary radioactive recoil 
atoms. Aluminum foils placed behind the catchers served 
to correct for radioactivities induced in them by neutrons 
or scattered protons. 

A 10-mil sheet of thorium, which was found to be thick 
enough to stop the proton beam and _ proton-induced 
fission particles, was oriented so that its normal made 
angles of 45 degrees with both the beam and the catchers’ 
normal. There are two such positions. In position I 
(Fig. 1B) the irradiated surface of the target faces the 
catcher; in position II, the other surface does. The amount 
of fission fragments caught in position II will be approxi- 
mately the same as that fraction caught in position I that 
is due to fission induced by the neutrons from the cyclotron 
and the target. After a bombardment of 5 Ah in 1.5 
hours, radioactivity was observed in a 1-mil aluminum 
catcher, whose decay could be followed for seven hours. 
A bombardment of 25 wAh in 8 hours yielded a sample 
which could be observed for over two and a half days. 
Identical bombardments with the target in position II 
gave no effect. When plotted on semi-logarithmic paper, 
the decay curves show the curvature characteristic of a 
mixture of many periods. 

By replacing the 1-mil aluminum catcher with a stack 
of 0.1-mil aluminum foils it was found that the maximum 
range of fission fragments from 9Pa**, the compound 
nucleus in question, lies between 0.5 and 0.6 mil of alumi- 
num. Reducing the proton energy in successive steps by 
inserting aluminum absorbers just in front of the cyclotron 
window gave several points of the excitation function. 
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The lowest proton energy for which we obtained an 
observable effect was 5.8 Mev. Had our thorium sheet 
been large enough to insure that no fraction of the proton 
beam would fail to strike it, we could have collected the 
background-induced fission fragments behind the target, 
without turning it into position II and repeating the 
bombardment. 

In the case of uranium, however, we proceeded in this 
more direct manner. Figure 1A shows the position of a 
second catcher behind the target. It also shows the elliptic 
perforated brass plate which was packed and coated with 
uranium oxide embedded in pyroxylin. The brass provided 
sufficient cooling, and the oxide-filled channels gave the 
surface coats a higher stability. Thus we were able to find 
in a single run that uranium gives an effect (fission of 
93°) of the same order of magnitude as that found in 
thorium, and that the background effect is negligible. 

Chemical investigations and experiments employing an 
ionization chamber for more accurate yield and range 
determinations are under way. 


iI. and N. O. Lassen, Rev. 58, 867 (1940). 
2D. H. T. Gant, Nature 144, 707 (1939) 
2 Niels Bohr and J. A. Wheeler, Phys. Rev. 56, 426 (1939). 


Note on the Normalization of Dirac Functions 


C. C. Lin 
Department of Applied Mathematics, University of Toronto, Canada 
April 18, 1941 


VERY simple method of normalizing the Dirac 
functions in the Kepler problem is given in this note. 
The normalization of Dirac functions in the Kepler problem 
was first done by Bechert' with the help of certain contour 
integrals in a complex plane. The calculation was, however, 
rather complicated. With the functions given in Infeld’s 
form, a much simpler method of normalization leading to 
a neater result is proposed here. The result is, as I verified, 
the same as that given by Bechert, but it appears in a 
completely different form. 
It has been shown? that the radial Dirac functions are 
given by 


with the following definitions and formulae (C being the 
normalization factor) : 


yi=(k—aZ)', (3) 


where / is a positive integer, or zero (only for «<0), 
and « is an integer, positive or negative, but not zero; 


0 
Fi and §;,, are the lowest two functions of a ladder of 
normalized functions defined by 
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1 


(l+y)(m+v7) 


m+y 


x! +Ye~ az/(t+7) (5) 


(x) = 


where m is an integer, 0< m< 1, and 


a= ead. (7) 


As C is the normalization factor and the functions 
o 
1, Fi. are normalized, we obtain from (1), (2) and (3) 


the condition for C to make (x2+x2*)dx=1: 


where 
oD 0 1 
1,0, 1)= J, Fi, 4 (x) Fr, y(x)dx. (9) 
It will be shown below that 
I,(0, 1) = — (ex? —y")'/aZe, (10) 
and hence, because of (8), we have 


C= e/2v?2. (11) 


This is the normalization factor required. 
To calculate J,(0, 1), let us consider 


m+1)= dx. (12) 
By writing 


a 


it is not difficult to verify that 
HO = + + (14) 


where Em, Mm, fm are constants given by 


(13) 


Ent am = m+7y—1 1’ Em =1, 


1 1 
(15) 


Furthermore, we have 


provided f and ¢ vanish properly at x=0 and x= ~. 
By making use of (13)-(16), we can easily see that (12) 
reduces to 


m+ 1) (Ent 1, m)+inh, (17) 


or 
Ii(m—1, m) a Ii(m, m+1) a 


We note that this means that the quantity on either side 
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of this equation is independent of m, and hence is equal to 
a/(l+y)?, (19) 


since J,(/,/+1) is evidently zero (as =0). Hence 
by using the value of A”*! given by (13), we have the 
value of J;(m, m+1). Putting m=0 in the result, we have 


1)= — (20) 


Equation (10) is then obtained by solving (4) for /-+y and 
substituting the result into (20). 

The author wishes to thank Professor Infeld for suggest- 
ing the problem and for his help. 

1 Bechert, Ann. d. Physik 4,7 me (1930). See also Bethe, Handbuch 


der Physik (1933), Vol. 24, 
2 L. Infeld, Phys. Rev. 59, P537 (1941). The notations used here are the 


same as those used in that paper. 


Note on the “Kepler Problem” in a Spherical 
Space, and the Factorization Method of 
Solving Eigenvalue Problems 


A. F. STEVENSON 


Department of : Mathematics, University of Toronto, 
Toronto, Canada 


April 30, 1941 


CHRODINGER! has developed an elegant ‘“‘factor- 
ization” method of solving certain eigenvalue prob- 
lems, an improved version of which has recently been 
given by Infeld.2 The only problem treated by these 
authors which had not previously been considered and 
solved by other more conventional methods is the ‘‘ Kepler 
problem” in a spherical space. As Schrédinger stated that 
he found this problem “‘difficult to tackle in any other 
way,” it may perhaps be of interest to indicate briefly 
how the solution may be obtained without too great 
complication by conventional methods. The differential 
equation can, in fact, easily be transformed into a standard 
type, but the nature of the singularities of this transformed 
equation makes the discussion a little different from usual 
(explicit use is made of the continuity, as well as the 
boundedness, of the solution). It may also be opportune 
to offer a few remarks on the applicability of the factor- 
ization method in general. 
The problem referred to leads to the equation ((4.3) of 
reference 1, or (7.1) of reference 2) 


7) +0 sin?y +2u sinx cosx—/(/+1)]y=0, (1) 
where X is the eigenvalue parameter, u a given constant, 
and /=0, 1, 2, ---.3 A solution which is bounded and 
continuous in the interval 0<x<_m is required. The 
substitution x =cotx transforms (1) into 


A+2ux 
and the fundamental interval is now —*x<€x¢+. 


Equation (2) has regular singularities at the points 
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x=i, —i, *, and since these are the only singularities, 
it may be solved in terms of hypergeometric functions.‘ 
A solution, valid for large |x|, which remains bounded 
at «, is 


y=(x—1) 
1+1+a—a*, 2/+2, —2i/(x-—i)), (3) 


where F denotes the hypergeometric function, @ is a root 


of the equation 
a®—a=(\—2ip)/4 (4) 


and a* is the conjugate complex of a. The second solution 
behaves like |x|'*! at «, and must be discarded. 

The series occurring in (3) is convergent if x lies within 
the circle C, defined by |x—i| =2. Because of the singu- 
larities at +7, however, the analytic continuation of (3) 
along the real axis inside C is not in general a continuous 
function. A sufficient condition for the existence of an 
acceptable solution is that the series in (3) should termi- 
nate. This requires the relations 


lta+ta*=—n’, n’=0,1,2,--- 


or, from (4), 
A=n?—1—p2/n?, n=n'+/1+1 


in agreement with the eigenvalues found by Schridinger 
and Infeld. To show that the condition is also necessary, 
we must examine the explicit form of the analytic con- 
tinuation of (3) inside C. This may readily be done, since 
relations between the various solutions of the hyper- 
geometric equation are completely known.‘ We thus find 
that the solution defined by (3) is, indeed, discontinuous 
at x=0 along the real axis, unless the above condition is 
satisfied. 
The (un-normalized) eigenfunctions may be written 


=sin'xe F(— nn’, 2/42, 1—e?**). 


These eigenfunctions are real in spite of their apparent 
complex form as may be shown by utilizing known rela- 
tions of the hypergeometric function. 

The general theory underlying the factorization method 
has been investigated by Coleman.® Although the analysis 
is not quite complete, it is very probable that the method 
is, in practice, restricted to certain equations of hyper- 
geometric or confluent hypergeometric® type, or equations 
reducible to these by simple substitutions. But this does 
not, of course, detract from the elegance of the method 
in those cases where it is applicable (which include most 
of the soluble problems of quantum mechanics). It is, 
indeed, remarkable that eigenvalue problems associated 
with equations whose complete analysis is fairly compli- 
cated, can be solved in such a simple manner. 

I wish to thank Professor Infeld for letting me see his 
paper in advance of publication, and for interesting 
discussion. I am also indebted to Mr. Coleman for com- 
municating to me his results before publication. 

' E. Schrédinger, Proc. Roy. Irish Acad. 46A, 9 (1940). 

2 L. Infeld, Phys. Rev. 59, 737 (1941). 

3’ The results hold, however, for any 1 >0. 

‘See, for instance, E. T. Whittaker, and G. N. Watson, Modern 
Analysis (Cambridge, 1927), Chapter 14. 


5A. J. Coleman, not yet published. 
Reference 4, Chapter 16. 
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Distinction between Longitudinal and 
Transverse Mesons 
M. Kosayasi 
Physical Institute, Faculty of Science, Osaka Imperial University, 
Osaka, Japan 
April 11, 1941 
OSMIC-RAY evidence shows that the cross section 
for the scattering of mesons by nucleons is much 
smaller than that calculated from the meson theory of 
nuclear force fields, while the interaction of mesons of 
unit spin with photons of high energy seems also to be 
too strong to account for the observed soft secondaries of 
the hard component. We want to propose a way of re- 
moving these discrepancies in the vector meson theory. 
If one adopts the single force hypothesis proposed by 
Bethe! and leaves out the contribution from the processes 
in which Heisenberg’s condition® 


(u being the mass of a meson, @ a number of order unity) 
for the applicability of the present quantum theory is not 
satisfied, the cross section for the scattering of the longi- 
tudinally polarized mesons by nucleons is at most of the 
order of 10-8 cm?, while the cross section for the trans- 
versely polarized mesons of high energy is roughly 


X 10°27 (2) 


Moreover, the transition probability for mesons from the 
longitudinal to the transverse state, which is mostly due 
to the process in which charged mesons are scattered in 
the electrostatic field of atoms, is very small, the cross 
section being given by 


2x (e?/uc*)*Z? log(137)?~Z? x cm’. (3) 


The above results lead to the conclusion that transverse 
mesons, if created, are largely scattered in the upper 
atmosphere, and that the greater part of the hard compo- 
nent found at sea level consists of longitudinal mesons. 
If we adopt the symmetrical theory of nuclear forces, the 
cross section for the transition of a transverse meson from 
the charged to the neutral state is also given by (2). 
This behavior of transverse mesons seems to favor the 
interpretation of the cloud-chamber experiments* made in 
the upper atmosphere. The very penetrating component 
which persists under thickness of matter might consist 
primarily of longitudinally polarized neutral mesons. 

We have also examined the radiative processes of 
longitudinal and transverse mesons separately. We have 
calculated the differential cross section for the collision of 
a longitudinal meson with a photon of equal and opposite 
momentum and found that its highest term in powers of 
the energy E of the incident photon is given by 


where @ is the scattering angle. This term has a very 
small value in the region where condition (1) is satisfied, 
whereas such a term as (1 —cos@)? does not appear in the 
corresponding expression for the transverse mesons or in 
the cross section for the creation of meson pairs by y-rays. 
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These circumstances seem to account for the small prob- 
ability of the emission of soft secondaries in mattert and 
the production of the hard component in the upper 
atmosphere. The cross section for a process in which a 
photon of very high energy is absorbed by a nucleon with 
subsequent emission of a meson is roughly given by 


ma? (ge /uc*)?~a? X 10°28 cm?, (5) 


and the meson thus created is the longitudinal one. By 
putting a3 in (5), we get the experimental value® for 
the cross section of the production of mesons by hard 
y-rays. In this connection, it is notable, as was shown by 
Sakata,® that the lifetime of the meson is also explained 
by making use of condition (1) and putting a3. 

Finally, we have calculated the interaction potential 


THE EDITOR 


between a proton and a neutron according to the single 
force hypothesis and the cutting-off condition (1), and 
found that for a<3 the force range is somewhat narrower 
than h/uc? and the spin-orbital coupling potential is more 
reduced than according to the cutting-off prescription in 
the coordinate space. We think, therefore, that the deu- 
teron problem will be solved in this way and the sign of 
the quadrupole moment may be explained by the precession 
of the nucleons. The detailed calculation of this problem 
is now in progress. 
H. A. Bethe, Phys. Rev. 57, 260 (1940). 

Heisenberg, Zeits. f. Physik 110, 2si 

3G. Herzog, Phys. Rev. 59, 117 (19 
Bostick, Phys. Rev. 59, 122 (1941). 

4M. Schein and P. S. Gill, Rev. Mod. Phys. 11, 267 (1939). 


5 M. Schein, W. P. Jesse and E. O. Wollan, Phys. Rev. 57, 874 (1940). 
6S. Sakata, Phys. Rev. 58, 576 (1940). 


. Herzog and W. H. 
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